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AB3TBACT 

!I!!he  basis  tor  current  propeller  design  methods  Is  lifting 
line  theory  supplemented  by  an  approximate  correction  for  lifting 
surface  effect.  Recent  studies  have  Indicated  that  this  correction 
Is  not  entirely  satisfactory,  amd  that  a  more  exact  lifting  surface 
theory  for  marine  propellers  Is  needed. 

In  the  present  work,  methods  are  developed  to  determine  pitch 
and  camber  corrections  for  propellers  vlth  arbitrary  blade  cutUne  and 
radial  load  distribution.  The  pitch  and  camber  Is  determined  by  the 
requirement  that  the  desired  load  distribution  be  obtained  vlth  the 
sections  operating  at  their  Ideal  angle  of  attack.  The  method  may  be 
used  both  for  homogeneous-flov  and  wake-adapted  propellers. 

The  method  Is  an  adaptation  of  the  vortex  lattice  method 
developed  for  vlngs  of  arbltreiry  shape  by  Ealkner.  By  replacing  the 
continuous  vortex  distribution  by  a  lattice  of  discrete  vortex  elements, 
the  singular  Integral  equation  occurring  In  lifting  surface  theory  Is 
replaced  by  a  set  of  linear  algebraic  equations. 

From  the  form  of  these  equations.  It  Is  shown  that  a  propeller 
with  symaetrlcgLl  blades  and  with  aiean  lines  idiich  cure  syametrlcal  about 
the  mid-chord  lias  no  pitch  correction  due  to  lifting  surface  effect. 

To  obtain  a  preliminary  check  on  the  accuracy  of  vortex  lattice 
theory,  methods  of  approximating  propeller  lifting  line  theory 
developed,  and.  numerical  results  obtained  with  an  IBM  JO9  Cosqputer  are 
given,  intese  results  agree  substantially  vlth  existing  lifting  line  data. 

Lifting  surface  results  obtained  with  an  IBM  7^9  aa  IBM  7090 
computer  are  discussed.  From  these  results  It  Is  tentatively  cosu:luded 
that  an  accuracy  of  +  2^  In  the  caid>er  correction  may  be  achieved  with 
reasonable  ccisputatlon  times.  The  saiqple  results  In^cate  that  lifting 
surface  corrections  are  dependent  on  such  variables  as  blade  shape  end 
circulation  distribution,  which  eure  act  taken  Into  account  In  current 
design  methods. 
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O 
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•  propeller  rotational  speed 


CHAPTER  I 


IMTROIUCTION 


Propeller  Design  Mettipdg^ 

Th®  basis  for  current  propeller  desigti  methods  is  lifting  Ifhe 
theoiy  supplemented  lay  aa  approximate  Correction  for  lifting  surface 
effect*  A  description  of  such  methods  may  he  fomil  in  recent  puhlloatlonc 
by  Iierbs^^^*  ,  Van  Maneil^^^^  and  Ecihardt  and  Mbrgan^^^#  Since 
the  historical  development  of  propeller  theory  ie  treated  extensively 
In  these  references,  we  vill  he  coaeemed  primarily  with  a  brief 
etumnaty  Of  the  assun^ptlons  end  general  methods  of  solution  Involved  In 
propeller  theoty  as  It  is  applied  et  the  present  time. 

Tn  lifting  line  theory,  the  propeller  blades  are  replaced  by 
straight  radial  vortex  lines*  A  free  vortex  sheet  extends  downstream 
from  each  of  the  lifting  lines  forming  an  approximately  helical  surfacft* 
She  propeller  Is  Assumed  to  he  rotating  with  eonstaht  Angular  ▼^eolty 
in  COK  axially  directed  stream  \dioBe  Velocity  may  be  a  function  rtP 
radius  only.  The  flow  will  then  be  steady  relative  to  a  aoordinato 
syetem  rotating  with  the  propeller.  Vie  Hov  la  tiie  nsi^hochood  of  ^ 
the  propeller  is  assumed  to  bo  imaffeeted  by  the  free  enrfhoa  or  by 
extraaecus  solid  boundaries* 

Even  -this  idesaised  model  cannot  be  eolvod  exactly  elneo  tho 
VeXoolty  ttLduced  by  the  vortex  sheets  end  the  posltioa  of  the  sheets 
ere  autually  dependent.  It  Is  therefore  assumed  tliat  the  ibduecd 
Velocities  are  small  compared  with  the  resultant  relative  velocittee 


5 

f 

i 

a-fc  the  liftifig  lines*  Ihe  elements  ©f  the  free  vort^  sheet*  eatt  fheix 
be  assumed  to  t)e  helical  lines  Gf  ConstaSlt  yadltls  an3  Jjltshj  vhere  th© 

Is  determined  by  the  angle  of  the  J'csuitant  flow  at  tha  liftlnft 
line  ix^cladlng  induced  ■\feiocitiea.  This  letter  refinement  complicates 
matters  scmewhat  sicca  the  pitch  ef  the  free  vortex  lines  and  the 
Velocities  induced  at  the  lifting  line  are  still  interdej)enaeat»  hwever, 

•S. 

•  fiolntiOR  may  readily  be  obtalxied  hy  Iteration  • 

She  Justification  for  neglecting  the  axial  deformatlOU  ct  the 
vortex  sheet  is  that  the  velocity  induced  at  the  lifting  line  by  aa 
element  ef  the  sheet  decreases  rapidlly  vith  distance  so  that  aa  errer 
in  the  ftssumed  position  of  the  sheet  becomes  less  criticsl  as  the 
distance  downstreeu  Increases - 

The  relationship  between  the  hound  vortex  strength  end  the 
induced  velocities  at  the  lifting  line  may  he  determined  by  the  Xerbs 
laduetlon  factor  methed^^^*  la  the  epeclal  ease  when  the  Inflow 

Velocity  le  constant  and  the  fitch  of  the  free  vortex  sheet  Is  lar 
deQpesdent  of  radius*  the  circulation  distribution  nay  also  be  detemoined 
by  means  of  the  Goldstein  factore^^^ .  These  methods  will  be  discussed 
ftarther  In  Chapter  4. 

Due  to  the  lov  aapeet  ratio  of  aost  marine  propeller  blades, 

the  nae  of  lifting  line  theory  results  in  unacceptably  large  errors 

unless  supplemented  by  a  lifting  surface  correction  of  some  kind.  Some 

early  attempts  to  e:i{fiiMLn  this  discrepancy  were  based  on  the  application 

(7^ 

Of  tvo-dlaenslonal  cascade  theoxy*  however,  as  pointed  out  by  Lerbs^*', 
this  applicatloa  was  not  Justified.  The  lifting  surface  Correction 
which  is  presently  used  was  first  developed  by  Lbdwi^  ajid  dlazel  la 

and  later  refined  by  Ginzel^^^^ 


.  Their  approach  was  to  find  the 


induced  flow  at  the  saiS-eliord  and  1»  t5b«  tMs  tO  €etermln® 

the  camber  o£  the  blade  sections*  The  plteh  vas  still  to  be  detexmljteS 
from  lif-y.ng  liaa  theety  by  the  requirement  that  the  seetlons  be  at 
zero  angle  of  attacls  relative  to  the  Induce!  fiow. 

Their  theory  is  linearised  to  the  esctexit  that  the  blade 
eorfaoe  Is  assume!  to  lie  in  the  neighborhood  of  a  tzne  helical  surface, 
the  Vortex  system  laid  the  point  where  the  indaeed  velocity  is  to  he 
determined  is  on  the  helical  surface  rather  than  on  the  blade  itself* 

Ole  curvature  of  the  flow  is  related  to  the  derivative  of  the  normal 
component  of  induced  velocity  in  the  chordwise  direction*,  or*  more 
briefly*  the  “dowawash  dejrlvatlve" .  They  assume  a  constant  circulation 
distribution  over  the  chord,  and  with  this  aicipllficatioa  it  Is  easy 
to  show  that  the  downwash  derivative  Is  equal  to  the  dowawash  produced 
by  a  "remainder**  vortex  system  consisting  of  a  line  vortex  representing 
the  blade  outline  and  a  set  of  chordwise  vortices  connecting  the  leading 

ir 

and  trailing  edge. 

Their  results  can  be  e^^ressed  In  terms  of  a  oaiaber  correction 
factor  k  which  is  defined  as  the  ratio  of  the  camber  required  In  three- 
dimensional  flow  to  the  camber  in  two-dimensional  propeller  flow  for 
the  seme  lift  coefficient  •  Vhile  the  theory  can  take  Into  account  the 
contribution  of  the  other  blades  to  the  dovniwaeh  derivative,  this 
affect  waa  neglected  to  simplify  the  eomputationa*  Their  results  shov 
that  the  camber  correction  factor  depends  principally  on  blade  sSTca 
(aspect  ratio)  and  on  the  radial  circulation  distribution* 

la  order  to  apply  their  results  to  propslher  seotlons  which 
do  not  have  a  constant  chordwise  circulation  distribution*  the  chord 
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lengtlvs  e.s'e  «iodtfied  tn  etJCTi  a,  that  the  actual  section  and  the 
Consteint-load  sectioa  would  have  the  same  total  lift  and  downwash 
derivative  In  two-dimensional  llow4» 

After  the  pitch  has  been'  determined  from  lifting  line  theoty 
«nd  the  ceunber  of  the  sections  from  the  ludweig  and  Glnzel  theory,  -the 
design  la  completed  by  auperimposlng  the  velocities  induced  hy  a 
ayimnietrloal  thickness  form  to  those  due  to  the  cambered  ciean  line* 
iv3  in  lljaearlaed  thla  alri'oil  theors',  the  velocities  due  to  the 
thlehness  form  contribute  to  the  local  pressure,  but  not  to  the  lift. 
Finally,  en  ellow3.nae  Is  jaiade  for  viscous  effects  by  adding  a  jproflle 
drag  force  and  by  adding  a  small  angle  of  attack  or  caiaber  Increment 
(or  both)  to  allow  for  the  loss  of  lii't  attributed  to  the  presence  of 
the  boundary  layer.  Both  these  corrections  and  the  velocity  Incrementa 
dae  to  thickness  are  determined  by  a  two-dltoenslonal  atrip  theory  based 
on  the  resultant  Inflow  velocity  from  lifting  line  theory. 

It  has  been  observed  that  propellers  designed  In  this  way  do 

fTT) 

sot  have  the  correct  pitch  In  oany  cases.  To  ejiplaln  this  herbs' *' 

considered  the  possibility  that  the  induced  curvature  siay  not  be 

constant  over  the  chord  and  that  e  pitch  eoi-reetion  might  be  neeessaiy 

(el) 

to  take  this  into  account,  'fo  do  this  the  tfclesinger  '  lifting 
surface  theory  was  applied  approxinately  at  one  point  on  the  blade,  la 
this  theory  the  bound  circulation  is  conceRtrated  at  tbs  l/h  chord 
line  and  the  doi^arv/ash  Is  detemtined  at  the  chord  line,  flba  pitch 
l0  then  adjusted  so  that  the  benodary  oonditlon  at  the  3/^  chord  line 
Is  eatiefted. 
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This  cdrrectlcn  Jts  used  in  the  design  mettiod  described 

Eckhardt  and  MOl’gan^^'^ .  HoVever^  ^aa  Manen  and  Crowley found  tlia^ 

this  Qcrraction  did  cot  seem  to  help  fn  bringing  their  theoretical  and 

experimental  results  into  agreement .  The  author  is  also  of  the  opinion 

that  the  approximations  Involved  In  applying  tlxl%  correction  aro  such 

that  It  is  queationahle  whether  It  can  seiVe  to  Improve  the  accttVac/ 

of  the  iidvrelg  -  Ginzel  theofjr.  tihis  was  illustrated  la  th«  paretent 

fl2^ 

author'©  discussion  to  a  paper  given  by  Morgan  in  1959  • 

Another  form  of  correction  which  has  been  used  prlnelpally 
at  the  JJetherlands  Ship  Model  3aain  la  an  empirical  modification  in 
the  ideal  efficiency  of  the  propeller,  which  results  In  a  change  in 
pitch.  This  Is  applied  principally  to  wake-ad&pted  propellers  and 
includes  the  effecte  of.un3tea(^  fiov^'^'^.  It  is  not  possible  to  say 
how  much  of  this  correction  is  due  to  errors  In  steady-state  propeller 
theory .  * 

Current  Research  in  Steady-State  Propeller  fhecry 

The  fact  that  current  design  methods  are  not  entirely 
reliable  has  resulted  in  e  recent  interest  la  propeller  lifting  surface 
theory.  There  are  many  possible  approaches^  some  of  which  will  be 
discussed  briefly  In  this  section. 

labile  the  Lud■^^leg-5^ixxael  theory  has  a  ta^er  Of  inherent 
siiigjlifying  assvunptlons,  it  is  still  by  no  means  being  applied  to  Its 
Itll  advantage  at  the  present  time.  For  example*  thejr  results  show  a 
very  strong  dependence  of  the  caidber  correction  on  the  radial  load 
distribution,  yet  this  fact  is  ignored  In  cuvrent  design  methods.  It 
appears  that  the  design  curves  given  by  Van  Manen 


are  for  an  optimum 


radial  load  distributioa^  Vhlle  those  appearing  in  Jtokhardt  and  Mbrgan^^  ^ 

ara  foj?  ^  reduced  circulation  at  the  CdteV  5)arl  of  tlie  "blade .  However, 

the  latter  Is  applied  to  propellers  with  both  optimum  ahd  fcon-optimum 

circulation  distribution.  iUtthermore*  the  modification  In  effective 

chord  length  due  to  changes  in  the  chord-load  distribution  is  not  taken 

into  eiccctxnt.  Finally,  the  effect  of  the  other  blades  which  was 

prigitmTTy  neglected  to  save  .numerical  work  ecold  easily  be  taken  into 

account  now  due  to  the  availability  of  high-speed  digital  computers. 

A  reeuialysis  of  the  Luiwelg  and  Glnzel  theory  has  Just  been  completed 

by  Cox^^^^,  and  it  is  possible  that  these  new  numerical  results  will 

result  In  better  agreement  between  theory  and  experiment. 

(l4l 

Following  another  approach,  Alef  '  '  has  been  working  on 

the  exact  application  of  the  Welssinger  theory  to  propellers,  eO-though 
to  the  author *8  knowledge,  no  numerical  results  are  available  as  yet. 

While  this  should  be  a  distinct  Improv’ement  over  the  approximate 
application  of  the  Welsslnger  theoxy^  it  is  still  subject  to  question 
whether  or  not  this  will  offer  emy  improvement  over  the  ludwelg  and  Gihzel 
theory. 

Work  is  also  in  progress  at^fietherlands  Ship  Model  Basin  by 

Sparenberg^^^^  on  a  more  rigorous  Hffclug  surfeuje  theory.  In  that 

* 

reference,  the  baelc  integral  e<;^ntioa  10  derived.  It  Is  understood 
that  work  Is  In  progress  to  solve  the  inte'gral  equation  for  the  special 

case  of  elliptic  blade  outllaas  VitSl  constant  elroulatloc  over  the 

* 


blade  surface. 


GeneraX  Method  of  Approach 

In  the  present  VorJs,  we  ©gnsidel*  the  Solution  o£  the  lifting 
aisjcsface  jjl’Ohlem  for  a  propeller  with  arbitraty  blade  oiStline^  pitch 
distribution  and  clrculatloa  distribution  flpw^atlng  In  au  ajCiaUy  • 
directed  velocity  field.  It  Is  asstimed  that  the  radial  circulation 
distribution  is  given  and  that  the  blade  surfaae  is  to  be  formed  froa 
a  known  mean-line  type  by  detetffilxting  the  eamber  and  pitch  at  each 
radius .  These  two  parameters  are  to  be  determined  by  the  requirement 
that  the  desired  radial  circulation  dletrlbutlon  is  obtained  with  the 
sections  operating  at  their  Ideal,  angle  of  attack.,  She  chordwlse 
circulation  distribution  will  then  be  determined  by  these  two  conditions; 
hy  the  boundary  condition  that  the  flow  be  tangent  to  the  blade  surface^ 
and  by  the  Kutta  condition. 

Ihls  approach  dlffero  from  any  of  the  theories  discussed 
In  the  preceding  eectlons  in  that  no  restrictive  assumptions  need  be 
made  as  to  the  circulation  distribution  or  blade  outline,  and  the 
results  may  be  applied  both  to  open  water  or  to  wake-adapted  propellers . 

The  procedure  is  similar  to  a  method  deirelopeA  by  Ihlkner^^^^^  Cl7)> 
to  determine  the  lift  distribution  of  wing^  Of  arbitraiy  shape,  the  con¬ 
tinuous  distribution  of  radial  and  helical  vortices  is  replaced  by  a 

.  ♦ 

lattice  of  discrete  Vortex  lines.  !Bie  lattice  can  be  considered  as 
foimed  from  a  number  of  "horshos"  vortex  elements  of  constaut  strength 
as  shown  schematically  in  Fig.  The  valocily  induced  at  en  arbitrary 

point  in  space  by  each  lattics  element  eaa  be  determizied  ty  integration 


*The  ideal  angle  of  attack,  or  condition  of  "shock  -free  entrance"  is 
defined  as  the  angle  of  attack  for  which  the  infinite  suction  at  the 
leading  edge  given  by  thin  airfoil  theory  yanishas* 


c 


acfi'Ordtn^  .  By  4etejcminlJi^  ^9 

VfSjo^SJty  at  &  ^rmibei*  of  oontfol  points  on  tKe  tlado  surfaoo  at  tli«^ 
}Qldr*polntsof  the  lattloo  a  sot  of  linear  aqiiatlons  may  ho  formod  arelatln^ 
*tha  strengths  of  th®  lattic®  elements  to  the  shape  of  th®  hlad®  stijface^ 

She  singular  Integral  equation  enooanterel  In  lifting  surfhoe 
thaorf  in  tiierefor®  replaeed  by  a  set  of  tlmultanaoua  linear  eqiiations* 
Slnoe  th®  prooess  is  rexar  largely  aomerloai,  it  ia  not  necassaxy  to  ma]^ 
the  vsual  sln^Xlfylng  assusgptlons  as  to  the  hlads  outline  end  elfcnlatlon 
dlstarfLhutlon. 

The  question  naturally  arises  ae  to  ^diether  the  lattice 
method  vill  converge  to  th®  solution  of  the  Integral  equation  aa  the 
spacing  is  made  sm^ler<  Obviottelyf  If  the  epaeing  Is  made  very  small^ 
the  coefficients  in  the  equations  will  hecoma  large,  dn®  to  the  proximity 
of  the  control  points  to  the  vortex  lines.  Consequently,  from  a 
computational  point  of  view  there  vUl  he  a  point  of  diminishing 
returns  after  'Vdiich  the  set  of  linear  equations  vill  he  too  nearly 
singular  to  he  solved.  The  question  of  vhether  a  etffflolently  accurate 
solution  cein  be  obtained  before  this  taJ»a  place  can  toe  settled  hy 
computing  special  cases  for  \dilcll  the  solntlcn  of  the  integral  equation 
Is  known  and  observing  how  the  error  depends  on  lattice  epaelng* 

This  was  done  by  Ealkaer^^^^  la  the  case  of  wings  of  various 
shapes  and  it  was  observed  that  errors  of  less  than  One  percent  could 
be  achieved  with  lattice®  of  reasonable  sice*  Slnoo  the  convorgonco 


*The  finest  speicing  used  twenty  vortices  over  the  semi-span  eind  eight 
over  the  chord. 


Jropejisies  of  the  lattiee  should  not  ha  altejiad  iajastdoally  by  going 
from  a  plane  to  &  helical  surface,  the  methol  should  he  eS^eoted  to 
Tioaeh  In  the  case  of  a  propelled. 

It  should  be  mentioned  that  this  approach  has  been  studied 
‘to  some  extent  by  Gullloton^^^^  and  Strscheietzl!y^^^\  However,  since 
their  vorh  vas  done  in  the  pre-digltal  computer  era,  it  Is  somewhat 
questionable  whether  a  numerical  solution  on  a  small  enough  scale  to 
he  done  by  hand  would  offer  any  advantage  in  accuracy  over  existing 
results.  Ihis  conclusion  is  based  on  the  results  of  the  present 
vorlc  In  vhlch  It  was  found  that  the  necessary  computations  were  far 
from  triVieJ.  even  for  a  large-scale  digital  computer  and  definitely 
heyond  the  capaelty  of  small  machines,  not  to  mention  humans. 

Basic  Assumptions 

The  assiamptlons  vlU  he  similar  in  part  to  those  made  in 
lifting  li:5e  theory  fis  described  In  the  beginning  of  this  chapter. 

Ihu  fluid  Is  assumed  to  ba  frictionless  and  incompressible  and  the 
flovr  in  the  neighborhood  of  the  propeller  is  assumed  to  be  unaffected 
hy  a  free  surface,  extraneous  solid  boundaries,  or  cavitation.  The 
Inflow  Velocity,  as  la  lifting  line  theory,  is  assumed  to  he  axial  and 
a  function  of  radius  only. 

The  free  vortex  sj'stem  is  assumed  to  lie  on  a  helical  surface 
vhose  pitch  is  determined  from  lifting  line  theory  with  the  same  radial 
load  distribution.  The  pitch  of  this  helical  referenee  surface  Bay  be 

o  « 

a  funetion  of  radius .  Ihe  blade  surface'  Is  assumed  to  be  approScimatsS^r 
on  the  helical  reference  sarfeec.  ^e  proV.^  Is  linearized  to  the 

extent  that  the  boundary  condition  is  applied  on  the  helical 
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guj;fac«  rather  them  on  the  blade  Iteelf  emd  the  Induced  velocities 

© 

ajff  assiuned  to  be  small  l^latlv#  to  the  resultant  inflow*  As  In  Ziftlnff 
line  theory,  the  flow  is  assiuned  to  lie  on  cylindrical  surfaces  concentiic 
Vith  the  propelled  axis  of  rotation.  This  assuiiQ)tlon  is  obviously  not 
Very  realistic  near  the  tip  of  the  blades,  but  should  be  reasonable 
elsewhere  for  moderate  propeller  loadings. 

It  Is  assumed  that  the  Kutta-condition  holds,  l.e.,  that  the 
boond  circulation  is  zero  at  the  trailing  edge.  It  is  also  assumed 

that  the  bocnd  circulation  is  *ero  at  the  blade  tip  and  at  the  hub 

» 

radius,  and  that  the  boundary  condition  of  zero  radial  velocity  at 

the  hab  eyllnder  can  be  disregarded.  These  last  two  assumptions 

conceralns  the  hub  are  by  no  neans  essential  to  the  vortex  lattice 
% 

method,  and  It  Is  beHeved  that  a  more  accurate  representation  of  the 
hub  effect  fan  be  added  at  a  later  time. 

Outline  of  Results 

In  order  to  apply  the  vortex  lattice  method,  the  velocity 
.Induced  at  an  arbitrary  point  In  space  by  a  set  of  helical  or  radied 
vortices  is  seeded.  Expressions  for  these  are  derived  in  Chapters  2 
and  3  respectively,  and  methods  of  con^iutation  emd  error  estimates  are 
discussed.  In  Chapter  h  vortex  lattice  methods  are  applied  to  solve 
the  lifting  line  problem,  both  lor  optimum  propellers  in  homogeneous 
How,  emd  lor  non-optimum  ©r  vahe-adapted  propellers.  This  is  inaluded 
to  tndteeto  to  some  extent  the  convergence  properties  of  the  lattice 
method  by  Comparison  with  known  results .  These  results  are  also  needed 
in  the  solution  of  the  lifting  surface  problem  for  symmetrical  blades. 


In  Chapter  5  a  iattlca  solution  ^s  developed  fotf  t)rqpellers 
ot  generally  arbf traiy  “blade  outline,  aeotioa  “type^  andi  l-adleuL  clscdlatlon 
distribution,  and  in  Chapter  6  these  results  are  specialized  In  the 

« 

case  Of  propellers  vlth  symmetrical  blades*  In  the  latter  case,  the 
jcesultlng  symmetry  greatly  simplifies  the  computations* 

finally,  in  Chapter  7  numerical  Results  for  camber  and  pitch 

Corrections  are  presented  and  compared  with  results  according  to  the 

(»«• 

ludicg  and  Glnzel  theory. 
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*  CHAPTBR  2 

THE  VBLOCITY  UHJCgD  BY  HKJCAL  VORTK  LUIS 

Introduction 

In  this  chapter  the  p:r6bleD.  of  determining  the  veloalt^ 

Induced  at  an  arbitrary  point  in  space  t>y  a  set  of  helical  vortices 
will  be  considered.  It  vUl  be  assumed  that  the  vortices  are  of  true 
helical  shape,  l.e.,  that  their  radius  emd  pitch  remains  constant^ 
and  that  there  will  be  g  vortices  of  equal  strength  synmetrically 
located  around  the  circumference.  The  axial  extent  of  the  set  of 
vortices  may  either  be  finite,  as  in  the  case  of  a  vortex  segment 
lying  on  the  blade  su.rfeu:e^  or  semi-infinite  as  in  the  case  of  the 
free  vortex  system  extending  downstream  from  the  trailing  edge  of  the 
blade « 

The  velocity  induced  by  a  vortex  line  of  eurbitrary  shape  may 

be  expressed  in  terms  of  an  integral  taken  ed.ong  the  vortex  line  by 

(19) 

means  of  Blot-Savart '  s  Law^  .  Expressions  for  these  integrals  in 

(2^)  f  ^  (22^ 

the  case  of  helical  vortices  have  been  derived  by  Bets'  St^heletsky'  ' 
and  others.  However,  since  the  derivation  is  very  short,  it  will  be 
■  included  here  for  convenience  since  these  references  are  not  widely 
available.  This  will  also  serve  to  establish  the  notation,  ^ich  is 
by  no  means  universal. 

Since  these  integrals  cemnot  be  solved  ^e;q>llcltly^  other  methods 
have  generally  been  used  in  the  past  to  obtain  the  induced  velocity 
CGo^nehts.  In  lifting  line  theory,  for  exanple,  the  velocity  Induced 
on  the  lifting  line  by  a  set  of  semi -infinite  helical  vortices  can  be 
reduced  to  the  two-dimensional  problem  of  finding  the  Velocity  Induced 
by  a  helical  vortex  of  infinite  axial  extent,  as  was  first  shown  by  Betz'  , 
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HU.!  can  be  treated  as  a  tvo-dlaenalooal  potential  prdblea  and  aolutlons 
for  the  case  of  a  set  of  helical  Tortez  lines  bare  been  obtained  by 
Leibs^^^  and  in  the  case  of  a  set  of  true  helical  surfaces  by  Gtoldstein^^^. 

Howerer,  in  a  vortex  lattice  approxlmtion  to  the  lifting-surface 
_pr6bleB,  the  velocity  induced  at  an  arbitrary  point  in  space  by  a  segnent 
of  a  helical  vortex  line  Bust  be  determined.  Since  this  is  nov  a  three- 
dlBKnsional  problen,  the  Biot-Savart  integrals  vould  appear  to  provide 
the  best  way  of  obtaining  the  induced  velocities# 

In  the  case  of  a  finite  interval,  the  integration  auy  be 
perforsed  by  nusKrical  methods  as  viH  be  discussed  later.  In  the 
seed.- infinite  case,  nuaierlcal  integration  aay  be  used  up  to  a  sufficiehtly 
large  distance  downstrean  at  vhlch  point  the  reaminlng  value  of  the 
integral  to  infinity  can  be  estlaiated*  Both  of  these  steps  Introduce 
errors  nonaaUy  defined  in  muaerical  analysis  as  “truncation  errors". 
However,  in  this  application  the  tern  "integration  error"  will  mean  the 
error  introduced  by  the  nuaurical  integration  fonaila,  while  "truncation 
error"  will  i^fer  to  the  estiasite  of  the  integral  to  infinity#  Both  of 
these  errors  will  be  considered  in  detail  later  in  the  chapter. 

Ihe  Induced  Velocity  Cosponents  Detemlned  by  Biot-Savart  *s  Law 

As  shown  In  Pig.  (2.l)>  a  rlj^t-handed  cartesian  coordinate 
systeai  is  located  with  the  x  axis  along  the  propeller  axis  of  rotation 
with  positive  direction  downstrean.  The  y  axis  pMses  through  the 
control  point,  i.e.,  the  point  in  space  where  the  velocity  is  to  be 
detemined.  A  cylindrical  systen  (x,  r,  9)  is  oriented  So  that  the 
line  X  s  0,  9  B  0  in  the  cylindrical  system  corresponds  to  the  y  axis 
in  the  caurteslan  system# 


FIG.  2.1  COORDINATE  SYSTEM  AND  NOTATION 
FOR  .  HELICAL  VORTICES 
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Tbere  will  be  g  bellcal  Yortlces  (one  firca  each  blade)  lAlcb 
baTe  the  foUovlng  properties: 

a)  !l3ie  vortices  all  start  vlth  the  same  axlsLl  coordinate 

X  .  radial  coordinate  r  ,  but  with  different  angnlar 
o  o 

coordlmites  9  »  9^,  P  •  1^  2,  ...  g. 

b)  The  vortices  are  of  constant  radius  r^,  and  constant 
pitch  angle  3^^. 

Blot-Savart's  Lav  aaiy  be  written 


r  r  dJ  x"? 

kJ 


(2.1) 


Q 

where  r  *=  vortex  strength  (ft  /sec) 

S  •  vector  distance  froK  vortex  element  to  control  point  (ft) 
dl  -  vector  element  of  distance  along  the  vortex  (ft) 
u  -I  vector  Induced  velocity,  (ft/sec) 

The  distance  3  has  the  following  x,  y,  and  x  ccajponents: 

3  »  [-x^  -  r^qp  tan  3^^^,  r  -  r^  cos  (9  +  <Pp)»  -  sin 

(9  +  9p)]  (2-2) 


where  9  Is  the  angular  coordinate  measured  from  9p  as  shown  In  Fig.  (2.1). 


The  vortex  element  di  Is 

di  *  l^tan  "lo'  -  sin  (9  +  9^),  cos  (9  +  9p)]  r^  d9 
The  cross-product  di  X  3  Is  as  follows 


(2.3) 


di  X  3  -  r^  d9 


1  J  X 

tan  3  -sin  (9^9  )  cos  (9*9p} 

-x^-r^9  tan  3  r-r^  cos  (9^p)  -r^  sin  (q^Pp) 
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r  -  r  cos  (9  9  ), 

tan  [slii  (9  t  9p)  V  cos  (9  +  9p)] 

-  cos  (9  i  9p), 

tan  3^0  [  r  -  cos  (9^-h  9^)  -  9  sin  (9  +  9p}J 

-  Bin  (9  +  9p)  (SX) 

axid  the  scalar  quantity  3“^  is 

5^  ■  j^(x^  ^io^^  ^  ^o  ^ 

(2.5) 

Substituting  (2.2)  tlrcxTugh  (2.5)  in  (2.1)  and  sumning  over  the 
g  blades  gives  the  foUoving  expressions  for  the  cucial,  tangential,  and 
radial  velocity  coiponents 

r  X*  ^ 

""a  =  IJr  ?  t^'o  "  ^  ^ 

rr 

'  "t  “  ^  4^  ?  ^ 

-  sin  (9  +  9^)  (x^  +  9  tan  3^,^)]  <ip  (2.?) 

Fr  ® 

-  Sr  1^4'^  [-  f'’o  ’P  tan  .  x^)  cos  (»  .  9^)  ■ 

♦  tan  r^  sin  (9  4  9p)J  ckp  (2.8) 

The  above  equations,  after  due  changes  in  nosienclature,  are 

S  "z.  (22) 

in  agreesient  with  Sti^heletalQr's'  *  formula  3^.  Furthexmore ,  in  the 
special  case  vhen  one  of  the  helix  starting  aixgles,  9^,  as  well  as  the 
axial  starting  points,  x^,  are  zero,  these  expressions  agree  with  those 

•  (yiS  (■\\ 

given  by  Betz'  and  Lerbs'  '.  ISiis  latter  case  corresponds  to  the 


the  velocity  components  at  a  blade  In  propeller  lifting  line  theory. 

Equations  -  ^2»&)  can1)e  made  non-dimensional  im  terns 

of  the  following  variables 


n  =  rjT 

t  -  x^r 

-  imru, 
u 


© 


(2.9) 


The  non-dimensional  induced 


velocity  components  u  can  then  be  written 
[t1  -  cos  (<p  +  9p)]  dtp  (2.10) 

[ten  1  -  “H  cos  (9  +  tp^)} 


-  sin  (9  +  9^)  {5  +  T1  tan  3^^  9}]  d9  (2.11) 

I ,  ^  [*“■  ho  *  ’’p)  - 

p=l 

-  {9  tan  3j^q  +  ?/T1)co8  (9  +  9^)]  d9  (2.12) 

where  the  denominator  in  each  of  the  Integrals  above  Is 

.  [(5  *  n  9  tan  +  1  +  Tjp  -  2  T)  cos  (9  +  SPp)J^^^ 


(2.13) 

The  non-dimensional  velocity  u  Is  related  to  the  Lefbs^^^ 
Induction  factors  1  by  the  relation 

u  -  5^  ^  (2.14) 

nie  reaison  for  selecting  a  different  non-dimensional  fozs  is 

s 

based  on  a  consideration  of  numerical  accuracy.  The  total  velocity  at 
a  control  point  is  to  be  obtained  by  sunning  the  velocities  induced 
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by  the  elements  of  a  lattice  system.  TSxe  velocity  induced  by  the  nearby 
elements  will  become  veiy  lajcge  as  the  lattice  spaK^lcg  becomes  small, 
so  that  these  must  be  caiQiuted  to  an  increasingly  large  number  of 
significant  figures  for  a  prescribed  accuracy  in  the  resultant  velocity. 

The  quantity  u  will  tend  to  infinity  as  (l  -  11)”^  as  I)  -*  1,  hence  , 

requiring  a  fixed  accur8K:y  in  u,  (say  three  decimal  places  correct) 
is  equivalent  to  requiring  a  higher  percentage  accuracy  as  the  imgnltude 
of  u  Increases* 

On  the  other  hand,  the  Induction  factors  remain  finite  due 
to  the  factor  (1  -  T|),  so  that  if  the  numher  of  decimal  pieces  in  the 
canputation  of  the  Induction  factors  is  sufficient  for  the  nearby 
elements  of  the  lattice,  the  induction  factors  for  the  distant  elements 
will  be  Tumecessarily  accurate. 

In  general,  the  velocity  coaQ>onent  noxmal  to  a  particular  boundary 
is  to  be  determined.  Let  (/,  m,  n)  be  the  (x,  y,  z)  coiiQ>onents  of  a 
vinlt  vector  nozmal  to  the  surfMe.  The  non-dimensional  normal  velocity 
is  then  given  by 

u=iu_+mu+nu.  (2.1^) 

n  a  ’  r  w 

Pbr  purposes  of  computation,  it  is  convenient  to  expreea  the  integral 
in  the  following  form 

(2.16) 

where  the  c's  and  d's  are  constants  in  the  integration,  but  depend  on 
the  blade  index  p.  From  (2.5)  these  constants  can  be  written  as 


(Cj  +  cos  y  +  8^  y  c^  (p  cos  ip  +  Cj  (p  sin  y)  dtp 

(d^  9^  +  dg  <p  +  d  +  d|^  cos  9  +  d_  sin  9)^'^^ 


ft 
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°3  “  “  ®3t 

.  ic^  +  «  +  n 


By  e^qwmding  sin  (9  -•■  9^)  and  cos  (9  ■*■  9^)  In  (2.10) 
coefficients,  the  foUovlng  e3q>res8lans  ere  obtained 

1 


(2.17) 

-  (2.12)  and  collecting 


°1mi  “ 

V  -  ’I  »p 


'6.  *  '7.  ■  ° 


axial  coB^ponent 


C|^^  «  Tf  tan  sin  9p  -  ^1  ?  cos  9^ 

c^y  «  tan  COS  9^  +  11  5  Bin  9p 

cgj.  =  71^  tan  cos  9^ 

Cyr  -  f  tan  P^^  sin  9p 


radial  coiQ>onent 


.  1)  tan  P^o 

-  -11^  tan  Pj^Q  cos  9p  -  1)  5  sin  9j 

-  11^  tan  Pj^Q  sin  9p  -  11  §  cos  9^ 

-  11^  tan  P^^  sin  9^ 

=  -11^  tan  P_  cos  9^ 


"tangential  ccagponent 


Ihe  coefficients  of  the  denoalnator,'  ^dilch  are  the  saae 


all  three  coaponests,  are 


d^  -  if  tan  2  Pj^ 
dg  -  2  9^  if  tan  ^ 
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S  “  TI^  +  1  +  Tip 

dj^  -  -2  T)  cos  9^ 
d^  «  2  T)  sin  9p 


(2.18) 


By  considering  the  non-existent  constants  c^,  c^,  dg,  and 


to  be  zero,  and  by  defining  a  function  (9)  as  foUovs 


F^  =  9 


^2  =  9 


^3-1 


F^  =  sin  9  m  cos  9 


cos  9 
9  sin  9 


(2.19) 


a  more  coiig)act  ejqpreasion  for  is  obtained 

[j 

n  p*»l  7 


(2.20) 


If  the  Integral  is  to  be  evaluated  by  an  I  point  integration  formula  vlth 


weights  V^,  (2.20)  may  be  written 


7 


u 


s4  ^ 


(2.21) 


where  means  P^  (9^).  ®iis  is  a  convenient  form  for  use  with  a  digital 

con^uter.  As  is  described  in  Appendix  (  A  ),  values  of  F^^  may  be  coaq^uted 

and  stored  in’  a  table  so  that  only  the  constants  c  cmd  d  need  be  computed 

n  n 

for  each  Integration.  Ihls  results  in  a  large  saving  in  conqputaticn  time, 
which  is  laq^ortant  since  the  evaluation  of  these  Integrals  r^resents 
the  major  part  of  the  numerical  work  in  obtaining  lifting  surface  solutions 
by  a  lattice  method. 

Ihe  velocity  cosgKment  nozmal  to  a  true  helical  surface  can 
be  detezvlned  by  substituting  the  components  of  the  unit  nozmal  in  (2.1$). 


Choosing  the  positive  direction  for  the  noroiBLi  to  he  directed  upstream, 
l.e..  In  the  direction  In  which  a  propeller  would  nonasUy  he  developing 
thrust,  there  follows 

i  «  -cos  m  •  0  n  =  +sln 

u  =  -5  cos  P.  +  u.  sin  P.  (2.23) 

Z1  A  X  V  ,  X 

where  P^  Is  the  pitch  angle  of  the  helix  at  the  control  point  radius  r. 
Integration  Error 

In  the  case  of  a  semi-infinite  vortex,  equations  (2.10)  -  (2.12) 

or  (2.22)  nay  he  solved  hy  numerical  Integration  up  to  some  angle  9^, 

and  the  remaining  contribution  from  9^  to  (O  estimated.  In  this  section 

the  error  introduced  In  the  numerical.  Integration  from  0  to  9^  will  he 

considered.  HheBe  results  nay  he  applied  equally  well  to  the  integration 

of  vortex  segments  of  finite  length  on  ^e  blades. 

To  get  some  idea  of  the  spacing  required,  the  error  in  the 

axial  coiiq>onent  will  he  derived  In  the  case  of  nunerlcal  integration  hy 

•  (2k) 

Slapson's  Buie.  She  e3g;>resslon  for  Slapson's  Buie'  ,  Including  the 
error  term,  Is 

J  f  (x)  dx  »  I  (fj,  +  ^^*^3) 

*b 

where  the  total  length  of  -the  Interval  Xg  -  x^  =  2h,  and  <  5  5  *2* 

Note  that  x  and  §  refer  to  the  variable  of  Integration  In  general,  not  to 
the  eoordiiiatfs  defined  Ih  Fig.  (2.l)* 

If  the  magnitude  of  t&e  maximum  allowable  error  In  one  revolution 
of  the  Integration  Is  c,  the  i^miber  of  Sinqpson's  Buie  elements  per 
revolution  Is 


-22 


and  the  maxlnim  error  per  element  Is 


sh/tT 


oIV 


(?) 


so  that  the  maximum  Integration  spacing  Is: 
90  € _ 


h  = 


n  (5) 


(2.25) 


(2.26) 


If  (;}  Is  Interpreted  as  the  maximum  value  In  the  Interval^  c  will 
be  an  upper  bound  on  the  error  for  a  spacing  h. 

The  fourth  derivative  of  the  Integrand  of  (2»10}  after  an 
elementary,  but  lengthy  calculation,  may  be  ergpraSsed  as  foUovs  In 
terms  of  the  notation  of  Fig.  (2.1).  > 


g 

(«P)  =  ^  ^  "  r  cos  Dg  +  r  sin  cp^  dJ  (2.2?) 

where: 

®2  “  ^^1  ^2^  *** 

^^1  ^2^  9  -  0^  cos  9 


=  59.0625  s*^  -  78.75  s’^  s^'  +  11.25  S"^ 

+  15.0  S'  S'"  -  1.5  s’^/^  S''' 

Cg  =  -22.5  3"^/^  s'^  +  9.0  3"^^^  S"  +  s"^/^ 

Cj  =  52.5  3"^/^  s'^  -  45.0  S'  s"  +  65"^/^  s’"  -  6  3"^/^  S 

# 

2 

S  =  d+  e9f9  +  g  cos  9  +  h  sin  9 
S'  a  e  -f-  2f9  -  g  sin  9  h  cos  9 
s"  s  2f  -  g  cos  9  a  h  sin  9 

♦  f 

S  a  g  sin  9  -  h  cos  9 

iv  .  ^ 

3  a  g  cos  9  h  sin  9 
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2  2  2 

d-x  +r‘^r‘^ 
o  o 

e  >  2  X  r  t«n  3. 

o  o  lo 

2  2 

f  -  r  tan  3. 
o  ,  lo 

g  «  -2r  r  608  CP 
o 

h  s  2r  r  sin  (p 
o 

Uxxfortunately,  many  of  the  tezma  in  the  above  ej^resalon  are 
of  the  same  magnitude,  ao  that  It  doea  not  seem  foaaible  to  obtain  a 
simple  upper  bound  for  f  (cp)  without  being  unreasonably  coSiservatlve . 

The  above  equations  were  therefore  prograiiiie4  for  an  IBM  6$0 
and  a  few  aaiqple  cuzvea  of  f (qp)  were  coiqputed. 

Ilg.  (2.2)  shows  a  8aiq>le  plot  of  for  a  three 


and  flve-bladed  propeller  with  T)  ■  2  and  3^^^  ■  20*^.  Frcm  (2.26)  this 
la  seen  to  be  Inversely  proportloxial  to  the  spacing  required.  This 
indicates  that  the  spacing  after  one  revolution  can  be  about  ten  tines 
the  Initial  spacing  for  constant  error. 

When  T]  Is  close  to  one,  the  fourth  derivative  Is  Initially 
veiy  large.  The  following  values  are  for  1)  «  .95#  “  20°,  -and  g  »  3 

f'"*  (t)  [f ''' 


0 

3.31  X  10^ 

2ho 

3 

8.99  X  lo’^ 

91 

6 

2.74  X  10^ 

72 

In  order  to  guarantee  an  error  of  less  than  .0001  per  revolution 
In  this  case,  an  initial  spacing  of  about  *0^  degrees  would  be  required, 

a 

while  for  T|  »  2  the  Initial  spacing  could  be  2.6  degrees.  After  one 
revolution,  a  spacing  of  around  30  degrees  would  be  sufficient,  regardless 
of  the  value  of  1). 
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1  ^ 

J  f  ax  -  y  Wv  f  (a^)  +  fconat)  f*  (.5)  J2.26) 

where  the  weights  and  ordinates  are  given  In  Table  A-2  In  the  Appendix. 
While  this  formula  would  be  very  cvunbersome  for  a  hand  calculation  due 
to  the  Irrational  wel^ts  and  unevenly  spaced  ordinates,  on  a  digital 
computer  this  would  teihe  the  same  length  of  time  per  point  as  Simpson's 
Rule  and  would  have  a  much  higher  degree  of  precision. 

As  a  result  of  calculating  a  large  number  of  Induced  velocity 
integrals,  It  was  observed  that  in  all  cases  a  larger  spacing  between 
points  could  be  used  with  the  5  point  Gauss  Rule  than  with  Simpson's 
Rule.  The  eidvantage  was  greatest  for  vailues  of  1)  near  vinlty  where 
the  Gauss  rule  spacing  could  be  five  times  as  large  as  the  Simpson's 
rule  spacing  for  equal  accuracy. 

As  a  result  of  these  saaiple  calculations.  It  was  also  noted 
that  \dien  |l  -  T||  was  It  was  not  necessary  to  decreue  the 

spacing  when  Integrating  the  blades  other  than  the  index  blade.  By 
using  a  wide  spacing  for  the  non- Index  blades,  a  significant  reduction 
In  conputatlcn  time  could  be  achieved,  particularly  for  five  or  slx-bladed 
propellers . 

Although  the  spacing  required  for  a  i>artlcular  aceuxacy  dq^ends 
on  g,  11,  tan  and  x^,  there  Is  very  little  to  be  gained  In  Including 
a  parameter  which  has  a  relatively  small  effect  on  the  required  spacing 
since  the  time  spent  selecting  and  manlphlatlng  blocks  of  stored  tables 
my  affect  any  time  savings  In  the  actual  Integration  process.  It  appears 
as  though  the  critical  parameter  la  ]l  -  T)!  and  that  the  effect  of 
g,  tan  3^^  and  x^  on  the  required  spacing  can  be  Ignored.  It  also 
appears  reasonable  to  divide  |l  -  T)l  Into the  following  three  regions: 
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> 

*02  < 

1  1  -  11 1  5  .10  ^ 

Fine  facing 

.10  < 

|l  -  1!  <  -25 

Hedlua  Spacing 

.25  < 

II-  Hi 

Coarse  Spacing 

(2.29) 

Values  of  ll  - 

1)j  <  .02  were  not  considered,  since  %hls  Is  the 

smallest 

value  tiiilch  would  be  obtained  with  the  vortex  lattice  syatens  anticipated. 
Table  (A>I)  In  the  Appendix  contains  a  list  of  angular  intexvcU.s  which 
\dien  divided  Into  3>polnt  Gauss  ordinates  will  produce  values  of  the 
Integrals  correct  to  3  declaml  places. 

Truncation  Brror 

An  uj^er  bound  on  the  error  Introduced  by  truncating  the 
Integration  at  scsse  angle  9^  can  be  obtained  as  follows: 

Hie  Integral  to  be  estimated  Is: 

•  e 

®  ^a  *  ^  J  ^  ^  ^  (2.30) 

The  denominator  can  be  sli^llfled  as  follows: 

=  [(S  +  71  cp  tan  +  1  +  7)^  -  27)  cos  (<p  + 

>  7)^  cp^  tan^  3^^  ‘  (2.31) 

Subsltutlng  (2.31)  In  (2.3O)  and  replacing  -cos  (9  +  V^by  1, 


s 

71  +  1  V 

g(7l+U 

0  U 

a 

71^  tan^  3^^  p«l  ? 

1 

1 

2ir  tan3  3^^  9^^ 

(2.32) 


Similarly  from  (2.11)  the  tangential  velocity  estimate  is 

g 

I 

r-i  *i>. 


6  ^4. 

< 

,  t 

tan^  P 


r  (i+T|)  tan  +  g+f)  tan  cp 


io 


g 


D 


1  + 


t  cp^ 

(1+Tl)  tan  +  g 


211  tan  p^^  9^" 


f2.33) 


11  tan  p  <p^ 

‘ii 

For  exas^le,  if  11  =  1,  tan  p^^^  *1,  g  ®=  0  and  g  =  3>  "the 
maximum  error  introduced  by  truncating  the  integration  after  a  revolutions 
(cp^  =  2Tm)  is  shown  in  Table  2»1. 


Table  2.1  Truncation  Error  Bound 


No.  of  Revolutions  n 

I«  m*. 

1  A  u^  max. 

1 

.0760 

.5500 

2 

.0190 

.2570 

3 

.0084 

.1650 

k 

.004? 

.1240 

5 

.0030 

.0965  . 

6 

.0021 

.0815 

13 

.0005 

While  this  estimate  is 

very  conservative. 

particularly  in  the 

case  of  the  tangential  velocity. 

it  illustrates  the 

fact  that  after 

2  or  3  revolutions  the  error  decreases  very  slowly. 

On  the  other  hand. 

after  a  few  revolutions,  the  value  of  the  integral  to  infinity  can  be 
accurately  estimated  as  follows? 


For  large  values  of  cp^; 
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tan^  8^ 
gin  ® 


6 

-5. 


cos  9 


cos  q)  d9 

■t  ^ 


(2.3i>) 


The  last  two  integrals  In  (2.34)  ^an  he  reduced  to  the  Sine  Ingetral 
j^Sl  (9)J  and  Cosine  Integral  fci  (9)J  which  are  tabulated  functions. 
However,  If  ♦he  blades  have  equal  angular  spacing,  the  suns  over  cos  cp 
and  sin  cp^  are  zero  so  that  only  the  first  term  remains.  In  this  case 
the  estljoated  value  of  the  Integral  becooMs: 


t  u  fia  ■  —  ■  I - - 

2T1  tan3  8^^ 


(2.35) 


Similarly,  the  approximate  value  of  the  tangential  velocity  is: 


4  u .  — y  ^ 


(2-36) 


An  upper  bound  on  the  error  Introduced  by  using  (2.35)  can 
be  obtained  as  follows: 

Assume  that  the  actual  value  of  and  the  approxlamite 
value  differ  by  the  factor  [l  +  e  (9)],  where  e  «  1.  Then 

T  1  r  1  ^  •  cos  (<p  +  V  )  4p 

IT)  -  cos  (9  +  9  )  d9  7  e  (Tj  -  cos  (9  9  )  d9 

-T - 5 - ? -  -11  - , - 5—-? -  +  •••■ 

*  V  w  P,„  v-’  5>.  T|^  W  Sj  9^ 


=  6  a  +6 

a 


(2.37) 


Where 


^  tan3  8^^ 


«  (T)  -  cos  (9  +  9  )} 


(2.38) 


-30- 


iB  the  error  in  the  approximation  6  u  .  If  e  is  the  maximum  value 

&  max 


of  e  («p)  in  the  interval  <  <p  <  •»  4  can  be  written: 

Is  (T1  +  1) 
j  6  I  <  gmx  ^  ' 

21?  tan^ 

!nie  quantity  can  be  estimated  as  follows: 


(2.39) 


(1  +  e)(‘r)^  tan^  =  [(5  ♦  11  9  tan  +1+11^ 

-271  cos  (9  +  <5Pp)p^^ 


Solving  for  c: 


[(5  +  11  9  tan  +  2  +  11^  -  211  cos  (9  + 

'  “  T?  9'  tan3 

[5^  +  211  5  9  tan  p  +  if  9^  tan^  810+  1  ^  ^  ♦  Sll]^/^ 

1*1  <  - - - - - —  -1 


9  tan'^  P.  +  0 


(2.1^0) 


O 

In  the  case  idien  ^  «  0  and  9  »  1,  the  3/2  power  in  the 
numerator  can  be  eiqpanded  giving  the  approximate  result; 

1*1  5  I  ^2^  ■'^^1  (2.41) 

'  '  2  if  9^  tan^  ^ 

The  maximum  value  of  t  is  when  9  9.  >  Substituting  this  in  (2.39) 

gives  the  result 


5  i 


Cl  .  If  .  2T|)(TI  .  1] 

ii  s  4 

11“  Un^ 


(Z.kS) 


Solving  for  9^ 


3/  (1  Tf  2T1)(11  +  1) 
4  6  11  tan^  p. 


(2.43) 


Taking  the  same  numerical  exa]i^>le  as  before,  if  11  •  1  tan  ^ 

5=0  g  =  3  and  6  ■  .0005,  (2.43)  gives  the  result: 

9^  «  13.8  radians  2  revolutions. 


According  to  Table  (2.1),  it  would  require  13  rerolutlona  to 
obtain  the  same  accuracy  If  the  munerlcal  Integratlona  were  used  entirely. 
Since  Table  (2.1)  represents  a  very  conservative  estimate,  the  actual 
saving  in  using  the  approximate  value  of  the  integral  from  <p^  to  *  Is 
sooeidiat  less. 

Equation  (2*^3)  a  similar  one  for  the  tangential  velccity 
could  be  used  to  determine  cp^.  However,  this  is  eU.so  a  little  conservative, 
so  that  It  is  more  efficient  to  use  a  more  en^lrlcal  way  of  deciding 
when  to  stop  the  numerical  Integration.  OMs  is  done  by  estimating 
the  value  of  the  integrals  to  infinity  from  (2.35)  and  (2.36)  after 
each  revolution  in  the  numerical  integration  has  been  congpleted.  When 
two  successive  estimates  agree  to  the  desired  tolerance,  the  approximation 
of  the  Integral  Is  assumed  to  have  converged. 

lumerlcal  Results 

In  order  to  check  the  preceding  results.  Induced  velocity 
coiig>onents  were  computed  corresponding  to  thtee  numerical  exaaqples  given 
by  Vrench^'^^.  The  velocity  ecsiponents  obtained  by  numerical  Integration 
converted  to  Induction  factors  by  (2.l4)  agreed  to  four  decimal  places 
with  Wrench's  values,  idilch  was  the  total  nuzaber  of  places  given.  Checks 
against  gross  errors  were  made  by  comparing  Induction  factors  over  a 
wider  set  of  parameters  with  the  tables  given  by  Morgan^^^  and  In  all 
cases  the  agreement  was  satisfactory. 

In  addition,  large  nuabers  of  coagutatlons  were  made  to 
determine  the  optlaam  integration  spacing  as  was  discussed  previously, . 

however,  since  these  results  are  of  limited  usefulness  once  the  spacing 

» 

criterion  has  been  established,  this  data  will  not  be  reported* 
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CHAPEIR  3 

THE  ygiOCITr  mXJGlD  BY  RADIAL  VOREK  LIWS 


!Ihe  velocity  induced  by  a  atltii^t  radial  vortex  segment  of 
constant  strength  can  be  obtained  by  integration  using  Blot-3avart's 
Lav.  While  ^he  helical  case  vas  somevhat  complicated  due  to  the  necessity 
of  using  numerical  integration,  the  expressions  obtained  for  the  radial 
case  are  very  3iJig)le  and  may  easily  be  Integrated  e^qplicltly. 

yi 

The  notation  to  be  used  is  shown  in  Fig.  and  is  substantially 

the  same  as  Fig.  2.1.  A  set  of  g  radial  vortex  lines  are  located  at 

angles  and  extend  from  r^  to  r^*  The  remaining  notation  is  the  same 

as  in  the  helical  case,  except  that  the  variable  of  IntiSgratlon  is  now 

r  Instead  of  m. 
o 

The  ccmponents  of  the  vector  element  of  vortex  line  dJ*  are 
d2  =  ^  0,  dr^  cos  <Pp,  dr^  sin  (3*1) 

and  the  distance  from  the  vortex  element  to  the  control  point  is 

3  «  [-x^,  r  -  r^  cos  9^,  -  r^  sin  9J  (3-2) 

Substituting  these  quantities  into  the  esqpresslon  for  Biot-Savsut ' s 
Lav  (2.1),  the  foUovlng  e:q>resslons  for  the  velocity  coiiQ)onents  are 
obtained 


r  ?  V  “‘•'Pp  ■* 


r  ^  -  2r  r  cos  cp  ) 
o  o 


372 


u  =  0 
r 


r  T  V  \  ’’p 


I, 


■,  ^«1  (x  ^  r^  ♦  r  ^  -  2r  r^  cos 
X  o  o  c  p 


(3.3) 


As  in  Chapter  2, these  can  be  e^gpressed  in  terms  of  the  non>dlmenslonsil  quantities 

l»TTr  .  ,  . 

(3.^) 


T)  -  xJt  5  »  xJt 


U  =s  u 


a 


FIG.  3,1  COORDINATE  SYSTEM 
FOR  RADIAL  VORTICES 
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resviltlog  in  the  following  expressions 
g 


sin  q> 


P 


g 


“t  ‘5  I 


COS  9 


p=l 


1  ^ 


L 


p  j 


372 


(3.5) 


where  the  denominator  is 
g3/2  =  ^  ^  ^  ^ 


Equations  (3*5)  can  be  integrated  to  give  the  fcllowiqj 


u 


where 


=  -y 


sin  cp  I 
P  P 


I  = 


T)  -  cos  cp 


^  (5^+oin^ 


\  =  ? 

T1 


G 


cos  cp  I 
p  p 


(3.6) 


f  +  sin^  cpp  0 


(3.7) 


I  - — - ? 

P  2  (T)  +  cos  cp^r 


2  2 
5*^  +  Sin  cp  -  0 

^  '♦'p 


(j.?J 


The  latter  form  corresponds  to  the  case  \dien  the  vortex  segment  coincides 
with  the  y  axia,  at  which  point  the  velocity  is  zero  as  can  be  seen 
from  (3.5) 


pitch 


>dilch 


As  in  Chapter  2,  the  velocity  normal  to  a  helical  surfeice  with 

smgle  0^  at  a  radius  r  is 

u  =  -u  cos  0 .  +  u.  sin  0. 

XI  GL  jL  V  JL 

in  this  case  can  be  written 
g 


%=  I 

P^ 


sin  ®  cos  3^+5  cos  ®  sin 


(3.10) 
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CHAPTER  4 


SOLUTION  OF  PROPELLER  LIFTING-LlNE  PROBLEMS  BY  VORTEX  LATTICE  METHODS 


Introduction 

m' 

Before  applying  vortex  lattice  methods  to  the  solvition  of  pro¬ 
peller  lifting  surface  problems,,  it  would  seem  adviseible  to  apply  similar 
methods  to  certain  lifting  line  problems  whose  solutions  are  well  known. 

In  particular,  this  would  provide  some  preliminary  information  on  the 
spacing  and  arreingement  of  control  points  necessary  to  produce  results 
with  sufficient  accuracy  for  design  applications.  As  will  be  shown  in 
Chapter  6,  it  is  also  necessary  in  the  lifting  surface  case  to  have 
lifting-line  results  obtained  with  aa  identical  radial  lattice  arrangement. 
The  two  problems  which  will  be  discussed  are: 

1.  To  find  the  radial  distribution  of  circulation  to 
produce  a  free  vortex  sheet  of  true  helical  shape  in 
homogeneous  flow,  i.e.,  the  optimum  propeller. 

2.  To  find  the  radial  distribution  of  circulation  to  produce 
a  free  vortex  sheet  with  a  specified  radial  pitch  dis¬ 
tribution,  in  an  axially  symmetric  velocity  field. 

Goldstein  Factors 

The  solution  of  the  first  problem  is  e^ressed  in  terms  of 
Goldstein  Factors  which  are  defined  as  follows: 


H  (r,  X^, 


(^.1) 


where:  h  =  Goldstein  factor  (non-dimensional) 

* 

p 

r  =  strength  of  bound  vortex  at  radius  r  -  (ft  /sec) 
r  =  Raditis  of  vortex  element  tinder  consideration,  (ft.) 
u.^  =  tangential  component  of  induced  velocity  at  the 
lifting  line  as  shown  in  Fig.  4.1  (ft/sec) 


g  =  number  of  blades 

=  r/R  tan  =  X  "tan  3^^ 

=  angle  of  relative  flow  at  the  lifting  line 
X  =  non-dimensional  radius  r/R,  where  R  is  the  (ft) 
radiu®  of  the  propeller. 

[phis  prohlW  was  first  solved  by  Goldstein^^^  in  1929*  If  the  contraction 

and  axial  deformation  of  the  fre«  vortex  system  is  neglected,  the  problem 

can  be  reduced  to  the  two-dimensional  problem  of  a  rigid  helical  surface 

moving  with  a  fictitious  displacement  velocity  2u*  as  shown  on  Fig.  4.1.* 

Groldstein's  original  paper  Included  numerical  resxilts  foas  two-bladed 

propellers  for  2  <  l/X^  <  10  and  for  four-bladed  propellers  for  l/x^  =  5 • 

C  27 )  (2d) 

Later  Kramer'  ’  and  Lock  and  Yeatman'  '  obtained  veLlues  for  propellers 
with  2-5  blades  over  the  same  range  of  X^.  These  were  recomputed  in 
195^  by  Tachmindji  and  Milam^^^^  by  a  more  accurate  method.  Goldstein 
Factors  for  g  =  2-6  and  1.5  <  l/x^  5  ^  vere  obtained  using  a  Univac 
computer  at  David  Taylor  Model  Basin,  and  those  resxalts  showed  that 
previous  values  could  be  off  by  as  much  as  Tachmindji  emd  Milam^^^^ 

and  McCormick '  extended  Goldstein's  theory  to  include  a  finite  pro¬ 
peller  hub,  however,  their  initial  assumptions  regarding  the  vsilue  of 
the  circulation  at  the  hub  are  not  the  same. 


*The  velocities  shown  in  the  figure  are  at  the  lifting  line.  At  a 
large  distance  downstream  the  induced  velocities  are  doubled,  hence, 
the  displacement  velocity  is  2u^  . 


LINE 


VELOCITY  COMPONENTS 
RELATIVE  TO  OBSERVER 
ON  LIFTING  LINE 


FIG.  4.  I  VELOCITY  DIAGRAM  -  OPT  IMUM 
LIFTING-LINE  PROPELLER 


Another  way  of  computing  Goldstein  Factors  is  the  induction 
factor  method  developed  by  Lerbs^^^.  In  this  method,  the  velocity 
induced  by  each  helical  vortex  line  forming  the  sheet  can  be  computed 
from  a  potential  as  discussed  in  Chapter  2.  The  velocity  induced  at  a 
point  on  the  lifting  line  by  the  entire  sheet  can  be  obtained  by  integrating 
over  the  radius.  The  resulting  singular  integral  can  be  solved  by  ex¬ 
panding  both  the  circulation  distributiofl  and  the  induction  factors  in 
a  Fourier  series  with  a  prescribed  number  of  terms.  Tlie  integral  is 

then  approximated  by  a  series  of  singular  integrals  of  the  Glauert  type 

(19) 

whose  value  is  known  from  wing  lifting  line  theory' 

To  obtain  Goldstein  factors  by  a  lattice  method  the  free  vortex 

sheet  is  replaced  by  a  finite  niimber  of  helical  line  Vortices  as  shown 

schematically  in  Fig.  4.2.  The  velocity  induced  at  a  point  on  the 

lifting  line  by  any  of  these  vortex  lines  could  be  computed  either 

(7 ) 

from  the  potential  given  by  Lerbs'' '  or  by  numerical  integration  as 
described  in  Chapter  2.  In  this  case  rJumerical  Integration  will  be 
used  since  this  can  easily  be  extended  to  the  lifting  surface  case, 

while  the  two-dimensional  potential  for  the  induction  factors  cannot. 

♦ 

By  computing  the  velocity  induced  by  each  element  of  the 
lattice^at  a  number  of  control  points  on  the  lifting  line,  a  set  of 
linear  equations  results  relating  the  strength  of  the  individueil  vortices 
to  the  resultant  slope  of  the  flow  at  the  control  points.  This  can  be 
considered  as  another  way  of  getting  around  the  singular  integral  which 
occurs  with  the  continuous  vortex  sheet.  The  equivaO-ent  step  in  the 
induction  factor  method  is  determining  the  Fourier  coefficients  of  the 
Induction  factors  which  are  obtained  by  one  of  the  tisuaii  methods  of 
harmonic  einalysis  from  the  induction  factors  evaluated  at  a  number  of 
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FI6.  4.2  SCHEMATIC  ARRA*NGEMENT  OF  VORTEX 
LATTICE  WITH  M=  10 ,  P  =  5 


distinct  points.  In  general,  tine  velocity  induced  at  some  point  on  a 
propeller  blade  will  be  due  to  both  the  free  vortex  system  and  the  bound 
vortices.  However,  in  lifting  line  theory  where  the  blades  have  been 
replaced  by  straight,  radial  bound  vortices  only  the  free  vortex  system 
need  be  considered.  Shis  is  because  the  resultant  velocity  induced 
anywhere  on  one  lifting  line  by  a  symmetrically  arranged  set  of  lifting 
lines  of  equal  strength  is  zero.  . 

To  proceed  with  the  specific  formulation  of  the  problem,  it 

is  first  assumed  that  the  strength  of  the  bound  vortex  representing 

each  blade  is  given  by  an  I  term  Fourier  sine  series 

I 

“  (p)  -  •  1  »i 

i  =  1 


where  G  is  the  non-dimensional  bound  vortex  strength  and  p  is  a  new 

•  •  * 
variable  which  is  zero  at  the  hub  radius  r^^  and  tt  at  the  tip.  The 

variables  p  and  x  are  related  by 


X  =  I  (-1  +  )^)  -  I  (  1  -  cos  p 


p  =  cos 


-1 


1  + 


\ 


(4-3) 


The  vortex  distribution  given  by  (4.2)  is  automatlceJJLy  zero 
at  the  hub  and  tip  for  any  values  of  the  coefficients  a^.  This  is  in 
accordance  with  the  assumption  made  by  Lerbs^^^  and  Tachmindjl  and  Milam 
that  the  circulation  falls  continuously  to  zero  at  the  hub.  However,  as 

f  •31) 

indicated  by  McCormick'"^  '  and  a  recent  unpublished  study  by  Tachmindjl, 


■*<^Thi8  is  not  the  ususil  non-dimensional  circulation  which  is  defined  as 
G'  =  r‘/2nRV  when  V  is  the  speed  of  advance.  In  the  present  work^  it 
is  more  ccal^enient  °’to  use  u-*-  as  the  non-dimensionalizing  velocity  so 
that  G  will  be  independent  of  loading. 


-4l- 

the  assiuraption  of  zero  circulation  at  the  hub  does  not  appear  to  be 
valid,  but  rather  that  the  value  at  the  hub  should  follow  from  the 
solution  of  the  boundary  value  problem. 

In  any  event,  to  take  tR©  Irab  into  account  using  vortex 
lattice  methods,  it  would  still  be  necessary  to  obtain  a  sudbable 
series  ejfpansion  for  the  hub  potential  whose  coefficients  along  with 
those  in  (4.2)  could  be  obtained  by  including  control  points  on  the 
hub  cylinder  as  well  as  on  the  blade.  However,  sinee  the  effect  of  normal 
size  hubs  <  .2)  on  overall  propeller  performeince  is  small,  the 
solution  for  the  hub  potential  will  be  considered  at  a  later  time.  In 
the  megintime,  the  hub  will  be  takea  into  account  only  by  requiring  that 
G  (xj^)  =  0  while  the  radial  velocity  boundary  condition  will  be  dis¬ 
regarded.  As  will  be  shown  later  in  the  numerical  examples,  the 
Goldstein  Factors  obtained  under  these  fairly  crude  assumptions  are  in 
reasonable  aigreement  with  the  values  givea  by  Tachmlndjl  and  Milam^^^^ . 

The  vortex  lattice  arrangement  is  shown  schematically  in 
Fig.  4.2,  while  the  actual  arreuigements  used  in  the  numerical  examples 
are  shown  in  Fig.  4.3.  The  Interval,  from  r  =  r^^  to  r  =  R  is  divided 
into  M  equail  spaces  and  the  radius  to  the  inner  end  of  the  m'th  space 
is  called  r^.  The  continuous  bound  vortex  distribution  G(r)  is  replaced 
by  a  stepped  distribution  whose  value  is  equal  to  that  of  the  continuous 
distribution  at  the  mid-point  of  each  interval. 


'  'Cl  *  (■'oUj 

(1  <  m  <  M  -  1) 

‘’l  •  Cl  *  '■h'] 

(m  =  1) 

G,,  =  G[i  CB  +  (r^y] 

(m  =  M) 

(4.4) 
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FIG.  4,3  LATTICE  ARRANGEMENTS  USED  IN 
NUMERICAL  EXAMPLES  * 
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O' 


i 


The  free  vortex  lines  originate  at  (r  )  where  the  value  of 

o  m 

changes.  Calling  the  free  vortex  at  Gjjj  there  follows 


G  =  G  -  G  , 
m  in  m-1 


t4.5) 

This  can  he  made  to  hold  for  m  =  1,  2,  . M  +  1  by  defining  the 

non-existent  vortex  segments 


G  =  G„^,  =  0 
0  M+1 


(4.6) 


It  should  be  noted  that  the  same  result  could  be  obtained  by  noting 
that  the  strength  of  the  continuous  free  vortex  sheet  at  a  radius  r  is 
dG/dr  and  replacing  the  derivative  of  G  by  the  first  order  central 
difference. 

Ihe  free  vortex  lines  can  be  considered  as  replacing  a  con¬ 
tinuous  vortex  sheet  which  extends  l/2  space  on  either  side  of  the  free 
vortex.  The  only  exception  is  at  both  ends,  vrhere  in  the  continuous 
case,  the  sheet  must  end  at  the  hub  and  blade  tip.  It-would  therefore 
seem  reasonable  to  move  tlie  end  vortices  in  l/8  space  so  that  they 
vrauld  be  located  approximately  in  the  region  which  .viould.  actually  ' 
be  occupied  by  the  sheet.  In  this  case,  the  free  vortices  are  at  the 
following  radii: 

(r  )  =  r.  +  '  h'  ' 


o'm 


(“■0)1 


2  <  m  <  M 


M 


t  1/8  -  ''h) 


M 


1  -  va 


(4.7) 


M 


The  velocity  is  to  be  computed  at  P  control  points  located  at  radii  r^, 

r_  . r  midway  between  free  vortex  elements.  There  is  no  restrictior 

on  how  many  of  the  availiible  control  point  positions  are  to  be  used. 


4 
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'.Uhe  non-d^imensional  velocity  con5>onents  induced  at  by  a  set  of  semi¬ 
infinite  belical  vortices  originating  from  each  blade  with  radius  r 


om 


are 


(u  )  =  (u  ]l 

'  a'n^)  '  a'mp 


(u+L  =  (u.  ) 
t  'mp  t  'mp 


(u  )  =  (u  ) 

'  n'lirp  '  n'mp 


4tt  r 


-E  = 


a' 


mp 


m 


u*  G 


m 


mp 


m 


u*  G 


m 


kn  r 


m 


^  ^^n^inp 

Om 


(Jt.8) 


where  u  is  the  non-dimensional  velocity  as  defined  in  Chapter  2,  u  is 
the  dimensional  velocity  and  the  subscripts  a,  t,  and  n  denote  aocial, 
tangential  and  normal  congjonents . 

The  recpairement  that  the  relative  flow  at  the  lifting  line 
be  of  constant  pitch  can  be  seen  from  Fig.  4*1,  tc  be 


u. 


u 


u 


u*  = 


sin  cos 


n 


cos  P , 


cos 


=  const 


(4.9) 


ejqjressed  in  terms  of  either  the  tangential,  social,  or  normal  components. 
These  relations  make  use  of  the  known  result  that  the  resultant  induced 
velocity  is  normal  to  the  helical  surface  forme(^  by  the  free  vortex  system. 

The  tangential  velocity  induced  at  Xp  "the  set  of  vortices  is 


Mfl 


/  \  u* 

‘^t^p  “  2  Xp 


u* 


m' 

Mfl 

V 


4*  )T»m  ^TTI 

o  mp  m 


.  I 


2  Xp  A  i4  "  Vl) 


-45- 


*:  u*  sin  3.  cos  3.  (4.10) 

which  follows  from  (4-. 2),  (4.5 )>  (4-. 8),  and  (4-. 9).  The  subscript  i  in 

3^  following  gsneraUy*  accepted  propeller  nomenclature  stands  for  "induced 

angle"  euid  is  not  to  be  confused  with  the  index  i  in  the  Fourier  series. 

Rearranging  (4-. 10)  euid  cancelling  out  u*  gives 
I  M+1 

)  a.  )  (u.  ) _  (sin  ip  -  sin  ip  ,) 

i  Zj  '■  t'mp  ^  ^m  m-1' 

i=l  m=l  m 


=  2Xp  sin  3^^  cos  (4.11) 

Substituting  the  geometrical  relations 

Xj  =  X_  tan  3j  sin  3^  =  ^i 

i  p  ip  ip  - 

2 

V  ^ 

cos  3^^  =  _  (4.12) 

/  2  ^  ^2 

V 

into  (4.11)  gives  the  set  of  linear  equations  for  the  luiknown  coefficients 
®i 

I  M+1 

/  a.  )  (u,  )  (sin  ip  -  sin  ip,) 

L  1  L  ^  t'mp  '  m  '^m-l*^ 

i=l  m=l 

»  ('*■13) 

p  =  1,  2/ . I 


By  selecting  I  control  points  as  indicated  above  a  set  of  I  equations 
for  the  unknown  coefficients  results.  The  Gtoldstein  factor  at  any 
radius  can  then  be  determined  in  terms  of  the  a's  from  (4.1)  and  (4.2) 

=  ^  I  a,  .in  IP 

2X^ 


i‘=  1 


K 


(4.14)“ 
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Since  the  induced  velocity  components  are  a31  related  hy  (4.9)> 
the  set  of  equations  for  a^  can  he  e:q)ressed  in  terms  of  the  axial 
component  * 


(sih  ip  -  sin  ip  , ) 
'  '^m  m-1 ' 


(4.15) 


or  in  terms  of  the  ifiormal  componeftt 
I  M+1 

I  “l  I  ^  "m  "  ^  ‘’m-l) 

i=l  m=l 


(4.16) 


Numerical  Examples 

Since  the  integral  for  the  axial  velocity  is  the  easiest  to 
compute,  equation  (4.15)  would  he  the  most  efficient.  However,  to  test 
the  computation  scheme  for  the  normal  coii5)onent  which  would  he  needed 
later  in  the  lifting  surface  case,  a  program  using  equation  (4.l6)  was 
also  prepared.  The  greatest  discrepeincy  between  the  results  using  the 
axial  eind  normal  velocity  was  found  to  he  .0001.  The  method  of  com¬ 
putation  is  discussed  in  Appendix  (a  ). 

U 

Figure  4.4  shows  the  GrO*.stein  Factors  for  3-hladed  propellers 

with  zero  huh  diameter  hy  a  lattice  arrangement  with  M  =  24  and  P  =  8 

shown  schematically  in  Fig.  4.3.  The  curves  shown  in  solid  lines  are 

(29) 

taJcen  from  Taciimindji  and  Milam'  '  while  the  points  and  dotted  lines 
(where  necessary)  are  the  values  obtained  from  the  lattice.  Fig.  4.5 
shows  a  comparison  of  five  different  lattice  arrangements  in  the  case 
where  g  -=  3  aiid  =  .5  which  is  the  value  of  X^  which  showed  the 
greatest  disagreement  with  existing  data.  Each  of  the  lattice  arrangements 
are  shown  in  Fig.  4.3. 


GOLDSTEIN  FACTOR 


0.6  0.7  0.8 

NON-DIMENSIONAL  RADIUS  X  ‘ 

FIG.  4.5  COMPARISON  OF  GOLDSTEIN  FACTORS  OBTAINED 
BY  SEVERAL  LATTICE  ARRANGEMENTS  WITH 
VALUES  FROM  DTMB  REPORT  1034 


It  is  evident  that  the  lattice  ^results  are  in  agreement  with 
existing  data  both  for  low  and  high  values  of  In  the  region  where 

the  agreement  is  not  as  go®d,  extreme  variations  in  lattice  arrangements 
produce  cheingea  of  no  more  than  .003,  while  the  basic  disagreement  with 
(29)  is  about  .010. 

A  possible  explanation  of  the  discrepancy  may  lie  in  the 
method  of  computation  of  the  Goldstein  Factors  in  (29).  The  solution 
of  the  potential  problem  involves  the  solution  of  an  infinite  system 
of  linear  equations  relating  the  coefficients  in  the  series  e;q)ansions 
of  the  potential  outside  and  inside  the  propeller  radius.  For  small 
values  of  X^,  an  approximate  solution  to  the  set  of  equations  may  be 
expressed  in  closed  form.  For  large  values  of  X^,  this  approximation 
is  not  sufficiently  accurate,  and  a  more  exact  solution  was  developed 
by  Tachmindji  and  Milam  for  values  of  X^  >  .667.  For  values  of  X^^  <  .4 
the  approximate  coefficients  were  used,  and  the  range  in  between  from 
.4  <  Xj^  <  .667  were  obtained  by  interpolation. 

Since  the  only  noticeable  disagreement  exists  in  the  in-between 
region,  it  would  seem  likely  that  the  lattice  values  are  more  ewcurate 
in  that  interval. 

0 

As  an  additional  check,  calculations  were  made  for  6-bladed 
propellers  where  the  approximate  coefficients  were  known  to  be  much 
more  accurate  than  for  3“hladed  propellers.  The  results  are  shown  in 
Fig.  4.6  for  Xj^  =  .2,.  .4,  and  .667  and  it  can  be  seen  that  the  agreement 
is  very  satisfactory. 

As  was  mentioned  previously  in  the  discussion  of  the  hub 
boundary  condition,  Goldstein  Factors  were  aalculated  for  g  =  3,  =  *2^ 
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FIG.  4.6  GOLDSTEIN  FACTORS  6  BLADED  PROPELLER  ZERO  HUB 
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and  =  0.2.  These  are  shown  in  Fig.  h.6  where  4t  can  be  seen  that  the 

consequences  of  neglecting  the  boundary  conaition  of  zero  radial  velocity 

at  th®  hub  are  not  too  serious. 

Finally,  since  fewo-bladed  propellers  were  not  included  in 

recent  re-calculations  of  Goldstein  Factors,  a  complete  set  was  obtained 

By  the  lattice  method  and  the  results  appear  in  Fig.  ^.8.  Shown  on 

the  same  plot  are  some  values  tedcen  from  Lock  and  Yeatman^^^^  which 

seem  to  be  in  reasonably  good  agreement  with  the  new  data.  These  results 

also  seem  to  agree  very  closely  with  results  appearing  in  Goldstein's 

(6) 

original  paper  \ 

Non-Optimum  or  Wake-Adapted  Propellers 

The  preceding  development  can  be  extended  very  easily  to 
the  case  ■vrtiere  the  pitch  of  the  free  vortex  system  is  arbitrary,  and 
the  axial  inflow  velocity  is  a  prescribed  function  of  radius.  It  is 
assumed  that  the  pitch  angle  of  the  free  vortess  system  (r)  and  the 
geometrical  inflow  angle  P(r)  =  tan’^  (V  /(ur)'is  known  and  that  the 
non-dimensional  circulation  G  is  to  be  determined.  In  this  case  it 
will  be  necessary  to  compute  the  normal  velocity  con^jonent,  since  the 
resultant  velocity  is  not  necessarily  normal  to  the  free  vortex  sheets. 

In  this  case  the  boundary  condition  may  be  written  as  follows. 

M+1  I 

(u  )  =  ^ —  y  (u  )  y  a.  (u*  sin  ip  -  u*  ^  sin  i  t  ) 

'  n'p  2y  L  '  n'^fip  Zj  i  '  m  m  m-1  fc-1' 

^  m=l  i=l 

=  u*  (cos  (4.1?) 

In  this  case  u*  is  a  function  of  radius 
u*  =  u)r  (tan  0^  »  tan  0) 


(4.18) 


GOLDSTEIN  FACTOR 
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HYDROOYNAMIC  ADVANCE  COEFFICIENT  X| 


F’IG.4.8  GOLDSTEIN  FACTORS* 
G  =  2  BLADES 
ZERO  HUB 


as  can  be  seen  frcan  Fig.  4.1.  Introducing  the  ratio 


-  u*  .-("tan  3.)  -  (tan  P)  r 

C__  =  m  =  I  i'm  'mi 

^  Iir-  L(tant3J  -(fanPTJ  “ 

p  i'p  '•  'p 


m 

_  r"” 
P  P 


(4.19) 


into  (4.17)  gives  the  result 
I  M+1 

y  a.  y  (u  )  (  C  sin  ip  -  C  sin  ip.) 

1  Zj  '  n'mp  '  rap  m  m-l,p  m-1' 

i=l  m=l 

=  2Xp  cos  (P'^)p  (4.20) 

For  an  optimum  propeller  in  homogeneous  flow 
-1. 

2  2 

and  (cos  P^)^  =  Xp/  'Xp  + 

so  that  (4.20)  reduces  in  that  case  to  (,4.l6). 

The  program  prepared  for  the  computation  of  Goldstein  Factors 

was  modified  to  accept  an  aribtrary  distribution  of  P  and  P^,  and  the 

results  were  found  to  be  in  agreement  with  the  standard  induction  factor 

(^2) 

method  in  use  at  the  David  Taylor  Model  Basin  ' ,  except  near  the  hub 
■vdiere  the  hub  boundary  conditions  are  not  the  same. 
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CHAPTKR  5 

LIFTIHO-SURFACl  SOLUTIOIS  K)R  BLADES  07  ARBITRAlg  SHAHB 


Introduction 

In  this  chapter  we  consider  the  problem  of  determining  the 
camber  and  pitch  correction  for  a  propeller  vlth  a  prescribed  blade 
outline,  mean  line  r8.dlal  load  distribution.  As  indicated 

In  Chapter  1,  the  pitch  and  cainber  corrections  arst  determined  by  the 
requirement  that  the  prescribed  radial  load,  distribution  be  obtained 
vlth  the  sections  operating  at  their  Idee^.  angle  of  attack.  The 
chordvlse  load  distribution  Is  unknown  Initially  and  will  be  detenalned 
along  with  the  pitch  and  caid^er. 

The  nomenclature  used  In  this  chapter  Is  basically  the  same 
as  in  the  lifting  line  case  except  that  an  extra  dimension  must  be 
added  due  to  the  chordvlse  load  distribution.  As  shown  In  Figures  $.1 
and  5.2,  an  ('x',  y^,  z' )  cartesian  coordinate  system  Is  fixed  on  the 
propeller  with  the  x*  axis  axial  and  the  y*  axis  passing  throue^  the 
tip  of  the  Index  blade.  Ihe  z*  axis  coo^letes  the  rig^t-handed  system. 
A  cylindrical  system  (x',  r',  0)  corresponds  to  the  (x',  y',  z')  system 
with  e  B  0  on  the  y'  axis  and  positive  0  clockwise  when  looking  In  the 
positive  x'  direction. 

A  movable  cartesian  system  (x,  y,  z)  and  a  correspondixig 
cylindrical  system  (x,  r,  9)  Is  oriented  with  the  x  axis  axial  and 
the  y  axis  (or  (p  »  0  line)  passing  through  a  particular  control  point 
on  the  Index  blade. 

lliere  are  F  x  Q  control  points  on  the  index  blade  vhere 
p  B  1,  2,  . . . .  P  Indicates  the  radial  position  and  q  b  1^  2,  ....  Q 
Indicates  the  chordvlse  position*  It  should  be  mentioned  that  all  pairs 


vr 


FIG.  5.1  UFTIN6  SURFACE  NOTATION 


of  coordinates  or  subscripts  referring  to  radial  and  ohordvise  directions 
are  given  adjacent  alphal|R»tic  sysibols  with  tJie  higher  synbol  (alphabetically) 
z*ef erring  to  the  chordvlse  direction. 

Qiere  will  be  P  x  Q  gosslble  positions  for  the  novable  system 
and  the  nothtlon  y  ,  for  exaiq>le,  means  the  y  axis  of  the  movable  system 

It  » 

corresponding  to  the  pq'th  control  point.  FoUovlng  this  notation,  the 
(]uantlties  6^^  and  (^Q)pq  displacements  of  the  movable  system 

measured  frcm  the  fixed  system. 

A  non-dimensional  radius  is  defined  as  x  *  r/R  vhere  R  ’Is  the 
radius  of  the  propeller.  To  distinguish  the  radius  of  a  control  point 
from  that  of  a  helical  vortex  line  (on  the  end  of  a  bound  vortex  segment) 
the  latter  is  given  a  zero  subscript,  dbe  xion- dimensional  qiutntlties 
Tt  -  r^/r  and  ^  »  x/r  as  defined  In  Chapter  2,  will  also  be  used. 

Finally,  a  curvilinear  system  is  defined  at  any  radius  by  the 
intersection  of  an  axial  cylinder  with  the  reference  helical  surface. 

The  origin  is  taken  at  the  mid-chord  line  of  the  blade  whose  angular 
coordinate  in  the  (x*,  r',  6)  system  is  9.  Jibe  s  axis  is  along  the  helix 
with  the  positive  direction  towards  the  trailing  edge.  The  n  axis  is 
pezpendlcular  to  s  and  lies  on  the  gyllndrical  surface  with  positive 
direction  upstream  as  shown  in  Fig.  5*2.  If  the  cylindrical  surface  is 
expanded  and  viewed  from  the  propeller  axis  out  towards  the  tip,  a  blade 
section  results  as  shown  In  Fig.  5.3.  Ibe  chord  length  of  the  e:q>anded 
section  is  i(r),  consequently^  s  »  -i/2  corresponds  to  the  leading  edge 
and  s  s  f//2,  the  trailing  edge.  The  angle  of  attack  of  the  section 
relative  to  the  reference  helix  Is  a  and  the  maxi  mam  camber  measured 
from  the  nose-tall  line  is  given  the  symbol  f . 


LEADIN9 

TRAILINt 

CD«E 

EDOE 

FIG.  5.3 

EXPANDED  BLADE  SECTION 
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Ihc  Reference  Helix 

As  stated  In  Chapter  the  blade  surface  is  Msuaed  to  be 
approximately  on  a  helical  surface  vhose  pitch  at  any  given  radius  Is 
detezmlned  by  the  angle  of  relative  flov  according  to  lifting  line 
theory  vlth  the  same  radial  load  distribution.  However,  this  does  not 
define  the  surface  coag>letely  since  so  far  nothing  has  been  said  about 
the  relative  orientation  of  the  helical  lines  forsilng  the  surface. 

Since  actueil  propellers  may  have  both  rake  and  skew,  an  a^urate  definition 
of  the  blade  surface  is  a  fairly  disagreeable  geometrical  problem.  It 
Is  also  possible  that  the  effects  of  soaie  geometrical  vaintlonJy4re  of 
the  saste  order  as  the  errors  Introduced  by  the  basic  assumptions,  such 
as  the  neglect  of  the  defonoatlon  of  the  vortex  sheets.  Consequently, 

In  the  present  work  It  will  be  assumed  that  the  reference  helix  passes 
throufi^  the  y*  axis.  If  the  helix  Is  of  constant  pitch,  any  radial 
line  will  be  contained  In  the  surface,  however,  this  will  obviously 
not  be  so  If  the  pitch. Is  a  function  of  rsidius.  In  the  latter  case 
It  IS  further  eissimied  that  the  bound  vortex  segsients  are  radial,  and 
that  the  axial  distance  between  a  control  point  and  a  vortex  element 
Is  the  same  as  If  the  helical  surface  were  of  constant  pitch  corresponding 
to  the  pitch  at  the  control  point  radius.  While  these  slnpllfying 
assxuptlons  are  not  essential  to  the  application  of  the  vortex  lattice 
method^  it  would  seem  that  a  more  exact  geometrical  treatment  could 
not  be  Justified  until  the  effect  of  the  principal  variables  have 
been  detennlned. 

Bound  Vortex  Distribution 

Ihe  bound  circulation  distributed  over  the  blade  surface  will 
be  esqpressed  by  a  trigonometric  series  in  the  variable#  p  and  a  idilch 
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are  related  to  r  aM  s  by 

r  -  I  (R  +  rj^)  -  I  (R  -  rj^)  eoa  p 

■  ■  -i/2  coe  a  (5»i) 

lAich  there  f oUove 

r  >  r^^  lAen  P  -  0,  r  >  R  when  p  »  rr 

8  ■-i/2  lAen  a  «  0,  s  ■  f/2  ^dien  a  ■  tt  (5*^ 

The  vortex  sheet  streogth  Y  can  be  comreirted  to  a  xion-dlmenslonal 
quantity  3  by  dividing  by  the  displacement  velocity  u*  as  defined  In 
the  preceding  chapter.  It  is  assumed  that  3  can  be  represented  by  a 
series  of  the  foxm 

'  I  I  J 

3  (0-  »)  =  ^  [  I  'io  ^  I  ♦  I  1 

1-1  i«l  >1 

c^j.sin  1  p  sin  j^crj  (5.3) 

The  second  part  Is  a  Fourier  sine  series  which  has  the  property  that 
3-0  along  the  edge  of  the  blade  for  any  values  of  the  constants  c^j. 

Ihe  first  tern  goes  to  zero  all  along  the  trailing  edge,  but  tends  to 
Infinity  at  the  leading  edge .  For  a  fixed  value  of  p  this  Is  the 
chordtflse  clrctilation  distribution  of  a  flat  plate  at  a  small  angle 
of  attack  In  two-dimensional  flow.  According  to  lineaorized  two-dimensional 
thin  airfoil  theory the  chordwlse  circulation  distribution  of  any 
mean-line  can  be  obtained  by  superimposing  the  flat  plate  distribution 
and  a  general  distribution  iAicb  Is  zero  at  both  the  leading  and  trailing 
edge.  The  an^e  of  attack  for  idilch  the  coefficient  of  the  "flat  plate” 
tern  Is  zero  it  called  the  Ideal  angle  of  attack. 
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Qm  radial  circulation  distribution  la  obtained  by  Integrating 
Y  over  the  chord  at  a  particular  radius 


r  (p) -fv  (s-i.) 

o 

An  tenw  of  non*»dlaensloxyQ  g^aatltle8 

TT 

»(P)-^  J  ^5.5) 

o 

tdiere  Q  Is  the  non-dimensional  circulation  defined  In  the  preceding 
chapter  as 


G  = 


ttDu* 


(5.6) 

II  IJLt 

Sof  A 

iteistltutlng  (^>3)  for  3  In  (5.5)  Integrating  gives  the  result* 

I 

9  (p)  »  ^  (2  c^^  +  Cjj^)  sin  1  p  (5.7) 


1«1 


If  ve  now  require  that  a  particular  radial  lo§d  dlstrlbutlcm  0  (p)  Is 

to  be  obtained  In  the  sections  operating  at  their  Ideal  angle  of  attack, 

there  follows  that  c.  *0  and  that  c.^  are  the  known  Iburler  coefficients 

lo  u.  * 

of  the  radial  clrcxilatlon  distribution.  The  remaining  coefficients 


r-1  »  1,  2,  ....  I-| 

°1J  Lj  =  2,  3,  ....  jJ 

irtilch  do  not  contribute  to  the  radlsd  load  distribution  are  to  be 

deteznined  by  the  boundary  conditions  on  the  blade  surface.  ?or  later 

.  use.  It  will  be  coxxvenlent  to  define 

I 

“I  ""ij  ^  ■  (5.8) 

1-1 


*nie  details  appear  In  several  aerodynamics  texts  such  as  "Ilieory  of 
Wing  Sections"  (33). 
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bo  that  (5*3)  becomeB 

J 

S  (p,  a)  bj  (p)  fin  J  <7  (5.9) 

provided  the  angle  of  attack  at  each  radius  Is  Ideal. 

Vortex  Lattice 

Ihe  contliuous  bound  vortex  sheet  Is  to  be  approximated  by 
a  finite  number  of  radial  bound  vortex  segments  ectch  with  constant 
strength.  At  the  ends  of  each  segment  a  free  vortex  of  the  same 
strength  must  be  shed  forming  a  "horseshoe"  voxi:ex  system  as  Shorn 
in  Figs.  ^.1  and  5.2.  Naturally >  parts  of  the  free  vortex  system 
originating  from  bound  vortices  at  the  same  and  Immediately  adjacent 
radii  coincide.  Although  this  fact  will  be  useful  for  computational 
purpo^(e8,  each  horseshoe  system  vlU  be  considered  logically  to  be  an 
Independent  \mlt. 

The  lattice  arrangement  Is  obtained  by  dividing  the  Interval 
between  the  hub  and  blade  tip  Into  M  equal  spaces.  Free  vortices  are 
shed  at  radii 

(ro)„  -  *  r,  (5.10) 

M 

except  at  the  ends,  where  they  are  moved  in  I/6  space  towards  the 

'A 

Interior  of  the  blade  (as  In  the  lifting  line  case).  There  are  N 
radial  vortex  elements  between  any  two  adjacent  values  of  r^.  Ibese 
will  be  centered  at 

f  fCr  )  +  (r  )  il  •  (5«11) 

m  2  L'  o'm  '  o'm  +  IJ  / 

and  will  be  located  by  dividing  the  chord  length  at  r^  Into  N  equal 

panels  with  the  bound  vortex  at  the  mld-polnt  of  each  panel.  The 
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cliprdihjie  position  relstlve  to  the  ■Id-chord  line  is  given  hy 

(  2n-  M  -  1)  (5.12) 

and  the  angular  coordinate  measured  clocfetfise  fron  the  y*  iusts  is 

e_  -  0_  +  “*■  ^i^m  (^  -  H  -  1)  (5.13) 

on  D  xsst  "■ 


Control  points  are  located  at  the  nl^^ints  of  the  panels 
fomed  hy  the  horseshoe  elements.  In  generalf  there  will  he  nany  more 
horseshoe  elesients  than  control  points,  and  it  is  coqpletely  axhltrary 
vhich  of  the  possible  control  point  arrangesMnts  are  to  be  used.  However, 
to  sliqillfy  the  coaqiutatlons  sosievhat,  it  vlH  be  assuaed  that  the 
chordvlse  arrangeoent  of  control  points  will  be  the  ssme  at  each  radial 
position  used.  The  niad>er  of  chordvlse  control  points  is  given  by  the 
expression 

I  -  2  C,  -  Cp 

Q - 2  i3.1h) 


where  is  the  nuaber  of  radial  vortex  eleamts  between  each  control 
point  and  is  the  number  of  uaused  control  point  positions  ^between 
the  leading  edge  and  the  first  control  point.  If  (5*1^)  is  a  fraction, 
only  the  Integer  part  is  to  be  retained.  Fig.  (^.4)  shows  a  nustber 
of  chordwlse  lattice  arrangasents  corresponding  to  various  values  of  if. 


Ihe  control  point  angles  are  then  given  by 


P<1 


®p 


Jtp  (cos  P^)p  [2CCj^  (q  -  1)  +  Cg  !)•  -  »] 


(5.15) 


Ibere  are.  a  total  of  P  radial  positions  used,  and  are  subject  only  to 
the  restriction  that  P  <  M.  The  total  number  of  control  points  is  P  x  Q. 

i/^ 

Relating  Conintuous  and  Lattice  Distributions 

Let  G _  be  the  non-dimensional  strength  of  the  bound  vortex 

JBQ 

located  at  0 _  and  centered  at  r  *  The  strengths  of  the  individual 

mn  n  « 


Q  y.  A 


elements  are  first  of  all  subject  to  the  requirement  that  the  radial 


load  distribution  be  the  same  as  In  the  continuous  case 

1 

naOL 

Ihe  remaining  I  -  1  requirements  will  be  that  the  lattice 
and  continuous  dlstrlbutloiis  induce  the  same  velocl'ty  at  each  of  the 
X  -  1  possible  control  po&ait  positions  In  2-dimenslonal  flov.  Iren 
thln-alrfoll  theorj'-  the  non-dimensional  velocity  Induced  at  thef  q'th 
control  point  by  the  1  vortices  at  a  particular  radius  r^  Is 


(u^)  in) 

u*  m 


N 


>  mn 

^=1  2  (n  -  q)  -1 


(5.17) 


\diere  u^  Is  the  dimensional  velocity  normal  to  the  vortex  sheet.  The 
velocity  induced  at  the  same  point  by  the  continuous  distribution  can 
be  shown  to  be 


J 

^  bj  cos  (j  +  1) 
J-0 


'  where  a  »  cos 
<1 

t 

Bquatlng  (5.1©)  and  (5.17) 
Is  obtained 


idilch  cGsdsined  vlth  (5.16) 

the  unlmown  0 _ . 

mi 


for  each  value  of  q  the  following  equation 
J 

•II  bjCO.  (J  *l)a^ 

*  J-0 

q  »  1,  2,  ...  II  -  1  (5.19) 

results  in  a  set  of  K  linear  equations  for 


Let  the  solution  of  this  set  of  e^stlons  he  expressed  In  the 


fora: 


an 


I 


(5.20) 


13ie  chord  load  factors  (i  .  are  constants  tdilch  can  he  conputed 


(17) 


and  hy  Van  Bom 


once  and  for  all.  Values  of  4  are  glren  hy  mUmer 
and  deYoung^  The  latter  values  are  slights  different^  the  author# 
stating  that  the  former  values  are  Ipcorrect.  Hovever,  on  re«calcttlating 
the  chord  load  factors,  it  vould  i^ear  that  TOUmer's  original  values 
are  correct  i!  Values  of  correct  to  6  decloal  places,  vere  re-coqputed 
for  M  •  2f  kf  6f  6  and  J  ■  0,  1,  2,  ....  X-1  using  an  IBM  65O  and  these 
results  appeal:  in  Appendix  (  C ) . 

* 

Velocity  Induced  hy  the  Lattice  In  3-Biinenslonal  flow 

Let  u _  he  .the  noxaal  cosmonent  of  the  non-dimensional 

nopq. 

velocity  induced  hy  the  conplete  horseshoe  system  0^  at  the  control 

* 

point  at  r  ,  0  .  The  subscript  n  for  "noxaal"  vlU  he  ooitted  in 

P 

this  section  slnct  only  the  noxaal  cosponent  will  he  considered.  As 
in  Chapter  2,  u  is  related  to  the  dimensional  velocity  u  hy 

(5.21) 


u 


ni^q 


u  hr  r 

flPipq  ] 

SSI 


vhlch  can  also  he  eigpressed  in  terms  of  the  non*dliienslanal  circulation 

(5-22) 


3  .UsssJJs 

nopq  u* 


8SX 


IMs  velocity  can  he  cooputed  hQr  a  procedure  vhlch  is  outlined  in 
Appendix  (  A)  using  the  results  of  Chapters  2  and  3. 


DeteiMinlng  tiie  Canber  and  Angle  of  Attoclc  ® 

Af  nftitloaed  preriously;  It  Is  SMimaBM.  tbat  the  l^^fide 
surface  is  to  be  formed  such  that  Its  esgpaoded  sections  mi^  all  be 
derived  from  a  single  mean-line  by  suitably  selecting  the  caober/length 
ratio  f/i  and  angle  of  attack  a  at  each  radius  t  33ie  angle  of  attack  Is 
to  be  measured  from  the  Induced  Inflow  angle  3^  detezmlned  from  lifting 
line  theory.  It  la  also  assuiaed  that  the  magnitude  of  the  resultant 
Inflow  velocltjr  Is  the  same  as  In  the  llftixig  line  case,  namely,  V*. 

The  value  of  '.f/innd  a  at  each  radius  are  determined  by  the  boundary 
condition  that  the  flow  be  tangent  to'  the  mean  line  at  each  control 
point.  The  slope  of  the  mean  line  relative  to  3^  at  a  particular 
chordwlse  station  is 

*  (5-23) 

\diere  h  is  the  slope  of  tlie  mean  line  with  unit  camber  ratio.  As  can 

<1 

be  seen  from  ?lg.  5*^,  the  boundmry  condition  can  be  written 

Op  - (f/i)p  -  Y  L  ^ 

assuming  that  the  Induced  angles  are  small.  Introducing  ($.22)  axid ' 
noting  that  3^  .  3  w  u*  cos  3^/V*,  there  follows 

M 

[uj  (CO.  ^ 

(5.25) 

It  is  now  convenient  to  ej^ress  vc^/V*  In  terms  of  the  lift  coefficient 


o  -  h  (f/i)  »  ^ 
P  q  '  'P  V  * 


Z  ^^mnpq  ®nnJ 


nml 
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of  the  section.  -  From.  Kutta-Joukov^kl’s 

dL  ■  P  V»  r  dr  (5*26) 

idiere  dL  Is  the  lift  force  acting  on  an  elesMnt  of  £ound  vortex  of  radius 
dr  and  P  is  the  fluid  auss  density*.  Ihe  lift  eoeffielent  Is 


®l"T 


dL 


|P  V**"  idr 


ST  /2ff  Q\/U*\ 

Ti*  ■ 


<5.27) 


Replacing  0  hy  (5*27)  and  conblnlng  with  (^.29)  and  (9*20),  *there 

foUovs 


i. 


(f/f)  -(i/D) 


M 


2TT 


[(CO*  ei)p  *  h;; 


V  J 


n^r 


Bmpq 


,^nj 


(5*26) 


vhere  C'  is  a  factor  ^Ich  takes  into  account  that  u*  may  be  a  function 

j^p  m  ^ 

of  radius  and  is  defined  by 


^  (tan  3^)^  -  (tan  3 

^mp  “  u*  *  (tan  3. )  -  (ian  P) 

Mr  *  Jkr  Mr 


0 


(5.29) 


For  optimum,  open  water  propellers,  u*  is  Independent  of  radiiis  so  that 
C  ■  1  and  may  be  oaiitted  in  (5.26). 

The  quantities  on  the  left  in  (^.26)  are  the  angle  of  attack 
and  caai>er  ratio  per  tialt  lift  coefficient  and  are  given  the  syid>ols 


a/c. 


(5.30) 


'L 

In  two-dimensional  flow,  these  are  constants  which  depend  only^  on  the 
type  of  mean  line.  The  zhtio  of  the  caid>er  required  in  three-dimensional 


*In  all  equations  except  (3*26)  and  (3*27)  the  syd>ol  p'  is  the  trans« 
foxmed  radial  coordinate. 
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to  that  required  for  an  equal  lift  coefficient  In  tvo-dlaenalomal 
flow  Is  the  camber  correction  factor  as  defined  In  current  propeller 
desl^  methods  However,  a  similar  definition  cannot  be  used 

for  the  pitch  correction  since  the  Idea^  angle  of  <&ttaeXf  of  many  mean 
lines  In  two-dimensional  flow  is  zero. 

Iqiuatlon  (5.26)  written  for  each  control  point  represents  a 
set  0^  linear  equations  for  a,  f  and  the  coefficients  of  the  non-llft- 
producing  part  of  the  circulation  distribution.  Rearranging  (5*^) 
to  put  the  unknowns  on  the  left  and  introducing  (5*8) 


M  K  I 


q  -  1,  2,  ...  Q  (5.31) 

If  the  nuniber  of  reidlal  tezms  1  In  the  Fourier  series  for  the  circulation  distri¬ 
bution  is  equeil  to  the  number  of  radial  control  point  positions  P,  and  If  the 
number  of  chordwlse  terms  J  is  one  less  than  Q,  the  number  of  unknowns  will  be 
2P  +  I  (J  -  1)  «  2P  +  P  (Q  -  2)  -  PQ  (5.32) 

which  equals  the  number  of  equations.  The  reason  that  J  s  q  -  1  is 
that  the  first  term  of  the  series  Is  determined  In  advance  by  specifying 
the  radial  load  distribution.  Consequently,  there  must  be  at  least  two 
chordwlse  control  points  in  order  to  determine  a -pitch  and  camber  correction. 

The  set  of  equations  represented  by  (5.31)  can  be  written  In 
matrix  notation 
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vhere 

r 


1 


•la  i 


. 


0^  .../«2p-l<2P-l 
•••'<-  2p  <  2P 

c  ...  i  -  ^  +  (1  -  l)(j  -  1)  +  j  -  1 


k  -  (p  -  1)  Q  +  q 


(5.33) 
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CHAPTER  6 

A  LiraiMG  SURFACE  30LUTI0H  FOR  PROFELURS 
WITH  SYMMETRICAL  BLADES 


The  Syaitetry  of  the  Velocity  Field 

In  th^t  case  \ihen  Both  the  blade  outline  and  the  mean 

line  are  syiimetrlcal  about  the  y'  axis^  an  Ingportant  sln^liflcatlon 

results  from  the  syianetry  of  the  integrals  determining  u  .  As  a 

nmpq 

result,  it  can  be  shown  that  within  the  limitations  of  the  assvunptlons 
outlined  in  Chapter  1,  a  propeller  with  syometrical  blades  has  no  pitch 
correction  due  to  lifting  snrfeuie  effect. 

First  of  all,  defining  9  as  the  angle  between  a  control 
point  and  a  radial  bound  vortex  or  an  elemeat  of  a  helical  vorti^x,  it 
is  evident  that  the  non-dimensional  normal  velocity  induced  hy  a  bound 
vortex  u^  is  an  odd  function  of  9,  while  the  nozmsQ.  velocity  induced 
by  an  element  of  hellceLL  vortex  6u^  is  an  even  function  of  9.  !Ihls  can 
be  seen  frc*  (3.9)  and  (3.10)  for  the  bound  vortices,  since  both  sin  9 
and  5  are  odd  functions  of  9.  The  fact  that  6Uj^  is  an  even  function 
of  9  can  be  deduced  from  (2.10)  and  (2.11). 

We  now  consider  the  velocity  Induced  at  three  synmetrically 
oriented  control  x>olnts  (labeled  L,  M  and  R)  as  sketched  in  Fig.  6.1. 

For  simplicity,  portions  of  three  horseshoe  vortex  elements  are  shown 

and  are  numbered  1,  2,  and  3  with  2  on  the  y*  axis  and  1  anA  3  syoBrntrically 

arremged  Vlth  respect  to  the  y'  axis. 

The  relative  strength  ot  the  n'th  bound  vortex  corresponding 
to  the  J'th  term  in  the  Fourier  sine  series  is  given  by  as  defined 
in  (5.19).  However,  it  is  'sufficient  to  note 'that  is  an  even 
function  of  n  anl  6  when  J  is  odd^  and  an’  Odd  function  of  n  and  6  idien 
J  is  even. 


FIG;  6,1  ILLUSTRATION  OF 
SAMPLE  CONTROL  POINTS  8 
VORTEX  LATTICE  ELEMENTS 
ON  A  SYMMETRICAL  BLADE 
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We  first  detezmlne  the  velocities  Induced  alf  M  by  the  vortices 
located  at  1  and  3  with  strengths  corresponding  to  ttie  first  term  In 
the  Fcarler  series^  \dilch  Is  the  only  term  contributing  to  the  radial 
load  distribution.  Since  the  strengths  of  1  and  4  are  e^pral  In  this 
case,  the  velocity  Induced  by  the  bound  vortices  cancels,  idille  the  tvo 
hellced  vortices  starting  at  1  and  3  are  equivalent  to  a  gingle  vortex 
of  twice  the  strength  starting  at  2.  Consequently,  the  velocity  at  M 
due  to  the  first  term  in  the  series  Is  the  same  as  In  lifting  line 
theory.  It  Is  also  evident  that  the  difference  between  the  velocity 
according  to  lifting  line  theory  and  the  velocity  induced  at  L  and  R 
Is  an  odd  function  of  9.  Therefore,  as  far  as  the  first  tens  In  the 
series  is  concerned,  the  mean  line  should  be  symuetrlcal  about  the 
mid-chord. 

lext  consider  the  even  terms  In  the  series,  J  4,  6  ... 

In  which  case  the  strengths  of  1  and  3  will  be  equal  and  opposite.  The 
velocity  induced  at  M  by  1  and  3  vUl  be  non-zero  since  the  effects  of 
1  and  3  vUl  add.  furthermore,  the  velocity  induced  at  L  and  R  will 
be  equal* 

finally,  we  consider  the  case  idien  ^  •  3)  T  •••  bo  that 
1  and  t  again  have  equal  strengths.  Using  the  same  synoetrf  argumants 

as  In  thd  case  of  J  «  1,  we  conclude  that  the  velocity  at  M  Is  the  same 

3 

as  If  1  emd  4  were  cooblned  and  located  at  2,  sad  that  the  difference 
between  the  velocity  ssssordlng  to  lifting  line  theory  and  the  velocity 
Induced  at  L  and  R  Is  an  odd  function  of  9.  However,  for  J  >  1  the 
total  strength  of  the  chordwlse  lattice  elements  must  be  zero  according 

a 

to  so  that  the  Induced  velocity  obtained  by  ccmjbining  all 

the  vortex  elements  at  2  must  be  zero.  Hence,  the  velocity  induced 
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at  M  li  zero,  and  the  velocitlea  Induced  at  L  and  R  are  equal  and  opposite. 
the  SlJBultaneoue  Equations 

We  next  eonsider  the  effect  of  this  synmetry  on  the  set  of 
equations  given  in  (5*32).  For  simplicity  it  will  he  assumed  that  P  >  1 
and  Q  B  however,  tl^  conclusions  will  be  vSHid  in.  ^Sjjgeneral«case^ . 
When  written  out,  the  equations  would  look  as  follows: 

•  *11  ®  *  “12  f  +  *13  “la  *  *13  *  *15  ■  '“1 

*U  “  *  *22  '  *  *23  *12  *  *21  *13  *  *25  *11  ■  ”2 

*11  “  *  °  *  *33  *12  *  °  *  *35  °ll  ■  ° 

*U  “  ■  *22  ^  *  *83  *12  ■  *2l  *13  *  *25  *ll  "  "”2 

*11  “  ■  *12  ^  ’  *13  *12  -  ‘11  *13  *  *15  *11  ■  •”1 

where  the  a's  and  b's  cure  elements  of  the  A  emd  B  matrices  respectively 
as  defined  in  (5.32).  The  unknowns  a  and  f  are  the  pitch  and  camber 
factors  defined  In  (5.29)  and  the  c's  are  the  unknown  coefficients  in 
the  circulation  distribution  defined  in  C5*3)«  'Oxe  symmetry  of  the 
coefficients  has  already  been  incorporated;  for  example,  a^^  has  been 
replaced  by  a^^ . 

Kliainatlng  c^j^  between  (6.1)  and  (6.2)  as  well  as  between 
(6.4)  and  (6.5),  a  reduced  set  of  equations  is  obtained 

d^  i  +  d^  f  +  d^3  c^  4  d^l^  c^3  -  e.  (6.6) 

^^a  +  0  +«33Ci2  +  0  -0  (6.7) 

d^  a  -  d^  f  4.  d^3  c^2  -  ^4  ®13  "  ■®' 

where  the  d's  and  e's  cure  related  to  the  a's  and  b's  as  follows 

^  “  ®11  “  ®11  “15/^25 
®1-  -  ^  '  ^2  *15/^25 


(6.1) 

(6.2) 

(6.3) 

(6.4) 

(6.5) 


The  uokQovuja,  ?  and  can  be  ellalnated  between  (6.6),  (6.7)  and 
(6.8)  to  give  the  following 

%  -  ‘33  '02  '  ° 

from  which  we  conclude  that  c^  mist  be  zero,  provided  the  constant  In 
parentheses  Is  non-zero.  However,  since  the  constant  Is  made  up  of 
Ind^endently  variable  geometrical  liqputs.  It  will  not  be  zero  In  general. 

Consequently,  It  can  be  seen  from  (6.7)  that  a  mast  also  be 
zero,  hence,  there  Is  no  pitch  correction.  Furthezmore,  It  Is  evident 
in  this  case  that  (6.6)  and  (6.8)  are  redundant. 

Equations  (6.1)  and  (6.5)  nay  now  be  re-wrltten  as  follows 

.  (6.11) 

(6.15) 


a,«  f  +  a,,.  c,„  +  a,^  c,,.  «  b. 


-12  -  ■  *14  ^13  "  “15  ^14  ‘'l 

-*oa  ^  -  *11.  '13  *  *15  'o'*  "  ■'’1 

showing  that  c^^^  0.  Following  the  same  proceddre.  It  cam  be  concluded 


that  c^j  mist  be  zero  for  all  even  values  of  J,  so  that  the  circulation 
distribution  must  be  an  even  function  of  6. 

By  removing  all  the  zero  terns  and  redundant  equations  from 
the  original  equations  (6.1)  -  (6.5),  the  following  equlwalent  set  of 
equations  is  obtained 


®12  ^  ■*■  ^14  *^13  “  ^1 
*22  ^  ^  *24  ‘'13  '  ^2 


(6.13)  • 


which  is  a  fairly  drastic  sliq>llficatlon. 

Modification  of  Preceding  Besults  for  Synnetrlcal  Blades 

The  development  in  Chapter  5*  vlU  ^w  be  modified  to  taka 
advantage  of  these  results.  iQie  continuous  vortex  sheet  strength  (3.3) 


Is  re-wrltten  as  follows: 
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I  J 

3(p,  ^  ^  ^  ^  ^  -  1)  ‘y  (6.14) 


irtiieh  Is  syanetrlcaLL  about  the  rntd-diord.  Control  points  vlll  be  dis¬ 
tributed  only  over  the  dovnstream  half  of  the  chord,  and  in  peurticular 
cannot  be  located  at  the  Md-chord,  since  this  will  result  in  A  being 
singular.  It  is  also  convenient  to  define  I  ca  the  nuniber  of  chordvise 
lattice  elesients  on  each  side  of  the  mid-chord,  so  that  the  total  nuaiber 
is  2K. 


The  angular  coordinates ’ of  the  bound  vortex  elements  are 

given  by  the  expression 

/  (cos  B . )  /.  \ 

^  1 _ -lis 

2”  % 


m 


(6.15) 


which  replaces  (5.13).  The  number  of  chofdwise  control  points  Q  is 
atm  given  by  (5.14)  since  J  has  been  re-defined.  However,  the  ex¬ 
pression  for  the  control  point  angles  (5*15)  follows 

I  (cos  3 .  )  r  -I 

(6.16) 

The  final  set  of  equations  is  practically  the  same  as  in  (5  >30), 
except  that  the  terms  containing  the  pitch  correction  a  are  no  longer 
present . 


Jwt  (b, )  V  h 


M 


2ir 


Ll  ^  Ll  ^.2  '‘"J 


sin  1  p 


M  2H  I 

2)(p  (CO.  6,)p  .  .in  1 


(6.17) 
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Finally,  the  location  of  the  matrix  elements  corresponding  to  (^•32) 
is  as  follows: 


f^k  »  4p  -  1)  Q  +  q 

ri/D)p 

i  -P 

1  M  2ir 

r*  - 

L  ( 

k  ■  (p-1)  Q  q 

1  .  P  +  (1-1)  (J-1)*  J-1 

M  SX 

I 

®K  "  %  ^  2l  i^ 

^  Pm 

k  -  (p-1)  Q  +  q 


?...  /■=p<P. 

Ci^  ...  i  -  P  +  (1-1)(J-1)  +  J-1  (6.1fi) 

There  is  one  Important  consideration  in  using  the  sli^llfled 
set  of  equations  given  in  (6.17).  In  the  case  of  Chapter  the  pitch 
angle  of  the  free  voztex  system  for  a  prescribed  radial  circulation 
distribution  .did  not  have  to  be  given  exactly,  since  small  errors  in  0^ 
could  be  absorbed  in  the  pitch  correction.  However,  In  this  case  any 
dlscrepamcy  between  G  and  0^  will  come  out  as  an  error  in  the  caoiber 
correction,  since  the  assumed  symmetry  will  not  actually  be  present.' 

A  slsg>le  way  to  avoid  this  difficulty  is  to  obtain  the  relationship 
between  0  and  0^  by  the  method  discussed  in  Chapter  4,  using  precisely 
the  sane  radial  lattice  arrangement  as  in  the  lifting  surface  case. 

Ohls  also  happens  to  be  convenient  since  the  Fourier  coefficients  of 
the  circulation  dlstrl'butlon  c^  are  obtained  directly  in  the  lattice 
solution  of  the  lifting  line  problem. 


This  procedure  vas  incorporated  in  the  con5>utation  scheme 
\<hich  is  outlined  in  Appendix  (A  ).  The  resulting  camber  correction 
feictors  are  shovm  in  Chapter  7^  together  with  the  results  for  asymmetrical 
blades  using  the  results  of  Chapter  3* 


CHAPTBR  7 

RESULTS  ARI)  COHCLUSIOIS 


Analyals  of  Lifting  Surface  Results 

There  are  two  prlxiclpal  questions  vhlch  need  to  be  ansvered 
In  determining  the  effectiveness  of  the  vortex  lattice  method.  First 
of  all^  it  is  in^ortant  to  detexmlne  hov  fine  a  lattice  speiclng  Is 
necessary  to  produce  results  with  the  desired  accuracy.  Obviously^ 
the  method  would  be  of  little  practical  value  If  the  required  spacing 
were  so  small  that  unreasonably  long  ccoputatlon  tines  were  needed.  »In 
addition;  extremely  small  spaclngs  would  require  special  measures  to 
avoid  the  loss  of  significant  figures  which  would  also  Increase  the 
cooQjutatlon  time.  * 

The  second  question  Is  whether  the  formulation  of  the  llftlng- 
surfaoe  problem  with  the  sliiq>llfylng  assus^tlons  introduced  in  Chapter  1 
is  an  adequate  representation  of  the  physical  situation. 

Considering  the  first  question,  the  convergence  of  the  lattice 
approximation  in  a  typical  case  was  studied  by  coagwtlng  camber  corrections 
using  six  different  lattice  spaclngs.  Gbe  characteristics  of  the  propeller 
and  the  lattice  parameters  are  given  in  Table  J.l.  The  blade  outline, 
in  this  case,  was  sysmetricej,  and  corresponded  to  the  Troost  B-8eries  (3$)* 

Table  7-1 

Data  for  Test  Calculations 

Propeller  Data 

Number  of  blades,  g  =  3  . 

Expanded  Area  Ratio  A^/A  ^0.6^ 

Mean  line  type  -  Parani^Ic 

Inflow  velocity  >  constant  (open  water) 

Circulation  distribution  >  optimum 
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Lattlce  Pfcruietere 

*  Teat  12  3  ^  5  6 
Badlal  lAttlee  •pace*,  N  6  8  Q  2k  ^  2k  2k 
Cbordvlae  lattice  apaeea,  2B  k  6  8  k  6  6 
Radial  control  polnta^  P  3  3  3  ^^^ 
Chordvlae  ocotrol  polnta,  Q  12  3  12  3 
^o^patatlon  tljae  (alnatea}  •  IBN  709  2  3  3  3  11  17 


The  Initial  reaulta  of  thta  teat  vere  fairly  erratic,  particularly 
near  the  tip  of  the  blade .  The  reaeon  for  thla  waa  that  too  nany  teraa 
la  the  Fourier  aerlea  for  the  circulation  dlatrlbutlon  vere  retained, 

X 

aa  can  be  aeen  froai  the  following  conalderatlona . 

The  nuaerlcal  reaulta  Indicated  that  the  nomal  yelocity  ccai- 

% 

poneat  Induced  by  the  kaova  part  of  the  circulation  dlatrlbutlon 
I 

aln  1  p  aln  a  (7«l) 

waa  alaoat  a  linear  function  of  the  chordvlae  dlatance  a.  It  waa  alao 

*  f* 

noted  that  the  Induced  velocity  flelda  obtained  fron  each  of  the  lattice 
arraageaenta  vere  in  good  agreenent,  the  'only  noticeable  dlfferencea 
occurring  with  the  largeat  apeclng  uaed.  Conaeguently,  the  erratic 
reaulta  could  only  be  due  to  the  way  In  which  the  hl^ier  coefflclenta 
In  the  circulation  dlatrlbutlon  vere  detexslned. 

Since  a  parabolic  nean  line  vaa  uaed  In  theae  exanplea,  the 
hltfier  tezna  In  the  Fourier  aerlea  would  be  zero  if  the  velocity  Indeed 
by  (7.1)  were  exactly  a  linear  flmctlon  of  a,  at  which  point  the  chord¬ 
vlae  load  dlatrlbutlon  would  be  the  aaae  oa  In  two-dlaMnalonal  flow. 

However,  In  thla  caae  additional  texiu  axe  required  alnce  the  velocity 

•  '' 

Induced  by  (7*1)  la  not  exactly  a  linear  function  of  a.  Theae  hl«^r 
texaa  Induce  velocity  flelda  which  vary  acre  or  leaa  alnaaoldally  over 
the  chotd*  Since  the  coefflclenta  of  theae  texaa  are  detcialned  only  bp 


CAMBER  CORRECTION 
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FIG.  7.1 


COMPARISON  OF  CAMBER  CORRECTIONS  • 
OBTAINED  WITH  SEVERAL  DIFFERENT 
LATTICE  ARRANGEMENTS 


the  boundaxy  condltlona  it^  a  fev  dlstlnet  poiats  ,  eoapletely  errcmeoua 
reeulte  are  obtained  unleae  a  eufficlent  mud>er  of  e^rdvlse  control 
points  are  used.  In  this  case  the  nuad>er  vas  insufficient,  so  that  the 
higher  texns,  uhlle  satisfying  the  boundary  conditions  at  the  control 
points,  Bade  natters  considerably  worse  ereiywhere  else. 

Consequently,  In  the  six  teft  rune  listed  in  Table  7.I,  the 
caaber  corrections  were  re^cosqputed  simply  by  deleting  all  of  the  terns 
in  the  circulation  distribution  except  (7<1)/  and  obtaining  the  eaad>er 
froB  the  average  value  of  9u^/ds  at  each  radius. 

The  casiber  factors  obtained  In  this  way  are  shown  in  Fig.  7«1* 

It  can  be  seen  that  the  results  obtained  fron  three  saallest  spacings 
(24  X  8,  24  X  6,  24  x  4)  all  agree  to  within  *  2$,  and  that  the  only 
large  error  occurs  with  the  coarsest  spacing  (8  x  4)  at  x  *  0.8^. 

While  the  characteristics  of  this  propeller  are  fairly  typical, 
this  one  set  of  testa  cannot  be  considered  as  establishing  the  convergence 
of  the  lattice  snthod  under  all  conditions.  Sowever,  fron  these  results 
It  In  tentatively  concluded  that  the  24  x  8  spacing  should  give  casiber 
corrections  which  are  within  *2$  of  the  values  tdilch  would  be  obtained 
fron  a  continuous  vortex  sheet. 

The  second  question,  nasnly,  whether  the  foxmalatloa  of  the 
llftlngisurface  problen  with  the  sli^llfying  assuaptlons  Introduced 
in  Chapter  1  is  an  adequate  representation  of  the  physical  problen.  Is 
soaethlng  which  Is  vexy  difficult  to  answer  due  to  the  laxge  nuad>er  of 
vaxrlables  Involved.  While  a  cosiparlson  between  theoxy  and  ej^eriaent 
nightie  succesiful  in  one  or  two  paztlcular  cases,  this  is  no  assurance 
that  agraenent  will  exist  in  general. 

Another  difficulty  results  fron  the  fact  tlxat  existing  experlnental 

•  * 

t 

data  Include  only  overall  neasiumients  of  thrust  and  torque,  so  that  it  is 
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Iji^sslble  to  determine  whether  the  desired  radial  iMd  distribution  has 
been  obtained.  The  first  saccessful  pressure  iseasureaents  on  a  rotating 
propeller  blade  were  sade  rfeen€ly  by  Aualaender^^^^  at  the  David  Taylor 
Model  Basin  with  fairly  elaborate  instruaentatlon;  however,  even  these 
results  contain  some  esqperlmental  scatter.  Evidently,  It  is  very  difficult 
to  locate  enough  pressure  taps  on  the  blade  to  detezrolne  the  lift  coefficient 
accurately.  The  transmission  of  pressure  readings  from  a  rotating  shSLft 
also  presents  a  dlfficxilt  Instrumentation  problem. 

In  the  present  work  caiii>er  corrections  ana  given  for  eight  propellers 
shoving  the  effect  of  a  few  of  the  many  possible  variables.  !Chese  propellers 
all  have  symetrlcal  B-3erles  blade  outlines  and  a  hub  radius  of  0.2.  The 
lattice  arrangement  Is  the  same  as  In  test  6  described  previously,  i.e.,  the 
finest  spacing  possible  with  the  current  program.  As  In  the  test  runs, 
the  higher  terms  In  the  clrc\ilatioa  distribution  were  deleted. 

Ihe  first  six  results,  shown  in  Elgs.  7*2  -  7*7,  are  for  optlmim, 
open  water  propellers  with  peurabolic  mean  lines.  These  results  Include 
a  limited  nuiiber  of  vcurlations  In  etpstnded  area  ratio  hydrodynamic 

advance  coefficient,  and  nusiber  of  blades.  Caid>er  correctlcms  given  by 
Van  Msnen^^^  and  Sckhardt  and  Morgan^^^  are  shown  on  the  same  plots  for  co^parli 

It  Is  evident  that  the  lattice  results  have  the  same  general  shape 
as  those  given  by  Van  Manen,  both  caiift>er  corrections  becoming  larger  near 
the  tip  of  the  blade.  The  Eckhardt  and  Morgan  results,  on  the  other  band, 
becoBK  more  or  less  constant  on  the  outer  regions  of  the  blade.  As  mentioned 
In  Chapter  1,  the  latter  corrections  are  derived  from  Ludwel^  and  Glnxel  results 

4l 

for  circulation  distributions  with  reduced  loading  at  the  tip.  Cosweqpiently, 
the  lattice  results  seem  to  substantiate  the  fact  that  a  large  caid>er  correction 
is  necessary  at  the  tip  In  order  to  achieve  an.  optimum  radial  load  distribution 
with  normal  blade  sh^es.  Bovever,  this  increase  near  the  tip  is  somewhat 
less  than  the  results  given  in  Reference  (3). 


0.a  0.4  0,6  *  0.8  •  1.0  0.4  0.6  0.8  1.0 

NON  DIMENSIONAL  RADIUS 


FiGS.  7.2 -7.5  CAMBER  CORRECTION  K  VS. 

NON  DIMENSIONAL  RADIUS  >? 


-66- 


NON  DIMENSIONAL  RADIUS 


FIGS.  76-7.9  CAMBER  CORRECTION  K  VS.. 

NON  DIMENSIONAL  .RADIUS  ^ 
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Fig.  7 <8  sliom  a  coaparlaoa  of  a  iiBke~adapt«Aa*>BA  opca-water 

propeller,  both  haring  the  eaae  adraaee  coefficient  at  x  "  0.7*  The  vake 

C  ^2) 

distribntian  ia  taken  from  the  nemerical  example  given  bjf  lecker'*^'.  The 
two  reaulta  azv  practically  identical.  Hbwerer,  the  vake  variation  in 
thia  exaagple  ia  fairly  amaU,  and  the  radial  load  dlatrlbution  la  almoat 
the  aaae  aa  In  the  open-water  eaae.  Ckunaequently,  It  ia  pbpaible  that 
■ore  extreme  vake  variatiooa  auch  aa  wouM  occur  vith  low  apeed  cargo 
ahlpa  al^t  affect  the  caid>er  correction. 

Finally,  the  effect  of  radial  load  dlatrlbution  ia  ahovn  ia 
Fig.  7 *9  for  two  open- water  propellera.  One  propeller  baa  a  reduced 
circulation  at  the  tip,  following  the  pitch  dlatribation  recoanwaded  by 
Bckhardt  and  Norgan^^  ^ .  The  other  la  an  optlmim  propeller  vith  the  aame 
adrance  coefficient  at  x  "  0*7 •  The  reaulta  ahov  that  a  reduction  in 
local  propeller  loading  tenda  to  reduce  the  cahber  correction,  and  vice  veraa. 

Ihe  reaulta  given  in  Figa.  7.8  and  7*9  were  obtained  with  a 
allghtiy  different  lattice  arrangement  conslatlng  of  aizteen  radial 
lattice  apaeea  with  additional  half  apacea  at*  the  enda .  Ihla  arrangement 
vaa  found  to  give  the  aame  reaulta  aa  with  twenty- four  equal  apacea,  but 
with  aomevhat  leaa  ccaiputation  time. 

To  teat  the  program  for  aaynmetrlcal  bladea,  two  propellera 
were  run,  one  with  a  aynmetrlcal  and  the  other  with  a  akewed  blade.  All 
other  characterlatica  ware  the  aame*  The  reaulta  ahowed  that  the  caaiber 
corx^ctlona  for  the  two  propellera  were  practically  identical.  loveWer, 
the  propeller  with  akewed  bladea  required  aa  additional  pitch  correction 
mf  about  2.5  degreea/unlt  lift  coefficient  near  the  tip.  Vhlle  thia 
correetlon  la  not  vef^  large.  It  indlcatea  that  a  pitch  correction  nl^t 
be  imcoxjocated  ia  tha  dealgn  of  propellera  with  a  large  amouat  of  akev. 
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ConclnBl&M 

On  th«  bMis  of  the  Halted  mniber  of  tnnierleal  results  described 
In  the  preceding  section,  It  atppears  that  the  rortex  Sinttlce  nethod  Is  a 
feasible  my  of  obtaining  lifting  surface  corrections  for  narlne  propellers. 
The  method  has  the  admatage  that  variations  in  blade  shape,  vake*  and 
clronlation  distribution  can  be  taken  Into  account.  The  nuaerlcal  exsaples 
given  Illustrate  the  fact  that  the  latter,  uhlch  Is  not  taken  into  account 
in  current  design  methods,  can  effect  the  lifting  surface  correction. 

It  is  therefore  recosnended  that  a  systematic  series  of  cal¬ 
culations  of  caaber  and  pitch  corrections  be  made  covering  a  vide  variation 
in  such  parasmters  as  nosiber  of  bledes,  pitch,  blade  shape,  and’ radial  load 
distribution.  These  results  smy  be  of  use  both  for  design  applications,  and 
to  determine  irtilch  parameters  cause  significant  differences  ,1a  the  lifting 
surface  correction. 

At  the  same  time,  these  results  vill  p>exBlt  cut  evaluation  of  the 
effectiveness  of  the  vortex  lattice  method  by  caBq)arl8oa  vlth  existing  ex¬ 
perimental  results . .  Hovever,  it  vould  also  be  desirable  to  build  and  test 
a  nuaber  of  model  propellers  designed  according  to  these  results.  These 
tests,  if  possible,  should  Include  pressure  distribution  meaisuraBeatm. 

However,  before  this  is  done.  It  Is  recommended  that  a  more 
accurate  treatment  of  the  hub  boundary  condition  be  Included  In,  the  lattice 
method.  As  Indicated  in  Chapter  the  lattice  method  developed  in*  the 
present  work  ttJces  the  hub  into  ascount  In  a  fairly  crude  way  slJi^ily  by 
requiring  that  the  aireulatloa  at  the  hub  be  sero  while  asglecting  the 
cooflltlom  that  the  radial  velocity  must  be  seno.  It  Is  believed  that  the 
presence  of  the  hab  cam  be  taken  into  account  by  a  discrete  source  dls- 
trlbtttlom  within  the  bub  cylinder*  The  strength  of  the  source  distribution 


•ad  the  ralae  of  the  olxeulatioa  at  the  hah  could  be  detexvlaed  hjr  iaeludlng 
control  points  oa  the  hub  cyllader  la  addition  to  those  on  the  blade  surface. 
This  added  reflneaent  should  not  greatly  Increase  the  cosvlexlty  of  the 
coiqputatlons  t  and  should  produce  aore  accurate  results  In  thS  Inner:  part 
of  the  %lade . 

It  Is  also  reeoaaeaded  that  the  lifting  surface  prograas  be 
aodlfled  to  accoaaodate  finer  lattice  spaclags  vith  an  Increase  la  the 
nuniber  of  chordvise  control  points  la  order  to  obtain  additional  teras 
In  the  Fourier  Series  for  the  circxilatlon  distribution.  This  would  also 
provide  as,  additional  cheek  on  the  accuracy  of  the  caailber  corrections 
obtained  with  the  preseat  prograas. 


A 


jp 
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j^CTDIX  A 
PROQKAM  DBCBIPriOia 


Introduction 

Digital  congniter  programs  vere  prepared  $o  obtain  numerlced  solutions 
of  the  following  ^hree  problems 

a)  Determine  the  non-d^menslonal  radial  circulation  distribution 
for  a  lifting-line  propeller  with  a  prescribed  distribution  Of 
tan  3  and  tan  3^. 

b)  Determine  the  camber  and  pitch  correction  for  a  propeller  with 
an  arbitrary  blade  outliae,  tan  3  and  tan  3^  and  mean-line  type. 

c)  Determine  the  camber  correction  under  the  same  conditions  eis 

* 

(b)^  but  for  the  special  case  of  a  symmetrical  blade  and  a 
mean-line  ^dilch  is  synaetrlcal  about  the  mid-chord. 

A  number  of  other  programs  were  prepared  to  test  Tarlous  features 

of  the  vortex  lattice  method,  however,  these  are  not  of  sufficient  general 

*  • 

Interest  to  be  reported. 

The  above  programs  were  prepared  for  Wsa  with  the  IBM  709  Data 
Processing  System  at  the  M.  I.  T.  Cooptation  Center,  and  were  run  using 
the  Fortran  Monitor  System,  nie  principal  source  program  language  was 
FORTBAH,  however,  some  of  the  programs  vere  written  in  FAP  in  order  to 
perform  certain  operations  not  within  the  scope  Of  the  FOBTBAN  language. 
Descriptions  of  these  systems  appear  In  References  (37 ),  (33),  (39)>  auxd 
(1^). 

Programs  (b)  and  (c)  werm  also  modified  for  use  with  the  IBM  7090 
Installed  at  the  David  Taylor  Model  Basin,  and  seme  of  the  results 
shown  in  Chapter  7  were  Obtained  there. 

Bach  of  the  three  p^ograas  oonslsts  of  a  number  of  specially  prepared 
subroutines  as  well  as  standard  lArary  routines.  In  some  cases  the  same 

a 

subroutine  cam  be  used  in  all  three  programs. 


Brief  descriptions  of  the  prlnelpi^  subrcutlnes  vUl.  be  given 
In  the  foUovlng  sections.  Hbvev’er>  these  sections  are  Intended  only 
^o  indicate  the  genscnl  node  of  ^eratlon  and  references  to  eongputer 
lAngoage  irill  be  avoided*  Idhstlngs  of  the  source  prograas  are  given 
In  Appendix  B 

Helical  Vortex  Integration  * 

fOie  hellccLl  vortices  are  divided  Into  two  parts;  'thi  j^art  on 

m 

the  blade  which  extends  between  the  boimd  vortex  elestents  closest  to 
the  leading  and  trailing  edges,  and  a  dawnstream  igtart  idilch  stix^ii  at 
the  bound  vortex  nearest  the  trailing  edge  and  extends  an  Infinite 
distance  dovnstreaii.  As  Indicated  In  Chapter  2,  the  velocity  Induced 
^  the  helical  vortices  on  the  blade  Is  obtained  entirely  by  nKUserleal 
Integration,  idiUe  the  Integration  of  the  downstrean  hellcei'  ls-^j^r> 
forsKd  by  nuaerlcal  Integration  up  to  a  sufflclsntly  large  '^tilii^'^df  9, 
and  the  ieaialnlng  contribution  estlnated.  *  » 

m 

It  Is  assdBted  that  the  nusKrlcal  Integration  can  be  truncated 

* 

within  the  first  six  revolutions  downstrean,  l.e.,  9^  <  12  tt.  Conssquently, 

It  will  be  sufficient  to  divide  the  Interval  from  the  bound  vortex 
nearest  the  leading  edge  to  a  point  six  revolutions  downstrean  into  a 
sequence  of  ^-potnt  Gauss  ordinates.  At  each  ordinate,  the  functions 
F^(9j^)  and  the  weights  defined  In  (2.19)  and  (2.21)  are  to  be 

cooputed.  Bach  Integration  nay  then  be  perfozned  by  coa^uting  the  con- 

¥ 

stants  c^  and  d^  defined  In  (2.l8)  and  applying  (2.21). 

For  the  downstream  Integration  It  has  been  found  eiqplrlcally 
that  the  angular  Intervals  shown  In  Table  A-1  when  subdivided  Into  5»polnt 
Gauss  ordinates  result  In  total  accusailated  Integration  errors  of  less 
than  *0009  In  the  non-dlmenslonal  Induced  velocities  defined  In  (2.10)  -  (2.12) 


ftible  A^l  gpaclag  For  Bbaerleal  Integimtlon 

(In  Degrees) 


let  Rerolution  -  Coarse  ^acjlng  -  .2$  <.  1  -  T|| 

0,  <>0,  50,1  90,  180,  270,  380 

1st  Revolution  -  Medium  gpaciig  -  .10  <  |l  -  T||  <  .25 

0,  5,  10,  20,  kO,  60,  100,  150,  200,  270,  360 

1st  Rerolutlon  -  line  Spacing  -  .02  <  |l  -  T||  <  .10 

0,  1,  2,  4,  7,  10,  20,  30,  50,  75  >  100,  150,  200,  250,  300,  360 

2nd  -  6th  Bevolutlon  -  .02  <  |l  -  T)! 

Table  A-2 

Welf^ts  and  Ordinates  for  Legmdre-Qauss  Integration  foriles 


K 

Weight,  W^ 

Ordinate, 

1 

.Il8it6h 

.0^<6910 

2 

.23931^ 

.230765 

*■  5  point  rule 

3 

.aBWfh 

.500000 

4 

k 

.23931^ 

.769235 

^  ' 

5 

.llBii6i» 

.953090 

1 

2 

•5 

.5 

.268675 

.711325 

1 

1 

^  2  point  rule 

-96- 


Ibe  velghts  and  ordinates  for  an  intezral  of  unit  length  is  given  in 
Table  A-2.  *  ^ 

From  these  tvo  tables >  a  set  of  values  of  9^  nay  be  obtained. 

For  each  9^  there  will  be  seven  elements  of  and  one  weight  so 
that  there  will  be  a  total  df  forty  numbers  associated  with  each  five-point 
Qauss  interval.  Hie  toted  downstream  integratioa  table  consists  of  l>dhO 
elements . 

The  portion  of  the  helical  vortex  on  the  blade  is  subdivided 

<e 

into  a  nuinber  of  elements  lying  between  bound  vortex  elements.  These, 
together  with  the  six  downstream  revolutions,  are  shown  schematically 
^  Fig.  ArX.  The  maximum  number  of  chordwise  bound  vortex  elements  is 
assumed  to  be  eight,  so  that  a  totad  of  fifteen  intervals  on  the  blade 
is  possible . 

The  emgular  Interveds  on  the  blade  depend  on  the  geonetzy  of 
the  blade  and  will  in  general*  be  different  at  each  radius .  Consequently, 
it  is  inqposslble  to  subdivide  these  Intervals  into  a  fixed  number  of 
Gauss  ordinates.  In  this  case,  the  minimum  number  of  Gauss  Intervals 
is  determined  such  that  the  spcusing  will  not  exceed  the  initial  speu:lng 
necessary  for  the  downstream  integration  for  each  of  the  three  ranges 
of  |l  -  T||  Since  it  is  possible  that  many  of  the  interveds  on  the 
blade  will  be  very  small  (such  as  11  and  I3  in  Fig.  A-l),  provision 

is  made  for  using  a  2-polnt.  Gauss  Buie  if  -^e  Interval  is  less  than  koft 

♦ 

of  one  ^-p6int  Gauss  Interved.  Finally,  if  the  Interval  is  leaS  them 
.556  of  a  5*polnt  interval,  the  Integral  is  approximated  by  its  mean  vadue. 

It  Is  obvious  from  gecxnetfieed  considerations  that  the  parameter 
|l  -  1)1  used  in  selecting  the  integration  spacing  is  applicable  only  to 
the  index  blade*  It  hem  been  found  that  the  integration  of  the  helices 
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on  the  other  blades  may  he  done  with  the  coarse  spacing  tgx  ill  raluii 
of  |l  -  7||  vlthodt  altering  Uie  final  result. 

She  integration  of  the  helical  vortices  requires  three  suhroutlnes. 
Ohe  downstream  integration  table  is  generated  by  a  subroutine  called 
and  this  needs  to  be  called  only  once  at  the  beginning  of  each  run.  The 
Instruction  CALL  HUlffiUO  (P^  L)  causes  the  1840  eleaents  of  the  table  to 
be  cosguted  and  stored  in  Increasing  neisory  locations  starting  at  P. 

Location  L  Is  the  first  element  of  an  "address  directory”  which  requires 
sixty-three  storage  locations  in  decreasing  munerical  order  starting 
at  L.  !nie  "Address  directory"  is  a  (21  x  3)  array  corresponding  to  the 
twenty-one  possible  integration  Intervals  shown  in  Fig.  A-1  and  the  three 
possible  spaclngs.  Each  element  of  the  array  ccatalns  the  starting 
address  of  the  integration  table  for  that  Inteival  as  well  as  the  nuober 

S' 

of  eulgles  9^  in  the  interval..  Subroutine  HUMBUG  fills  in  only  the 
first  (6x3)  elements^  which  correspond  to  the  downstream  part  of  the 
integration. 

Subroutine  LIST  does  more  or  less  the  same  ^iBg  for  the 
Intervals  on  the  blade.  Ihe  calling  sequence  is  * 

CALL  LIST  (irSPACE^  ANGLES^  L) 

where  I3PACE  is  the  number  «f  sx>aces  on  the  blade  (which  cannot  exceed 
15);  AEGUS  is  the  first  eleomnt  of  a  list  of  angles  defining  the  limits 
of  each  Interval,  and  L  is  the  "address  directory"  idilch  Is  -Uie  same  as 
In  the  calling  sequence  for  HUMBUG.  The  list  of  angles  starts  at  the 
bound  vortex  element  nearest  the  leading  edge,  and  Is  stored  In  decreasing 
mesioiy  locations.  These  are  all  angles  in  radians  relative  to  the  emgle 
of  the  trailing  edge  element,  and  will  consequently  all  be  <  0.  Subroutine 
LIST  determines  the  numiber  Of  Integration  spaces  4ji  each  Interval,  ccxnputes 
the  functions  and  and  stores  them  immediately  fbllowlng  thf.  functions 


'4^ 


generated  ty  BDHBUO  and  conpletea  the  stxtytthree  element  address  dlreetozy. 
If  the  jn^tloD  table  being  generated  beglnf  to  exceed  the  slse  of  core 
storage,  an  #x^i;or..stop  resu^s*  Ibis  subroutine  Is  called  once  at  each 
lattice  radius. 

Ihe  actual  Integration  is  perfonsed  in  subroutine  HILIX  vhleh 
is  called  as  foUovs: 


CALL  HELIZ  (XIA,  lASBlO,  lAXBZ,  COSBI,  FHIZ,  IQ,  MSIACI,  L^  UV) 

idiere  the  following  arguments  are  as  defined  in  Chapter  2: 

% 

Kd  -  H  lAlBIO  -  tan  TAIBZ  -  tan  0^ 

COaBI  -  cos  0.  IHIZ  -  0)  IG  -  g  01  -  u. 

i  O  u 

The  arguments  ISFACl  and  L  are  the  same  as  in  LIST.  Ihe  angle  is 
measured  from  the  particular  control  point  to  the  start  of  the  downstream 
helix,  in  accordance  with  the  notation  of  Fig.  2.1.  Ihe  syid>ol  01 
denotes  the  first  of  a  sixte^-element  array  stored  in  decreasing  memory 
locations. 

HILIX  starts  by  coasting  the  constants  c  aM  d  .  The  integration 

M  21 

is  then  performed  according  to  (2.21)  using  the  "address  directory"  to 
locate  the  pre'-co^puted  functions  and  to  determine  the  nuiiber  of  points 
In  each  In^rral.  The  first  downstream  interval,  designated  by  1  in 
Fig.  A-l  is  coiqnited  first.  If  |l  >  1||  >  .25  all  g  blades  are  Integrated 

simultaneously.  If  |l  •  T||  <  .25  all  but  the  index  blade  are  Integrated 
using  the  coeurse  spacing,  and  the  index  blade  is  then  cosguted  using 
the  medium  or  fine  spacing.  After  each  downstreem  revolntlon  has  been 


coi^eted^  the  integral  to  infinity  is  estimated  from  the  relation 
g  cos  (tan  Ian  0  -  1|) 


u 


ta^  0 


(A.l) 


io 


vhlch  is  obtained  from  (2. 35),  (2>36)  and  (2.23).  ^Aien  tiro  8uc<ie88iye 
e8tlnate8  agree  to  within  .000^,  the  downatream  integral  ie  aasumed  to 
have  converged.  If  the  number  of  spacea  on  the  blade  la  zero,  aa  would 
be  the  ease  in  lifting  line  theory,  the  integration  la  complete.  Otherwlae 
the  interved  cloaeat  to  the  trailing  edge,  dealgnated  by  7  A*1 

la  Integrated  uaing  the  funetiona  coasuted  by  LIST.  Thia  proceaa  ia 
repeated  tor  all  the  remaining  Intenrala  up  to  the  bound  rortex  neareat 
to  the  leading  edge.  Ihe  reault  of  the  preceding  Interwal  la  added  to 
each  new  Interval,  ao  that  the  z«Bult  la  a  table  Of  the  Integral  from 
(AIQIJI)^  to  ».  Thia  la  atored  In  decreaaing  memory  locatlona  ataztlng 
at  UV.  The  flrat  element  of  U1  containa  the  value  of  the  Integral  frcm 
the  ^eauUng  edge  bound  vortex  to  infinity. 

Hit  time  required  to  pez^orm  the  helical  Integration  dependa 
on  the  pitch  angle  and  the  number  of  blsdea.  fOr  a  three-bladed  pro¬ 
peller,  the  downatream  integration  talcea  rcughly  1  -  I.5  aeconda  on  an 
IBM  709.  The  Intcgratioa  on  the  blade  ia  mich  faater,  and  a  typical 
average  time  including  both  downatream  and  on-blade  Interveila  la  0.2^ 

It- 

aeconda  per  Interval  for  a  three-bladed  propeller.  Thia  includea  a  pro¬ 
rated  amount  of  the  time  apent  In  the  data-generatlng  ‘aubroutlnea  ffiJMBUG 
cmd  LIST.  A  alx-bladed  propeller  would  take  a  little  leas  than  twice 
aa  long.  Listings  of  HUMBUG,  LIST  AID  HBUX  appear  in  Appendix  B. 

General  Lifting  Line  Program 

This  program  foxs>s  and  solves  the  set  of  e(iuatlons  given  in 

(i.20),  using  the  helical  Integration  aubreutinea  previously  deacflbeft. 

•  • 

The  input  data  conalate  of  a  list  of  nine  values  of  the  non-diaenalonal 
radius  x>  vlth  corresponding  values  0^  tan  0^  and  tan  0.  The  remaining 


data  coniiits  ot  the  nuidser  of  g,  the  xmaber  of  lattice  qpaces 

It,  the  noaibftr  of  control  points  P  and  a  list  of  the  yalues  o^  M  containing 
control  points.  If  the  pitch  of  the  free  yo:i^ex  ^stea  Is  constant,  the 
first  elesKnt  In  the  list  of  tan  3^  ssiy  he  replaced  with  the  adyanee 

,  t 

coefficient  and  the  remlnlng  elements  of  tan  3^  and  tan  3  l<Kf^  blank. 
The  result  in  either  case  Is  a  table  of  the  non>dlaensional  circulation 
G  defined  in  (4.2)  as  veil  as  the  Fourier  coefficients  qjT  G«  In  addition, 
If  tan  3/^0,  the  circulation  Is  also  expressed  the  fora 


arrinr. 


tan  3^ 
tan  3 


-  1 


(A.2) 


in  accordance  vlth  the  definitions  In  (l)  and  (5)>  If  X^  is  glyen, 
the  propeller  is  assuaeA  to  be  optlmia  and  the  Goldstein  factors  k  are 
ccavuted  froa  (4.l4). 

Since  the  input  data  Is  not  necessarily  at  '^e  saae,  set  of 
radii  as  required  for  the  lattice,  the  required  values  of  tah  3^  *nd. 
tan  3  aire  obtained  by  three-point  Lagraogian  Intezpolatlon.  In  addition, 
_8lnce  the  conyerslon  froa  the  actual  radius  r  to.  the  transfozaed  zadlus 
p  according  to  (4.3)  occurs  very  frequently  In  both  the  lifting  line 
and  lifting  surface  programs,  the  tsansfozaatloa  Is  performed  in  a 
subroutine  called  HAP.  Finally,  the  printed  output  froa  this  program 
Is  controlled  by  a  subroutine  celled  HAITXR. 

The  cosputatlon  time  in  minutes  on  an  IBH  709  can  be  approxlaated 
by  the  following  relatlcm 


T  -  ^  (.7  +  .2Ai) 


(A.3) 


A  listing  of  the  progPaas  apd  a  saaqle  set  gf  resists  appear  in  Appendix  B 


Lifting  Surface  Programs 

Two  lifting  surface  prograas  vere  prepared,  one  corresponding 

to  the  general  caise  covered  In  Chapter  and  the  other  fof*  the  special 

0 

case  of  a  syBButrlcal  hlade’  as  discussed  In  Chapter  6.  Since  both  ^ro~ 
grass  are  pm^tlcally  the  saiie,  the  general  discussion  la  tills  section 
will  apply  to  both  unless  specifically  indicated  otherwise. 

The  Input  Includes  a  list  of  nine  vedues  of  x  together  with 
corresponding  values  of  teui  aad  tan  3  as  In  the  lifting^  line  case. 

Ill  addition,  the  chord  lengths  I/D  at  each  value  of  x  required  as 
well  as  the  chord  load  factors  defined  In  (^.19).  In  the  general 
p;rograai,  the  mid-chord  auigles  6  shown  In  Fig.  ^.1,  and  the  radial  load 
distribution  must  be  given  at  each  value  of  x*  latter  may  be  given 

in  the  form  of  Goldstein  factors  k,  or  either  non-dimensional  circulations 
(f  or  G' . 

In  the  ‘synmetrlcal  blade  program,  the  mid-chord  angles  are 
zero  by  definition  and  need  not  be  given.  The  other  difference  is 
that  the  Fouriar  coefficients  of  G  are  given,  rather  than  G  itself.  This 
avoids  the  Inaccuracies  introduced  by  interpolation,  since  the  total 
strength  of  the  bound  vortex  elements  at  a  particular  radius  will  be 
exactly  the  tame  as  In  tha.  lifting  line  case  with  the  same  radlaX  lattlffe 
arrangement .Finally,  the  slopes  of  the  mean  llna  with  unit  caoiber  ratio 
h^  defined  In  (5-22),  the  camber  ratio  for  unit  lift  coefficient  In 
tvo-dlmeilnloiBl.  flow  and  the  constants  defining  the  lattice  and  control 
l>olnt  arrangement  must  be  given. 

♦ 

In  either  case  a  main  program  reads  the  data  and  cos^tes  the 
various  geometrical  properties  associated  with  the  lattice  8u:raDgement . 
SilUk  angles  and  chord  length#  at  each  of  '^e  lattlbe  radii  are  obtained 


by  parabbllc  Intei^latlon.  In  the  general,  program  a*  sijbroutlne  chUed 
AKL  ccagnites  and  solves  the  set  of  equations  glfen  In  (^>32).  In  the 
synmetrlcal  blade  case  a  similar  subroutine  called  CAMQW  congnites  and 
solves  tlii  equations  given  In  ^^.liB). 

The  only  elements  in  (^.32)  and  (6.16)  which  require  any 

flgnlflcant  amount  of  cooqputatlon  are  tlte  velocities  Induced  the 

horseshoe  elements,  u _  .  As  can  be  seen  from  Jigs.  S,1  or  5*2,  these 

’  mnpq  -w  ^  w  • 

consist  of  two  semi* Infinite  helical  vortax  segmentsooonnected  by  a 
radial  bound  vortex.  The  velocity  contribution  of  the  bound  vortex  inay 
be  obtained  explicitly  by  evaluating  equations  (3 >9)  end  (3 <10),  and  this 
may  be  dene  very  easily  In  a  subroutine  ceilled  BOUND.  The  velocity  induced 
by  the  helical  segments  may  be  obtained  from  the  subroutine  HELIX  described 
.previously.  However,  connecting  the  right  helical  segment  to  the  right 
horseshoe  requires  a  little  bit  of  bookkeeping  since  the  order  In  which 
the  radial  vortices  Intersect  a  pe^lculeu:  helical  vortex  from  above  and 
below  depencLs  on  the  outline  of  the  blaide. 

She  cooptation  time  required  in  minutes  on  an  IBM  709  can  be 
estimated  by  the  following  relation* 

T  »  .62  +  .0033  (P<3M  (9  +  *)  )  (A. 4) 

idiere  the  symbols  are  as  defined  in  Chapter  5*  This  equatigo  holds 
for  both  the  general  and  wynBetrlcal  blade  programs  pz^ided  N  is 
Interpreted  as  the  total  number  of  chordvlse  vortices.  A  listing  of 
the  programs  .fox  coopting  the  syometrlcal  blade  case,  and  a  sample  set 
of  results  appears  in  Appendix  B.  The  programs  for  the  gaaaial  case  are 
very  similar,  and  will  therefore  not  be  included. 
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TABLE  B.l  LIFTING  LINE  MAIN  PROGRAM 


OIMENSICN  FfLLdiBOOO)  »X(9)  •XTriI'P)  iXTR(9)  .DUMMY  (  9  *3  I  (  25  ) 

1  *TANRZ(?5) .R(24) .RH0(?51 »TANBI (24) .TrFTA(24) iCOSRI  ( ?4 ) *R( 8) 

2  »ZETA(  8  .;>?))  •U(e*?5)iA(B.R)»P(16)  .GAMMAIR)  tGDTMB(9)  *  ANSI  9.5)  » 

3  MC(8).LZ(70) 

COMMON  FILLtPZ.LZ.ANS  . RZ » TANRZ »R . RHO . T ANB I . TBETA .COSBI .B.ZFTA » 

1  U» AiEtG»PUlZ»ALAM*RH»ZT.TFMP.AMT.DFLM»HDELM.Y.AI .SNl .SN2.TDEL 

2  »Tri .CBI »TbZ.FTa»WN.DFI .MC *NST0P .MT »NPT » NG .NTMl .MOPT 

FOUi valence  ( X. dummy »ANS ) . (XTRI .ANSI  10) ) . I XTr . ANSI  1 9 ) ) . I  GAMMA. ANSI 
128) ) .(GDTMfi.ANS(37) ) 

CALL  OCTALS 
CALt  STOMAP 
CALL  HUMBUGIPZ.LZ) 

1  CALL  CL0CKI2) 

READ  INPUT  TAPE  4.101.NSTOP 
■I'di  FORMAT!  ID  * 

IFINSTOP)  14.24.14 

24  READ  INPUT  tAPF  4 . 1 00 . 1  X  I N ) . N  =  1 .9 ) . I XTR I  I N ) . N= 1 .9 ) . I XTR I N ) .N  =  l .9 )  . 

1  MT.NPT.NG.IMCIN) .N=1.8) 

100  F0RMAT(3(9F8.6/')D  14) 

MAX  =  MT  +  1'* 

G=NG 
NTMl  *0 
PHIZsO.O 

ALAM=X(6)*XTBI (6) 

RH=X(1). 

MOPT=0 

1F(XTBI(2))  7.2.7 

2  ALAMsXTBKl) 

MOPTsl 

IFIRH)  3.4.3 

4  X(l)=.01 

3  DO  5  N=1.9 

XTBI  IN  )  =  AtAM/X(li()  * 

5  CONTINUE 

DO  36  M=1.MAX 

DO  36  1  =  1. NPT 

ZETAI I .M)=1.0 
36  CONTINUE 

7  DO  15  N=l»9 

XTBI IN)=XTBI (N)»X(N) 

XTP(N)=XTR(N)#X(N) 

15  CONTINUE 

DO  16  M=1.3 

DO  16  N=1.4 

K=10-N 

TEMPxDUMMYIN.M) 

DUMMY  I N . M ) =DUMMY I K . M ) 

DUMMY(K.M)=TEMP 

16  CONTINUP 
AMT=MT 

DELM=I l.-RH) /AMT 

HDFLM=.5*DELM 

AM=RH-HDFLM 
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RZ( 1 )=RH+.25»HDELM 
TEMP=RZ(1) 

CALL  1NTERP(TEMP.Y.X.XTRI *3*9) 

TANBZ(1)=Y/RZ( 1) 

DO  9  M=1*MT 

R(M)=AM+DELM, 

AM=R(M) 

'  TEMPsAMe 

■  CALL  MAP(TFMP*RH)  „ 

RHO(M)=TEMP 
RZ(M+1 )=R(M)+HDELM 
IF(M-MT)  19*10*19 
10  RZ(M+n=RZ(M+l)-.25»HDELM 
19  TEMP=RZ(M+1) 

CALL  INTERP(TEMP*Y*X*XTBI *3*9) 

TANBZ(M+1)=Y/RZ(M+1) 

TEMP=R (M ) 

CALL  INTFRP) TEMP*Y,X*XTR1 *3*9) 

TANRf (M) =Y/R(M) 

CALL  INTFRI»tTFMP*Y,X*XTP*»*9» 

TRETA(M)=Y/R(M) 

COSBI (M)=1./S0RTF(1.+TaNRI (M)**2) 

9  CONTINUE 

DO  6  I»1*NPT 

MS=MC( I ) 

TDEL=R (MS)#( TANBI (MS ) -tret a (MS ) ) 

TBI=TANBI (MS) 

C3I=C0SBI (MS) 

R( I )=2.*R(MS)*CBI 
DO  6  M=1*MAX 

IF(MOPT)  18*18*17 

10  ZFTA( I *M)=R(M)«(TANRI (M)-TRFTA(M)  )/TnFL  • 

17  TBZ=TANBZ(M)  .  . 

ETA=RZ(M)/R(MS) 

CALL  HELIX(  FTA*TbZ*TRI  *CRI  *PH  1  z', N(5  . NTMl  *LZ  *  WN  ) 

U( 1 *M) =WN 
6  CONTINUE 

DO  8  I=1*NPT 

DO  8  •  ♦  K=1*NPT 

A( I .K) =0.0 
8  CONTINUE 

RH0(MAX)=0.0 

DO  34  I=1*NPT 

AI-I 

SN1=0. 

D()  34  M  =  1*MAX 

IF(M-l)  31*31*32 

31  J='m 

GO  TO  33 
3'2  J=M-] 

33  SN2=SINF(AI*RHO(M) ) 

DO  30  K=1*NPT 

A()C  *  I  ),  =  A(K*  1  )+U(K*M)*(SN2*ZETA(K*M)-SN1»ZETA(K*J)  ) 

30  CONTIflUE 
SN1=(SN2 

34  CONTlAjUE 

WRlTp  C(UT.PUT  TAPF  2  *  1  0?  *  (  (  A  (  K  ♦  I  )  *  I  =1  *NPT  )  .K  =  l  *NPT  ) 
102  format (4G15. 8) 
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DET=1.0 

ME=XSIMEOF(e.NPT.l .A*B.DET.E) 

GO  T0( 12.11 .11 ) tME 

11  CALL  ERR0R(20H  ERROR  IN  XSIMFQF) 

14  CALL  EXIT 

12  DO  13  M=1  .9 

GAMMA (M> *0.0 
TEMP*X!M) 

XTBI (M)=XTBI (M)/TEMP 
XTB(M)=XTB(M)/TEMP 
CALL  MAP(TEMP.RH) 

00  20  1  =  1. NPT 

AI  =  I 

GAMMA! M)=GAMMA(M)+S INF (AI*TEMP)»A(  1,1) 

20  CONTINUE 
IF(MOPT)  14.22.21 

22  GDTMR (M)=( (XTBI(M)/XTB(M) ) , 0 ) *GAMMA ( M ) 

GO  TO  13 

21  GDTMb(M)=( (X(M)»*2+ALAM*»2)/(2.*X{M)**2*ALAM) )*GAMMA(M)*G 

13  CONTINUE 
CALL  WAITER 
GO  TO  1 
END 
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Table  b.2  main  program-  lifting  surface-  symmetrical  blade 


DIMENSION  FILL! 8000 ) . A ( 5g  *  86 ) ♦ AMS ( ?4 ) > R  (  86 ) »RUG( 8) *  CHORD ( 24) » 

1  COSR I (;4 ) •COFFZ ( 8 ♦?) »D( 74 ) » F ( 86 ) » F ( 86)  »HMU( a *7  )  tHI 7  )  * 

?  PSr(]6)»P(l6)»PSrR(8)«RZt28)*R(?4) *RHO( 74 ) .SNRHOI R)*TANBZ(25)* 

8  YSffM  (  74)  *TrETA.(74)  >T1L(  74)  *THFTA  t  8.7)  .  U  (  24 . 8 . 7  )  » WN  (  I  6  )  »  X  (  9.)  » 

4  XC0RD(9) .XTni (9) .XTR(9) .XGAMla) »XRH0(9) *MC(8) »NFLIP(16) » 

5  LZ(68) »DUMMY(9.8) 

COMMON  FlLL.PZ.LZ.A.ANS.B.nUG.COEFZ .COSBI . CHORD .D . E . F .H .HMU 

2  .P.PSI .psib.r.rho.rz.snrho.tanbi .TANBZ.TBAR.THETA.TIL.U.WN.XGAM 

3  .AL AM. ANT. AMT, AW. A  I .ANGLE.CBI .DELM.DET »ETA .G.GNZL.HOELM.PHI? »  ,> 

4  RH.RBl .RB2 .TBI .TBZ.TFMP.UB.W.Y.ZETA.MC.NFLIP, JIN.JOUT. JT .XTEST  .' 

5  KlOl .MS.NBOTH.MT.NT.MPT.MZi .N22 .NG.NOT.NTT ,NTM1 .N I P . NF , TBETA 
EOU  I  VALENCF.IX,  FILL.  DUMMY)  .  (  XCORD ,  E  .  DUMMY  (  1 0  )  )  .  (  XTB I  .  DUMMY  (  19  )  )  . 

1  (XTb,DUMMY(28) ) ,fXRH0,DUMMY(37) )  ® 

CALL  OCTALS 
CALL  STOMAP 
CALL  HUMBUGIPZ.LZ) 

CALL  CL0CX(2) 

JIN=4 

J0UT=2 

?0  READ  INPUT  TAPE  J I N , 1 00 , ( X ( N ) , N= 1 . 9 ) . ( XCORD ( N ) . N= 1 , 9 )  .  ( XTB  I  ( N  )  ♦ 

1  N  =  1 ,9) . ( XTB( N)  ,N=1  .9  ) 

READ  INPUT  TAPE  J I N , 101 , K T FST , MT , NT , NPT ,NZl . NZ2 . NG . ( MC ( N ) . N= 1 . 8 ) , 

1  GNZL 

NOT=(NT+NZl-NZ2-2 ) /NZl 

G=NG 

JT=NOT 

NB0TH=NT+NT 

NTT=NBOTH  +*  NBOTH 

ZETA=0, 

ALAM=0. 

NTM1=NTT-1 

RH=X(1) 

READ  INPUT  TAPE  J I N , 102 . ( ( HMU ( N . J ) , N=1 . NT ) . J=1 . JT ) . t H ( N ) ,N= 1 , NOT ) . 
1  (XGAM(N) .N=1.8) 

DO  51  N= 1.8 

C0EFZ(N.1)=XGAM(N) 

DO  51  J=2.7 

COFFZ(N,J)=0. 

51  CONTINUE 

IFtXTBKZ))  7.2,7 
2  ALAM=XTBI ( 1 ) 

IF(RH)  3.4.3 

4  X(1  )=-,.0l  • 

3  DO  5  N=1  ,9 

XTfll (N )*ALAM/X(N) 

5  CONTINUE 
ZETA=-1.0 

7  DO  15  N=l»9 
TEMP=X(N) 

CALL  MAPITEMP.RH) 

XRHO’iN  )  =  TEMP 

XTBI (N  )  =  XTBI (N)»X(N) 

XTB(N) =XTB(N)*X(N) 

15  CONTINUE 
DO  16 


O 


M=1  #5 
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DO  16  N=l»4 
i<:  =  io-N 

tfmp=dummy;n»M) 

DUMMY (N*M)=DUMMY(<.M) 

DUMMY( K»M) =TEMP 
16  CONTINUE 
ANT=NT 
AMT=MT 

DELM=( l.-RH) /AMT 

HDELM=.5*DELM  „ 

AM=RH-HDELM 

RZ ( 1 ) =RH+.25*HDELM 

TEMP  =  RZ( 1 ) 

CALL  INTERP( TEMP»Y,X.XTRI *3,9) 

TANBZ( 1) =Y/RZ( 1) 

DO  9  M=1*MT 

R(M)=AM+DELM 

AM=R(M) 

TEMP=AM 

CALL  MAP(TEMP.RH) 

RHO(M) =TEMP 
RZ (M+1 )=R(M)+HDELM 
IF(M-MT)  19*10.19 
10  RZ(M4l  )  =  RZ(M4-1)-,25*HDELM 
19  TEMP=RZ(M+1) 

CALL  INTERP(TFMP.Y,X.XTB1' .3*9) 
tANBZ(M+l )=Y/RZ(M+1  ) 

TE.MP=R  (M) 

CALL  INTERP(TEMP,Y.X.XTBI  .3..9) 

TANBI (M)=Y/R(M) 

CALL  INTFRP(TEMP.Y,X,XTB.3 *9) 

T.8ETA(  M)  =Y 
TEMP=RH0(M) 

CALL  INTERP( TEMP.Y.XRH0.XC0RD.3*9 ) 

CHORD(M)=Y 

COSBI (M)=1./S0RTF(1.+TANRI (M)**2) 

9  CONTINUE 

DO  6  N:*1.NPT 

M=MC(N) 

D( N)=0. 

DO  6  I=1»NPT  ‘ 

AI  =  I 

D(  N)=D(N)+XGAM( I )»SINF(AI*RHO(M) ) 

6  CONTINUE 

DO  30  .  N=1.NT 

K=NBOTH-N+l 
DO  30  ‘  J=1.JT 

HMU(K» J)=HMU(N.J» 

30  CONTINUE 

write  output  TAPE  JOUT ♦ 1 03 .NT .MT . NPT .NZ 1  * NZ2 . i MC ( N ) . N= 1 ♦ 8 ) .NG ♦ 
1  ALAM*RH*GN2L 

write  output  tape  'jOUT. 104  .  (CH0RD(  N)  *N=1  .MT  ) 

WRITf  output  Tape  JOUT.105*(TANBI(N) .N=1 *MT) 
write  output  TAPF  JOUT.106*(TBFTA(N»  .N=1*MT) 

WRiTf  output  tape  JOUT,107*{D(N»,N=1  .NPT) 

write  output  tape  JOUT.108»( (HMU{N,J),N=1.NT) .J=1,JT) 

write  output  tape  JOUT.109.(H<N) .N=1 .NOT) 
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CALL  CAMBER 
CALL  CLOCK) 1» 

GO  TO  20 

100  FORMAT (9F8. 8) 

101  FORMAT (1514*F8. 6) 

102  FORMAT(7F10,7) 

103  FORMaT!6HO  NTsn.SH  MT=I2.6H  NPT=I1*6H  NZl  =  Il»6H  NZ2»n*5H 

lC=8I4t5H  N(*=I1.7H  ALAM  =  F6.4*5H  RH=F!f.3»7H  GNZL=F6.4) 

104  FORMAT(8HO  CH0RD=10F10.6) 

105  F0RMAT(8H0  TANBI=10F10.6) 

106  FORMAT(8HO  TBETA=10F10.6 ) 

107  FORMAT(8HO  GAMMA=10F10.6) 

108  FORMA^CaHO  HMU  =l€)F10.6) 

109  FORMAT(8HO  H  =10F10.6» 

END 


Table  b.3 
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WAITFR  SURROUTtNF 


SUBROUTINE  WAITER 

DIMENSION  FILLlflOOO)  tX(9»  *XTP,I  (o)  ,XTP(9)  *  DUMMY<9  »3  )  »RZ  (  ?5  ) 

1  .TANRZf  25) fR( 24> .RHO( 25) *TANBI (24) .TbFTA(24) .COSBI ( 24) »B(8) 

♦S(8»25)»U(8»25)»A(8*8)fF()6) fG#MMA(9) »GDTMB( 9) .ANS(9»5 ) *MC( 8) 

3  *1.2(70) 

COMMON  f ILLtPZ.LZ  tANS  . Ri * T ANRZ *R . RHO * T ANB I ♦ TBET A .COSB I  * R*S «U, A  * 
E*G*PHIZ .ALAM.RH.ZETA*TEMP.AMT.DELM.HDELM*Y.AI ♦SN1*SN2*TDEL 

2  ♦TrI tCRI * iRZfETA  *WN*nFT .MC*NST0P*MT.NPT*NG.NTM1 »M0PT 
equivalence  ( X*(*MMY*ANS) . (XTBI *ANS( 10) ) * ( XTB*ANS( 19) ) * (GAMMA*ANS( 
28 ) ) * ( GDTMB.ANS( 37 ) ) 

write  output  TAPE  2  *  lOO  *  NG  *X  (  4  )  *  ALAM  tMT  .NPT  »  ( MC  (  N  )  » N=  1  t.VPT  ) 

WRITE  OUTPUT  TAPE  2*i01 

write  output  tape  2*102*(A(N*1  )  .N»1*NPT) 

WRITE  OUTPUT  TAPE  2.103 

IF(MOPT)  1.1.2 

1  WRITE  OUTPUT  TAPE  2.104 

WRITE  OUTPUT  TAPE  2.105 

GO  TO  3 

2  WRITF  OUTPUT  TAPE  2*106 

WRITE  OUTPUT  TAPE  2.107 

3  DO  4  M-1.9 

k:»io-m 

WRITE  OUTPUT  TAPE  2 . 1 08 . ( ANS ( K  .  N )  .Nr  1 . 5 ) 

4  continue 
*  RETURN 

100  F0RMaT(25H1  NUMBER  OF  BLADES  '  G- 11 . 1 7H  LAMDA  I  AT  X«F4.2.4H 

IIS  F6.4/22H0  LATTICE  SPACES  MrI2,6H  I1.21H  CONTROL  POINTS 

2AT  Mr8l3) 

101  FORMaT(40H0  FOURIER  COEFFICIENTS  OF  G  A(l)  ) 

102  F0RMaT(5H0  4F10.6) 

103  FORMAT (25H0  G-GAMMA/TWO  PI  R  U»  ) 

104  FORMAT (2?M0  GBAR=6AMMA /TWO  PI  R  VA ) 

105  FORMAT (51U0  X  TAN  BETA  I  TAN  BETA  G  GBAR/) 

106  format (27H0  KAPP A=GOLDSTF I N  FACTOR) 

107  format (51H0  .  X  TAN  BETA  I  TaN  BETA  G  KAPPA/) 

108  format (F9.2 .F10.3.F1 1.3tF10.4.F12t4) 

END 


TABLE  B.4 


MAP  SUBROUTINE 


SUBROUTINE  MAP  (TEWP*RHt 
IF( TtMP-.999)  1»]»2 

2  TEMP=3. 1415926 

GO  TO  19 

1  CN=(1,+RH-2<t*TEMP)/(1.-RH) 

IFIABSFiCNi-.OOOOl )  17,17»18 

17  TEMP=l. 5707963 
GO  TO  19« 

18  CTN  =  SQRTF(  l.‘*CN**2  ) /CN 
TEMP=ATANF(CTN) 

IF(CTN)  20*19*19 

20  TEMP=TEMP+3,1415926 

19  RETURN 
END 
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TAftLS  B.5  (;4MBER  SUBROUTiN^ 


SU|iROUT?NE  CAMR&a 

fJT  MENS  ION  F  ILK  8000)  *A(  56*56  )  *AMS(  JA)  .R(56)  .ftUGl  8  )  ♦QHORDI  2A|» 

1  COSRI  (  24 1  »COFFZ<8  *7)  ♦D<  ?4)  *  E(  56  )  *  F  (*^6  )  tHMU  (  8  «7  )  .HI jf* 

2  PSlil6)*Pil6).PSIB(8)»RZ(25)*RI2A) .RHOI 24) .SNRHO(f) .TANRZI 25) . 

8  TANB! 1 24) .TRFTaIjA) »TIL (24) .TWETAI 8.7) .U ( 24 » 8 » 7  )  » WN I  16 )  .  X ( 9|  ♦ 

4  XCORD(9) .XTBI (9) .XTp(9) .XGAMIB) »XRH0(9) *MC(8) .NFLIPI 16J . 

5  LZ(69) .DUMMYI9.5) ♦GAMMA(24) 

COMMON  FILL  .PZ  .LZ  » A  » ANS  .  B  *BLIG.COEFZ  »C06Bi  . CHORD .O.E.F.rf. HMD 

♦  P.PSI  .PSIB.R.RHO.RZ.SNRHO.TANBl.TANBZ.fBAR.THETA.TIL.U.WN.XGAM* 

3  »  AtAM.AWT.AMT.AM.AI  .ANGLE.CBl  .DELM.DET.ETA  .G  »GNZL  .fiDELM  » PHI  Z  . 

4  RH.RBl .RB2.TBI .T RZ . TEMP .UB . W . Y . ZET A .MC . NFL  I P » J I N . JOUT . JT .KTEST . 

5  KlOl »MS.NBOTH.MT*NT.NPT.NZi .NZ2  »NG .NOT .NTT .NTMl  .N I P . NF . TBETA 

6  .GAMMA 

EOU I valence  I  X. FILL .dummy ) . I XcORD . P.DUMMY ( 10 ) ) . (XTBI .DUMMY ( 19) ) * 

1  (XTr. dummy (28) ) . (XRHO. DUMMY! 37) ) 

NIP=2*(NZ1-NZ2-1 ) 

NP  =  1 

DO  1  M=1*MT 

TIL  (M  )  =CH0RD(M)**C0SR  I  (  M  )  /  (  2.  *ANT*R  (  M  )  ) 

IF!M-MC(NP))  1.2.1  *  . 

2  DO  3  NQ=1.NQT 

TEMP  =  2**NZ1**N0-N1P  , 

THETAI NP.NO) =TIL(M)#TEMP' 

3  CONTINUE 
NP  =  NP,+  1 

1  CONTINUE 

K101=NPT*N0W 

DO  4  K=1.K101 

R(K )=0. 

DO  4  L=1.K101 

A(K.L)=0. 

4  CONTINUE 

DO  5  NU=1.NTM1.2 

TEMP=NU-NBOTH 
PSI  (NU)  =  TIL<  1)**TEMP 
PSI (NU+1)=PSI(NU) 

5  CONTINUE 
TBZ=TANBZ( 1 ) 

DO  6  N=1.NTT 

P(N)=PSI(N)-PSIINTT) 

6  CONTINUE 

CALL  LIST(NTM1*P.LZ) 

DO  38  NP=1.NPT 

MS=MC(NP) 

ETA  =  RZ(1  )/R(MS) 

TBI=TANBI IMS) 

CBI=COSBI(MS) 

DO  38  NO=l.NOT  *  *  ’ 

PHIZ=PSI INTT>-THETA(NP.NO) 

CALL  HEL IXIFTA.TRZ.TRI .CBI .PHIZ.NG.NTM1.LZ.WN) 

IF  (KTEST)  60.61.60 

60  WRITf  output  tape  JOUT.101.(WN(N) .Nsl .NTT) 

61  M=2 


© 
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DO  38  N=l*NBOTH 

tJ(N*NP»NO)=WN(i^> 

M  =  M+^, 

38  CONTINUE  Q 

*  -DO  M=1,MT 

IFIKTPST)  72 » 73 *72 

72  WRIIE  OUTPUT  TAPE  |JOUT.104*M 

73  RBl=RZtM) 

RB2=RZ(f<+l ) 

GAMMA(M) =0,0 

19  DO  21  I=1*NPT 

AI  =  I 

SNRHO(  I  )=SINF(AI*RHO(M) ) 

GAMMA( M) =GAMMA(M) +SNRHO( I )*XGAM( 1 ) 

«  21  CONTINUE 

N=1 

DO  22  NU=] fNTMl .2 

TEMP=NU-NBOTH 
PS  I  (NU)  =  TIL<M)»TEMP 
*  PSIB(N)=PSI(NU) 

IF(M-MT)  23*24.23 

'  23  PS  I (NU+1 )=TIL( M+1 )*TEMP 

GO  TO  25 

24  PSI (NU+1)=PSI (NU) 

25  IF(  PSI  (NU+1  )-PSnNU)  )  26*27*27  * 

26  NFLIP{NU)=0 
NFL IP<  NU  +  1 )  =  8 
AM=PSI (NU) 

PST(NU)=PSnNU+l) 

PSI (NU+1 )=AM 

GO  TO  84 

27  NFL1P(NU)=8 
NFL IP( NU+1 )=0 

84  N=N+1 

22  CONTINUE 

DO  8  NU=2*NTM1*2 

IF(PSI (NU+1)-PSI (NU)  )  9.8*8 

9  AM=PSI(NU) 

NF=NFLIP(NU)-1 
PSI (NU)=PSI (NU+l ) 

NFLIP(NU)=NFLIP(NU+1 )+l 
PSI (NU+1 )=AM 
NFLIP( NU+1 )=NF 

8  CONTINUE 

TBZ=TANBZ(M+1 ) 

DO  7  N=1*NTT 

P(N)=PSI(N)-PS|(NTT) 

7  CONTINUE 

CALL  LIST(NTM1 *P*LZ ) 

DO  14  NP=1*NPT 

Jl= (NP-1 )»NOT 
MS=MC( NP) 

6TA«l?Z(M  +  l)/fi(f1S) 

Taj=TANBI (MS) 

CBl'fOSBIIMS) 

IF(ZETA)  50*51*51 

91  ZfTA*>4ITAMnitM)-Tl^FTA(M)  ) / (  TfIT-rpPT 4(M&)  }  ) 

50  ALAM=R(MS)*Tni 
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IF(M-n  80*80*81 

80  BUG(NP)=?.*R(MS)*C0SRI(MS) 

81  00  14  NQ=i*NQT 

K=J1+N0 

if(M-MS)  83*82*83 

82  A(lf  »NP  )=12.^66575#GNZL*GAMMA(MS)»R(M5)#H(NQ)  /CHORD(MS) 
B(<)=B(K)+BUG(NP) 

B(K1  )=B(K1  )+B'JG(NP> 

8ft  F!K)=0. 

00  41  N=l*NeOTH 

U(N+16*NP*N0)=U(N.NP*NQ) 

41  CONTINUE 

PHiZ=PSI (NTT)-THETA(NP*NQ) 

CALL  HELlXIETA.TBZ.TBftCBI .PH I Z * NG * NTMl ♦LZ • WN ) 

IF  (KTEST)  62.63*62 

62  WRITE  OUTPUT  TAPE  JOUT . 1 0 1  * ( WN « N )  .  N= 1  *  NTT ) 

63  00  36  NU=1*NTT 

N=NFLIP(NU)+(NU+l)/2 
U(N.NP.NO)=:WN(NU) 

36  CONTINUE 

00  37  N=1.NB0TH 

ANGLE=PSIB( N )-THETA (NP.NQ ) 

CALL  BOUND ( RBI. RB2.ET A. AL AM. ANGLE* NG.UB) 
W=UB+U(N+8.NP.NQ)-U(N+16.NP.NQ) 

IF(KTEST»  64.65.64 

64  WRITE  OUTPUT  TAPE  JOUT . 1 02 *NP *  NO . N , UB * W 

65  IE(JT-?)  9‘*»40.40 

40  DO  35  I=1*NPT 

00  35  J=2*JT 

L=NPT+( I-l )« ( JT-1 )+J-l 

A(K *L)=A(K.L)-SNRH0( n*W*HMU(N*J)»AB5E(ZETAl 
35  CONTINUE 

99  F(K  )  =F(K  )+W*HMti{N#l  ) 

§7  CONTINUE 

B<K)=B(K.)+F(K)»GAMMA(M)*ABSF(ZETA)  ,  • 

14  CONTINUE 

write  OUTPUT  TAPE  JOUT. 107 

write  OUTPUT  TAPE  JOUT . 1 03  * ( ( A ( K . L ) . L  =  ! * Kl 01 ) . K  =  1 .K 1 0 1 ) 
WRITE  OUTPUT  TAPE  JOUT.108 

write  output  tape  JOUT. 103 .( n!K)  *k:=i  *K101  ) 

DET=3.0 

ME=XSIMEQF(56>K101 .1 .A.B.DET.E) 

, GO  10(68*69.69) .ME  » 

69  'call  ERR0R(20H  ERROR  IN  XSIMEOF) 

CALL  EXIT 
68  N  =  1  • 

DO  90  K=1*NPT 

MS=MC(K) 

ANS(N)=R(MS) 

ANS(N  +  1)=A(IC*1) 

ANS(N+2)=l./A(IC.l) 

N=N  +  3 

90  CONTINUE  ‘ 

J»3*NPT 

WRITE  OUTPUl  TAPE  JOUT. 109 

WRITE  OUTPUT  TAPE  JOUT .  1 1 0  . »  ANS  C  (^ )  ♦  N  =  I  *  J  ) 

K»NPT+1 
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DO  9i 

1=1 .tPT 

ir?  Ji-2) 

92.93.93 

93 

QO  94 

J=2. JT 

CO&Fli  I  1  J) 
K=K  +  1 

=  A(K.l  ) 

94 

CONTINUE 

92 

CONTINUE 

flUTPLfT  TAPE  JOUTilIf 

WRIT5  OUlPUT  TAPi  JOUT. 112. (  (COEF2:(*I  tJ)  *J=1*7)  .  1  =  1.8) 

RETURN 

101  F0RMAT(8H  WN=aF8.3) 

102  F0RMAT<5H  P=I3»5H  0=I3.5H  N=I3»^H  UB=F8.3.5)i  W=F8»3) 

104  FORMAT(15HO 

107  FORMaT(T3HO  coefficient  matrix  A(<.L)//) 

103  F0RMAT(8E15.5) 

108  F0RMAT{2B.H0  right  hand  S  I  OF  B(K)//) 

109  format (54H0  RADIUS  CAMBER  FACTOR  1C  CAMBER  FACTOR  1/K) 

no  FORMaT(7HO  F5.3».«H  F7.3.17H  F7.3) 

111  FORMaT(48H0  circulation  DISTRIBUTION  COEFFICIENTS  Cd.J)) 

112  FORMAT ( 7E15. 5) 

END 
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SUB  ROt*T  I  NS  BOUND  (  RR 1  ♦  RR?  ♦  F  T  A  ♦  AL  A'M  *  A  NGL  F  .  NG  ♦  UR  ) 
G=NG 

DEtBU=6.2831 853/G 

S»0. 

R=RB2/ETA 

RtAM  =  R/SOHlTF(R*-<f2+ALAM**2  ) 

A  =  R*#2  +  (  AUM#ANGLE  )  **2 

PHI=ANGee 

DO  1  N=1,NG 

T  =  0. 

CP  =  COSF(PHn 
SP=SINF(PHI) 

B=-2.>^R*CP  *  *• 

c=alam«*2*angle»cp+r-w-*?*sp 

X  =  RB1 

D  =  B^<-*2-4.#^ 

DO  2  1=1*2 

IF(  ARSF(  DI-. 00011  3t3»4. 

4  Y=-2,*(2.^fX  +  P)/(D*SORTP(  A+R»X4-X**2)  ) 

GO  TO  5 

Y=-l./(2.*(X  +  .5*B)*i*-2) 
rF(T-n  6»7»6  * 

T  =  T-Y 
X=R02 
GO  TO  2 

6  T=T+Y  *  ® 

2  CONTINUE 

S=S+T*C 
PHI=PHI+OELBL 
1  CONTINUE,*, 

UB=S*RLAM 

RETURN 

END 
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TABlE  B.7  HUMBUG  SUBROUTINE 


*  f^AP 

COUNT 

ENTRY 

BSS 

^UMpUG  SXD 
SXD 
SXD 
AXt 
CLA 
ADD 
STA 
T1  X 
CLA 
AOO 

,  StA 

AXT 
CL  A 
STO 
T  IX 
CLA 
STD 
AXT 

NUREV  SXD 
CLA 
PDX 
STA 
AXC 

NILE  SXD 
LXO 
TXI 
LDQ 
FMP 
STO 
LDQ 
FMP 
FSB 
.  STO 
•SXD 
.  AXT 

PINTO  LDO 
FMP 
-FAD 
STO 
STO 
XCA‘ 
FMP 
STO 
CLA 
5TC 
CLA 
\3X 
N’fi 
Pir 
BTi 


176 

HUMQI^'G 

3 

*-3*1 
*-l.»2 
*^3  »A 

56*! 

mI-8*1 

1  .4 
M+8  *  1 
*►3*1*1 
2*4 
ONE 
*  +  3 
48*1 
L+1  *1 
0.1 

*-2*1*1 
ONE  - 
XRIA 
8*2 

XR2A.2 

M+8.2 

0*1 

*+l  . 

0*2 

XRIB* 1 
XRIA* 1 
*+1 *1 *-l 
PHI  *1 

=.017453?P3 

X 

PHI+1.1 

=.017453293 

X 

DEL 

XRlA.l 

5*1 

DELT/+5*! 

DEL 

X 

X  . 

1*2 

X 

0.2 

=  1,0 

2*2 

X 

$C05*4 

*+? 

HUMPUC-I 


M 
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L 

M 


ONE 

PHT 


XCA 

FMP  X 

STO  5*2 

CLA  X 

1%X  SS1N«4 

NIR 

HUMBUG-1 
STO  4|2 

XCA 

FMP  X 

STd  692 

LDQ  GAIiSS+5»l 

FMI»  DEL 

STO  7f2 

TXI  *+l»2*-e 

TIX  PINTO*!?l 

LXD  XR1B»1 

II)?  NILE»1*1 

LXD  XR1A»1 

TXI 

SXD  XRlAfl  ^ 

LXD  XR2A»2 

Tf  X  NtJREV»2  f  1 

LXD  HUMBUG-3fl 

LXD  HUMBUG-2»2 

LXD  HUMBUG-1 #4 

TRA  3#4 

PZE  1720t0tl5 

P2E  1600»0»15 

PZE  Il480»0*15 

PZE  1360.0.15 

PZE  1240.0.15 

PZE  .  0.0.75 
BS5  15 

PZE  1720.0.15 

WE  1600.0,15 

PZE  1480.0.15 

PZE  1360.0,15 

PZE  1240.0.15 

PZE  600.0.50 

BSS  15 

PZE  ■  1^20.0,15 

PZE  1600.0.15 

PZE  1480,0.15 

PZE  1360.0.15 

PZE  1240.0.15 

PZE  1000.0.30  '« 

PZE  0.0.15 

PZE  600.0.10 

PZE  1000.0.6 

PZE  ,  1240.0.3 

PZE  1360.0.3 

PZE  1480.0.3 

PZE  1600.0.3 

PZE  1720.0.3 

PZe  1.0.1 

DPC  0.  .U  #2#  .4#  .7.  .If^#  .20..30,  .5ff.  .75«  f  lOO.  #150.  ♦  2^Tf)*  .?50# 


-120- 


> 


DE<? 

DEC 

OPfe 

Deg 

DELfA  DEg 
^AUoo  ngg, 
XRIA  PZE 
XRIR  PZE 
XRJA  PZE 
DEL  PZe 
X  pzr 
PND 


309.9360.  »(D.  010.  VIC.  « AO.  «60.  « 100. *1^0. » 200.  >27^. 

760.  *©.♦?, 0.  *50.  *90»*190.  »360.  .360.  .480.  .KOO.itiO. 

tiQ.  .fl40*.960..  lOftO.  .l9S0.  ♦  1320.  .I44e,.  1440..  1580» 

1680.  .1800,  .18(^0.  .1920.  .2040,  .2160. 
.046910..163865*.26923S*.269?3S».183855 
.116464.. 239314 #,284444. .23941 4.. 11 8464 


m 


Of 


s: 


sT 


lAftLE  B*8  *  LtSi  S<i<BROUTiN§ 


fAf> 

COtJNT 

176 

ENTRY 

LIST 

BSi 

3 

SXD 

*-3f  1 

sno 

^XD 

*-3#4 

CLA* 

1  *4 

STD 

NUH 

ft  A 

2*4 

STA 

AS  +  2 

APD 

=01 

6TA 

A5 

CLA 

3  9  4 

ADD 

=  01 

STA 

A9 

SUB 

=0* 

ST% 

*+l 

CIA 

** 

STO 

LAST 

STA 

^^  +  1 

4XC 

**94 

SXD 

A794 

AXT 

1  9  1 

SXD 

M91 

LXD 

NUH91 

CLA 

=06000000 

STD 

n 

SXD 

HUh 

LXD 

N92 

TXI 

*+19291 

SXD 

N92 

CLA 

**9l 

STQ 

X 

CLA 

**•1 

FSR 

X 

STO 

DFL 

LXD 

M9I 

FDP 

EPSLN+39r 

STQ 

D 

CLA 

D 

FSH 

=  •4 

TPL 

G5 

FAD 

=  •399 

TPL 

G2 

CLA 

=01000000 

STO 

e 

STO 

H 

CLA 

=010000000 

STD 

A2 

STD 

KL 

TRA 

46 

CLA 

=07000000 

-122< 


ST«) 

h2 

CLA 

^t)6QOOOO®, 

STD 

XL 

CLA 

=01000000 

S?D 

H 

CLA 

=92000000 

STD 

B 

TRA 

A6 

05 

CLA 

D 

UFA 

=0211001000000 

AffA 

=0000777000000 

5TD 

M 

LDO 

H 

MPY 

=(J5000000 

ALS 

IT 

STD 

B 

CLA 

=OSOOCOOO 

KL 

CLA 

=05900000 

STD 

A2 

CLA 

H 

LRS 

18 

ORA 

=0233000000000 

FAD 

=02^3000000000 

STO 

TEMP 

CLA 

DEL 

FDP 

TEMP 

STO 

DEL 

A6 

CLA 

M 

SUB 

=01000000 

XCA 

MPY 

=025000000 

ALS 

17 

ADD 

N 

PDX 

LDO 

LAST 

MPY 

=010000000  ' 

ARS 

1 

ADf) 

LAST 

STA 

LAST 

CtA 

B 

STD 

LAST 

CLA 

LAST 

A9 

STO 

STA 

AXC 

TXL 

ERROR»2t316 

A7 

TXH 

ERR0R»2 

*  LXO 

Htl 

A3 

SXO 

Ptl 

LXD 

KLfl 

Al 

LDO 

DELTA+1 #1 

FMP 

DEL 

FAD 

A  ^ 

ST® 

lf2 

TSX 

$C0S»4 

NTR 

*+2 

PZE 

LI5^T-1 

•1*3^ 

sra 

XCA 

FMP 

1tf2 

STO  ® 

5f2 

CLA 

lf2 

TSX 

$SIN*A 

NTR 

PZE 

.  LIST««1 

Si® 

XCA 

FMP 

1 

STO 

6f2 

LDO 

1*2 

FMP 

1  *2 

STO 

0*2 

CLA 

II 

• 

o 

STO 

2*2 

LDO 

GAUSS41 *1 

FMP 

DEI 

STO 

7*2 

TX  I 

*  +  l  *2  ♦-‘8 

TX  T 

*-4-1  *1*1 

A2 

TXL 

AI  *1 

CLA 

X 

FAD 

DEL 

STO 

X 

tXf> 

Ptl 

TTX 

A3*l*l 

LXD 

TIX 

A4*l *1 

LXD 

M*1 

TX  I 

*+1*1*1  ‘ 

TXL 

A8*l*3 

LXD 

LIST-3*1 

LXD 

L!ST-2*2 

LXD 

LIST-1*4 

TRA 

4*4 

ERROR 

TSX 

$MrST*4 

NU 

PZE 

®  *  NUH 

PZE 

N 

PZE 

M 

PZE 

*  LAST 

PZE 

*  A- 

X 

PiE  * 

DEL 

PZE 

* 

EP5LN 

OEC 

. 01745*. 087P7..3AQ07 

D 

PZE 

B 

PZE 

. 

H 

PZE 

KL 

PZE 

TEMP 

PZE 

DEC 

.  5  ».  866025*.  ?886  7  5*.  95^090*.  769;>35*’.  5*.  230765 

•* 

DELIA 

DEC 

.046910 

DEC 

CAtiSS  OEC 

.118464 

P 

PZE 

END 

• 

table  B.9 


HELJX  *Ui^R0UTIN6 


♦  FAP 

COUNT 

46f 

6NTRY 

HELIX 

HELIX  SXD 

HELIX-2f% 

SXA 

RESTOf 1 

SXA 

RggTO^l »2 

REM 

THIS  IS  THE  ^TART  OF 

LOO* 

5*4 

rMP'» 

114 

XCA 

FMP# 

2*4 

STO 

XZ 

CLA* 

6*4 

ARS 

18 

STA 

BLADS 

ORA 

=0233000000000 

FAD 

=0233000000000 

STO 

GFLO 

CLS* 

4*4 

STO 

CONST 

LDO* 

1  *4 

FMP# 

2  *4 

STO 

E% 

LDQ 

E4 

pMp 

E4 

STO 

FTP. 

LOO* 

1  *4 

FMP* 

1  *4 

STO 

El 

LOO# 

FMP# 

3*4 

FSB* 

1  *4 

FDP 

ET8 

FMP 

*-♦012665148 

fdf# 

2  #4 

FMP 

GFLO 

STO 

TRUNX 

CLA 

*01000000 

STO 

TEMP 

CLA 

*1.0 

FSB* 

1  *4 

SSP 

FSB 

=  .25 

IFL 

rouCh 

FAO 

=  .l^ 

TPL 

MED 

CLA 

TEMP 

ADO 

=02000000 

tl^A 

•  +  f 

MED  CLA 

Temp 

ADD 

*01000000 

STD 

TEMP 

ROUGH  CLA 

TEMP 

IHE  PARAM  fART  (GETS  CONSTANTS 


♦ 


♦ 


1 /A*PI 


125 


STD 

M 

AXC 

C8UG»2 

SXA 

Cf2 

AXQ 

D0UG»a 

SXA 

D*2 

axt 

Of  1 

CLA 

=  1.0 

STO 

J^LOOP  CLA 

GFLO 

F$B 

K 

FDP 

gflo 

FMP 

=6.283185 

FAD* 

5f4 

STO 

phik 

TSX 

SCOSfA 

NTR 

*+2 

PZE 

HELIX-2 

'STO 

CPK 

CLA 

PH  IK 

TSX 

$SINf4 

NTR 

*+2 

PZE 

HELIX^2 

STO 

SPK 

LXD 

HELIX-2f4 

LXA 

Cf2 

STZ 

•6  f  2 

STZ 

-5f  2 

LDQ* 

1  f4  • 

*  FHP 

CPK 

STO 

E2 

LDO* 

If4 

FMP 

SPK 

STO 

i3 

LDO 

•  E2 

FMP 

E4 

STO 

E5 

LDO 

E3 

FMP* 

E4 

STO 

E6 

*  LDO 

XZ 

FMP 

E2 

CHS 

FAD 

E6 

STO 

E7 

LDO 

E3 

fMP 

tz 

FAD 

E» 

STO 

E8 

LDO* 

3f4 

FMP 

E4 

CHS 

FAD 

El 

STO 

-4.2 

CDO* 

a<4 

PMP 

E% 

FSB 

E2 

STO 

-3*2 

LDO* 

m 


i* 


FMP 

UJ 

CHS 

FAD 

E3 

5T0 

-2*2 

LDO# 

3*4 

FMP 

E6 

STO 

-1  *2 

LDO# 

^  3*4 

FMP 

E5 

SiO 

0*2 

LXA 

D*2 

LOO 

ETB 

STO 

-6*2 

FMP<^ 

5*4 

XCA 

STO 

TE-MP 

FMP 

“2*0 

STO 

-5*2 

LOO 

TEMP 

FMP^^ 

5*4 

STO 

TEMP 

LDO* 

1  >4 

FMP* 

1  *4 

FAD 

=  1.0 

FAD 

TEMP 

STO 

-4.2 

LDO* 

1*4 

FfAP 

=  2.0 

XCA 

'*V 

STO 

TEMP 

FMP 

CPK 

CHS 

STO 

-3*2 

LDO 

TEKP 

FMP 

SPK 

STO 

-2*2 

STZ 

-1*2 

STZ 

0*2 

TXI 

*+1*2*7 

SXA 

D*2 

LXA 

C*2 

TXI 

*+1*2*7 

SXA 

C#2 

CLA 

K 

FAD 

=  1.0 

STO 

K 

TIX 

<LOOP*l.l 

REM 

START  HELTX  PART 

CLA* 

7*4 

STD 

NT 

CLA* 

6*4 

STD 

NO 

CLA 

ADRC 

STA 

A1 

CUA 

ADRD 

STA 

A2 

jCLA 

8*4 

PERFORMS  INTEGRATION 
NO  OF  ON  BLAOe  INTERVALS 
Flo  Of  PLADES 


9 


! 


ADO  *01  • 

£TA  A14 

5TA  /15 

CLA  9»4 

STA  All 

ADD  =01 

STA  All+1 

CLA  NT 

0  ARS  18 

SSM 

ADD  9*4 

STA  A13 

CLA  =0100000C 

STD  N 

STl  XNEW 

NUBLD  CLA  M 

STO  MBUG 

STZ  X 

SXD  NTBUG»4 

NUREV  CLA  MBUG 

SUB  =01000000 

XCA 

MPY  =O?50«C000 

ALS  IT 

AD®  N 

PDX  0,1 

A14  CLA 

5T0  LNM 

CLA  N 

PDX  0,1 

A15  CLA  ,  *#,1 

STO  LNl 

CLA  LNl 

ADD  .  =07 

STA  A3 

STA  AS 

SUB  =02 

STA  A4 

CLA  MBilG 

ARS  1 

ANA  =01000000 

SSM 

ADD  NG 

STD  NGBUG 

LXn  LN1,1 

CLA  A9 

STA  AlO 

A7  LXD  NGPUa»2 

SXD  XRBUG»2 

CLA  *  +  3 

STA  *+l 

Aft  AXT  0,2 

AXT  5,4 

STZ  T1 

STZ  T2 

LDO  **,2 

FMP  **,4 

FAD  T1 


L  +  l 


USED  10  CH&CK  CONVERGENCP 
INTEGRATION  SPACING  FACTOR 
SAVE  ORIGINAL,  H 


SELECT  DATA  TABLE 


21 (M»l )+N 

L+!  being  BACKWARDS  STORAGE 
ADDR  OF  P,0,N0  OF  POINTS 


t+1  SELECT  DATA  TABLE 

FOR  M=1  SPAflNG 

SET  UP 

ADDRESSES 

FOR  FIRST 

POINT  IN 

INTERVAL 

GET  NO  OF  BLADES 
IM  FIRST  GROUP 

MGBUG=NG  if  M«1 

ngbug=ng-i  if  'M  not  1 

POINTS  PFR  INTERVAL  COARSE 
TIX  A7,] ,1 

NO  OF  BLADES  IN  FIRST  GROUP 


7»K-2  FINDS  C  AND  D 
5  TERMS  FOR  ONE  POINT 
SUM  NUMERATOR  HERE 
SUM  DENOMINATOR  HERE 
C+7*NG 
P( N,M)+fi*J-l 


A1 

A3- 


STO 

ri 

A2 

LDO 

##*2 

•  D+7*NG-? 

A4 

FMP 

P(N«M)+B»J-3 

FAD 

T2 

STO 

T2 

TXI 

TIX 

A1  fA*! 

DOING  1  POI^IT  FOR  1  RL'ADE 

LDO 

T2 

DONE  1  POINT  FOR  1  BLAOg 

FMP 

T2 

XCA 

FMP 

T2 

TSX 

$S0RT»4 

NTR 

^  +  2 

PZE 

HELlX-2. 

STO 

U 

DEN0M**^/2 

CLA 

T1 

FDP 

T2 

A5 

FMP 

** 

P (N»M)+8»J-1*WE IGHT 

FAD 

X 

STO 

LXD 

’XRBUGf2 

TNX 

PTANOf2fl 

SXD 

XRRUGf? 

5ET  UP  for  same  point 

A8 

CLA 

A6 

NFXl  BLADE 

ADD 

=  07 

STA 

A6 

TRA 

A6 

PIANO 

CLA 

A3 

NEXT  POIMT  1ST  BLADE  6R0UP 

ADD 

=  010 

SIA 

A3 

SET  UP 

STA 

A5 

ADDRESSES 

SUB 

=  02 

FOR  NEXT  POINT 

STA 

A4 

IN  IfITFPVAL 

AlO 

TIX 

♦*♦1 ♦ 1 

A6  OR  A7 

CLA 

MBUG 

1ST  GROUP  DONE 

SUB 

=01000000 

TZf 

fJUINT 

IF  M  =  '!  all  PLAOES  have  BEEN  DONE 

CLA 

A8 

STA 

AlO 

TIX  A6.1.1 

CLA 

LNM 

ADD 

=  07 

SfT  UP  FOR 

STA 

A3 

MEDIUM  OR  FINE 

STA 

A5 

spacing  ON 

SUB 

=  0? 

Tndex  blade 

STA 

AA 

CLA 

A6 

Pirif  UP 

ADD 

=  07 

INDEX  BLADC 

STA 

A6 

CLA 

=oioooooc> 

riBUG 

Take  m 

LtO 

LNM*l 

NON-COARSE  spacing 

TRA 

AA  . 

aACK  TO  OO  i>Tr4.X  BLADE 

NUINO 

CLA 

n 

NfXT  iNTFliVAL 

SUP 

=07000000 

TMl 

^  +  2 

• 

PLADL 

CLA 

XNFW 

.129- 


STO 

XOLD 

LDO 

N 

MPY 

N 

ARS 

1 

ORA 

=0233000000000 

FAD 

=0233000000000 

STO 

Tl 

TEMP  STORAGE 

CLA 

TRUNK 

Tl 

XCA 

@ 

FAD 

X 

STO 

XNEW 

FSB 

XOLD 

SSP 

FSB 

=•0005 

allowable  truncation  error 

TM ! 

CONVR 

CLA 

A9 

STa 

AlO 

TIX  A7.1.1 

LXD 

N,1 

TXI 

*+l f If  1 

SXD 

Nf  1 

TXL 

NUREVflfS 

DO  MAX  6  (iEVS  DOWNSTREAM 

CONVR  LDO 

XNEW 

Pmp 

CONST 

A13  STO 

** 

WN  DOWNSTREAM  HELIX  DONE 

CLA 

NI 

NO  OF  INTERVALS  ON  BLADE 

TZE 

RESTO 

NO  INTEGRATION  ON  BLADE  RETURN 

CLA 

=07000000 

STD 

N 

DTD 

NTf4i 

SXD 

NTBUGfA 

TRA 

NUBLD 

BLADE  LXD 

flTRUGf  4 

Lt)0 

X 

FMP 

CONST 

All  FAD 

**f4 

WN 

STO 

♦  *•4 

WN+1 

LXf^ 

Nf  1 

TXI 

*  +  l f If  1 

SXD 

Iff  1 

TIX 

NUBLDf4f 1 

RESTO  AXT 

♦*f  1 

AXT 

♦♦#2 

LXD 

HELIX-2f4 

TRA 

10f4 

A9  PZE 

A7 

N  PZE 
NT  PZE 
NG  PZE 
NTBUG  PZE 
NGRUG  PZE 
M  PZE 
MBUG  Pit 
XRRUG  PZE 
LNM  PZE 
LNl  pze 
t  PZE 
D  PZE 


130. 


ADRC  PZE  CBUG+1 

AD?<0  PZE  DBUG-1 

trunk  PZE 
CONST  PZE 


SPK 

PZE 

CPK 

PZE 

K 

PZE 

PHIK 

PZE 

GFL(5 

PZE 

TEMP 

PZE 

etb 

PZE 

X 

PZE 

XOLO 

PZE 

XNEW 

PZE 

T1 

PZE 

Ti 

PZE 

BSS 

41 

CBUG 

PZE 

BSS 

41 

OBUG 

PZE 

El 

PZE 

E2 

PZE 

E3 

PZE 

64 

PZE 

E5 

PZE 

E6 

PZE 

E7 

PZE 

E8 

PZE 

XZ 

PZE 

END 

.131' 


UBU  B.io  •  iutho-uii  noasut  suau  oonor 


OFTZMM  Om-WXIR  FROnUBR 


NUHBER  OF  BLADES  G=3  LflHOfl  I  AT  X=0. 70  T?  0. 3333 


LATTICE  SPACES  M=24  4  COMTROL  POINTS  AT  M=  4  10  16  22 


FOWIER  COEPFICIEHTS  OF  G  A<I> 


0.1 396 VO  -0.008823  -0.000341  -0.000429 


G=GRMMA/TMO  PI  R  U* 


KPPPArGOLOSTEIN  FACTOR 


X  TAN  BETA  I  TON  BETA  G  KAPPA 

- 0720 - r76<57 - -TT. - -IT. - - 

0.30  1.111  -0.  0.0823*  0.8330 

- OTTO - (TTWJl - =ir. - 07TTT7 - 0;'8d70 

0.50  0,667  -0.  0.1320  0.8577 

O7G0  07556  =Xn - “TJ71'405 - O'.  8275^ 

0.70  0.476  -O.  0. I39I  0,7716 

0700  U.4I7  =07 - 077T32 - T5 . 6T7T 

0.90  0.370  -0.  0.1006  0.5146 

- rroo - 0,333 - =07; - O70TJO0 - 0.0000 

THE  DATE  IS  MAV  4.  1961. 

THETIAE  I^  1 31  9. 5 - 


KUaWUMFEED  FRORLLBR 


NUMBER  OK  BLADES  G=3  LAMOA  I  AT  X=D.7C  IS  0. 3332 
LATTICE  SPACES  M=24  4  CONTPOL  FCINTS  PT  M=  4  10  16  22 

FOURIER  COEFFICIENTS  OF  G  ACn 

0. 1 42769  -0. 009572  -0. 000906  0 . OOf  1  '9f 

G=GAMMA/TklO  PI  R  U* 


gbar=gamma/two  pi  R  VA 


X 

TAN  BCTA  I 

T«N  BETH 

G 

OF  ftp 

0.20 

1.415 

0.910 

0. 

0. 

0.30 

1.006 

0.691 

0.  03*42 

0.  C3!?4 

0.40 

0.737 

0.563 

0. 1 153 

0.C45:.' 

0.50 

0.646 

0.475 

0. 13h2 

0.1:405 

0.60 

0.543 

0*41  g 

0*143? 

o.ro 

f>.i61 

t.€4^r? 

o.eo 

0.421 

0.322 

0. 1 .320 

0.‘  406 

0.90 

0.377 

0.290 

0.  10S2* 

0-*  30'» 

t.oo 

0.342 

0.265 

O.OC'OO 

0 .  C  C  Ou 

THE  DRTE 

IS  MP,V  4, 

1961. 

THE  time  IS  2341.1? 


suu  >.u  -  umm  i 


NAME  ORIGIN 

entry 

NAME 

ORIGIN 

SUBPROGRAM  STORAGE 
ENTRY  NAME  ORIGIN 

MAP 

ENTRY 

NAME  < 

origin 

ENTRY 

NAME 

ORIGIN 

ENTRY 

(MAIN) 

00144 

00163 

MAP 

01374 

01402 

OCTAl.S 

01522 

01531 

HUMBUG 

02130 

02137 

LIST 

02454 

02464 

MIST 

02764 

02773 

HELIX 

03032 

03037 

IN T ESP 

03754 

0375  7 

BOUND 

04221 

04230 

AKL 

04573 

04610 

(FPT) 

07013 

07021 

EXIT 

07435 

07463 

EXITM 

07435 

07441 

(T5HM) 

07475 

07520 

(TSH) 

07473 

07531 

( CSH  > 

07475 

07511 

(F2EF) 

07761 

10156 

FTNPM 

07761 

10023 

(F2PM) 

07761 

10017 

(PCUP) 

14236 

14341 

ISTPC) 

14236 

14240 

(PRLTI 

14616 

15032 

(RSLT) 

14616 

14753 

(SVLT) 

14616 

14633 

(BCD)) 

16005 

16010 

ICMPR) 

16027 

16031 

(DCDR) 

16114 

16116 

(FLO)  ) 

16175 

16204 

(FIX)  ) 

16175 

16202 

(  ILSC) 

16262 

16264 

(MOVE) 

16324 

16326 

(NBLX) 

16404 

16407 

(OCT)  ) 

16467 

16471 

(OCTO) 

16520 

16522 

tOPCO) 

16575 

16577 

(PRNT) 

17153 

17156 

(PSTN) 

1735.3 

1735.7 

V74J4 

11556 

ISEtX) 

174  74 

17477 

(SPOT) 

17701 

1  7703 

STOMAP 

17732 

17737 

(CSHM) 

20161 

20164 

(SPH) 

20165 

2020S 

(STHM) 

20165 

20211 

(STH) 

20165 

20236 

ISPHMl 

20636 

20641 

(WTCI 

20642 

20705 

(WER) 

20642 

20654 

( f;DC ) 

20733 

20764 

(RER) 

20733 

20742 

IRTN) 

21000 

22450 

(FID 

21000 

22437 

(  lOH) 

21000 

21002 

(TCO) 

22601 

22675 

(TEF) 

22601 

22674 

(RCH) 

22601 

22673 

(ETT) 

22601 

22672 

(REW) 

22601 

22671 

(WEF) 

22601 

22670 

(SSR) 

22601 

22667 

(WRSI 

22601 

22666 

(RDSI 

22601 

22665 

(  lOS) 

22601 

22606 

(TRC) 

22601 

22676 

(  lOU) 

22732 

22735 

ATN 

22750 

22752 

ATAN 

22750 

22752 

SIN 

i3Q/i7 

.210.72 

CQ5 

23067 

,  21071 

SQR 

^3242 

23246 

SORT 

23242 

23246 

(TESJ 

23337 

23337 

( EXE ) 

23340 

23344 

RECOUP 

24317 

24322 

ERROR 

24325 

24331 

LOUMP 

24501 

24504 

TIME 

24510 

24562 

CLOCK 

24510 

24515 

GETTM 

24644 

24646 

XDETRM 

25005 

25526 

XSIMEQ 

25005 

25376 

MOV  I E ) 

2  562  7 

25627 

THE  DATE  IS  APRIL  17» 
THE  TIME  IS  1508.4 


NT«3 

MT»  8  NPT 

*3  N21»l 

NZ2-0  MC 

-  3  5 

7  0 

0  0  0 

0  N6»3 

ALAH-0.2423 

CH0RD= 

0.233404 

0.259304 

0,282426 

0.300560 

0.311472 

0.313185 

0.298909 

0.22821S 

TANBI- 

0.969200 

0,692286 

0.538444 

0.440545 

0.372769 

0.32*3067 

0.285059 

0.255053 

TBETA« 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

GAMMA« 

0.114469 

0.126195 

0.106211 

HMU  » 

0.109375 

0.182292 

0.208333 

0.189887 

-0.002532 

-0.107355 

H 

-1.333333 

-2.666666 

COEFFICIENT  MATRIX  A(K.L) 


-0,24452E“00 

0. 

O.17507E  01 

0. 

-0.80969E  00 


0. 

-O.SieSZE  01 
0.84853E  01 

0. 

-0.86532E  01 


-0.41693E  01 

0. 

-0«404e4E-00 
0.12562E  02 


-0.52090E  01 
0.26334E  01 
-O.7O60OE  00 
-0.85590E  01 
-0.77153E  01 


-0.43664E  01 

0. 

0. 

0.11601E  02 


0.25272E  01 
-0.35300E-00 
-0.74031F  01 
-0.57495E  01 


-0.48904E-00 

0. 

0.2?1??F  01 

0. 


“0.743??E  01 
0.83? 95F  01 

0. 


RIGHT  HAND  SIDE  B(K} 
“0.27468E-00  -0.56039E  00 


RADIUS 

0.450 

0.650 

0.850 


CAMBER  FACTOR  K 
1.169 
1.295 
1.682 


‘0.44718E-00  -0.905n5E  00 

CAMBER  FACTOR  1/K 
0.856 
0.  //2 
0.594 


-0.6994'9E  no  -0.13839E  01 


CIRCULATION  DISTRIBUTION  COEFFICIENTS  C(I»J) 


0,12837E-00 

-0,44394E-03 

0. 

0. 

0. 

0. 

0. 

-0,49850E-02 

-0.13892E-02 

0. 

0. 

0. 

0* 

6, 

0,25870E-02 

0.10633E-02 

0. 

0. 

0* 

0. 

0. 

-0, 

0. 

0* 

0. 

0« 

0. 

0* 

-0, 

0. 

0. 

0. 

0. 

0. 

0. 

“O. 

0. 

0. 

0. 

0. 

0. 

0. 

-0. 

0. 

0. 

0. 

0. 

0. 

0. 

-0. 

0. 

0. 

0. 

0. 

0. 

0. 

THE 

THE 


date  is  /PRIL  17. 
TIME  IS  1511.3 
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APFBHDIX  C 


TkBlM  OF  CEORD>LQAD  FACTORS 


n  J-l  J  =  2  J  =  3 


N«:2 

in 

2 

.5  ! 

.5  : 

1 

1 

Fe4 

1 

2 

3  i 

4 

.195312 
.304888  ‘ 
.304688 

.195312 

.292969 

.152344 

-.152344 

-.292969 

1 

.109375 

.182292 

.189887 

2 

.182292' 

.182292 

-.002532 

V-6 

3 

.208333 

.069444 

-.187355 

4 

.208333 

-.069444 

-.187355 

5 

.182292 

-.182292 

-.002532 

6 

.109375 

-.182292 

.189^7 

1 

.126010 

.146058 

2 

.123367 

.154209 

.068072 

3 

.147079 

.110310 

-.065429 

M»8 

4 

.157547 

•039387 

-.148701 

5 

.157547 

-.039387 

-.148701 

6 

.147079 

-.110310 

-.065429 

7 

.123367 

-.154^ 

.06^2 

8 

.072007 

-.126010 

.146058 

OEFmD  IM  ($.20) 


-4  J-5  J-6  J-7 


I 

I 


.134187 

-.184823 

-.131896 

.131896 

.184823 

-.134187 


.129591 

.080010 

.010426 

-.052692 

.069119 

-.154856 

-.124451 

.004068 

.153961 

-.054423 

.118564 

.147254 

.076562 

.129269 

.108879 

-.098430 

.076562 

.12^9 

-.108879 

-.0^30 

.153981 

-.154423 

-.118564 

.147254 

.069119 

-.154856 

+.124451 

.004068 

.129591 

-.129591 

.080010 

-.052892 
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